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ABSTRACT 

Condensation heat transfer in a horizontal rectangular duct was experimentally 
and analytically investigated. To prevent the dripping of condensate on the film, the 
experiment was conducted inside a horizontal rectangular duct with vapor condensing 
only on the bottom cooled plate of the duct. R- 113 and FC-72 (Fluorinert Electronic 
Fluid developed by the 3M Company) were used as the condensing fluids. 

The experimental program included measurements of film thickness, local and 
average heat transfer coefficients, wave length, wave speed, and a study of wave 
initiation. The measured film thickness was used to obtain the local heat transfer 
coefficient. The wave initiation was studied both with condensation and with an 
adiabatic air-liquid flow. The test sections used in both experiments were identical. 

Experimental results showed that the average heat transfer coefficient increased 
with increasing inlet vapor velocity. There was a significant increase in the heat 
transfer after the appearance of interfacial waves (when the inlet Reynolds number 
was approximately greater than 1,000,000). The local heat transfer coefficient 
decreased with axial distance of the condensing surface. The local heat transfer 
coefficient decreased rapidly with axial distance near the leading edge of the 
condensing surface but was nearly constant towards the trailing edge. 

It was observed that the condensate flow along the condensing surface 
experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional 
wavy flow. The change in the flow pattern depended on the vapor velocity and the 
difference between the saturation temperature of the vapor and the condensing surface 
temperature. The wave length decreased with axial distance and the inlet vapor 
velocity, while the wave speed increased with vapor velocity. 

An analytical model simulating the condensation process was formulated by 


employing the universal conservation laws. The resulting equations were solved 
numerically. The heat transfer coefficients predicted from the model are within +20 
% of the measured values. 

Using the present experimental data and analysis, correlations for the average 
heat transfer coefficients in the annular flow regime were developed. The average 
deviation between the predictions and the experimental values is within ±18%. 

Also, predictions from some correlations selected from the literature for the heat 
transfer coefficients of condensation in the annular flow regime were compared with 
the experimental data. It is found that the correlations of Shah (1979) and Soliman et 
al. (1968) yielded satisfactory predictions. (The average is within ±21 % for Shah’s 
correlation and ±22 % for Soliman’s correlation.) 

The effects of air velocity, liquid flow rate, and the liquid viscosity on the 
initiation of interfacial waves were studied with an adiabatic air-liquid flow. The 
liquid viscosity was varied by varying the mass fraction of glycerine in the water- 
glycerine mixture. It was observed that the interfacial waves were initiated closer to 
the leading edge when the air velocity or the liquid flow rate was increased. 

It was found that with condensation, the flow appeared more stable in the 
adiabatic liquid flow; condensation appears to have a damping effect on the initiation 
of interfacial waves. 

Several criteria for the instability in two-phase flow were examined and compared 
with each other. It is found that none of them can be confidently used in the present 
experiment. Wit^ the modification of one of the investigated criteria and the use of 
the present experimental data of the air-liquid flow, a criterion for predicting the 
instability of the two phase flow is tentatively given. This criterion provides a 
reasonable prediction for the initiation of interfacial waves with air-liquid flow, as 
well as with condensation. 
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1. INTRODUCTION 


In applications involving condensation of a vapor, the condensate is drained by 
the gravitational forces in external condensation or by the shear force if condensation 
is inside horizontal tubes. In space applications under microgravity conditions, for 
draining the condensate mechanisms that do not depend on gravity are needed. One 
possibility is condensation inside tubes. All the characteristics of condensation inside 
horizontal tubes are not fully understood, and those that are understood under 1-g 
conditions may occur differently under microgravity conditions. 

With condensation inside horizontal tubes it has been established that the major 
part of the heat transfer occurs in the annular and annular-wavy regimes. The results 
of an experimental and analytical study of condensation in the two regimes are 
presented in this thesis. The results can be used as a base for comparison with results 
of condensation under reduced gravity conditions. 

With condensation inside horizontal tubes interfacial instabilities lead to the 
annular-wavy regime. Amplification of the waves may lead to undesirable flow 
instabilities. Thus, there is a need to study the role of various parameters that affect 
the stability of the condensate film, and the effect of the interfacial waves on the heat 
transfer rate. 

In the annular flow regime in horizontal tubes under 1-g conditions there is the 
possibility of the condensate from the upper half of the tube dripping on the lower 
half leading to waves in the condensate. To avoid the effect of such dripping 
condensation, a rectangular horizontal duct with vapor condensing on the bottom 
cooled surface of the horizontal duct was chosen for the experimental and analytical 
study. This dissertation presents the results of condensation heat transfer in the 
annular and annular-wavy flow regimes with the condensate drained by the vapor 
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shear, and predictions from a simple analytical model. Refrigerant 113 (R- 113) and a 
fluorinert electronic fluid (FC-72 developed by the 3M company) were selected as the 
condensing fluids. The experiments yielded the condensate film thickness at several 
axial locations, the local and average heat transfer rates, and the wave length and 
wave speed of interfacial waves when they appeared. Observed conditions for the 
incipience of interfacial waves are also reported. 

To better identify the parameters relevant to the initiation of the interfacial waves, 
and to relate the effect of condensation on the initiation of such waves another series 
of experiments with an adiabatic air-liquid flow were performed. These experiments 
were conducted in a test section which was identical to the test section used in the 
experiments. 

The condensation was simulated by an analytical model. The resulting equations 
were solved numerically yielding the condensate film thickness, The interfacial shear 
stress, and the local and average heat transfer coefficients. The validity of the model 
is examined by comparing the predictions with the experimental results. There is 
considerable differences between the predictions and the experimental results. The 
difference is attributed to the lack of a proven model for the interfacial shear stress. 

The model prediction are good at moderate inlet velocities of the vapor but differ 
significantly from the experimental results at low and high inlet velocities. 
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2. LITERATURE REVIEW 


2.1. Introduction 

Condensation heat transfer (with condensation inside ducts) is influenced by the 
flow regimes and usually by the orientations of the condensing surface, if the vapor 
shear stress is not high enough for the gravity force to be negligible. Fig. 2.1 
illustrates the flow regimes occurring during condensation inside a horizontal tube. 

As the condensation proceeds along the tube, the liquid flow rate increases and the 
vapor flow rate decreases; the flow regime changes along the tube. Unless the vapor 
velocity is very low, the flow begins as an annular flow and is followed by wavy, 
slug, and plug flows (or stratified flow when the inlet vapor mass flow rate is low). 
The annular flow regime is known to be the dominant flow regime, existing over 
most of the condensing length both inside a horizontal tube and a vertical tube 
[Soliman et al, (1968)]. In the annular flow regime, condensate film forms on the 
periphery of the tube; the flow of the liquid film is driven mainly by the momentum 
and the shear stress of the vapor flow in the core of the tube. With increasing vapor 
velocity, waves appear on the surface of the liquid film. The annular flow regime 
with interfacial waves is termed annular wavy flow regime. 

In this review, emphasis is on a survey of the the analytical and experimental 
studies of condensation inside a horizontal tube. A brief review of condensation of 
vapor flowing parallel to a cooled surface is also included as there are some 
similarities between the two. 

The annular flow regime also includes the annular wavy flow regime. As the 
interfacial waves have a significant effect on the heat transfer rate, a survey of the 
literature of the interfacial instability is also given. 
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Fig. 2.1 Illustration of flow regimes of condensation in horizontal tubes 

[This figure is reproduced from Sardesai et al. (1981)1 

2.2. Condensation 

2.2.1. Analytical Studies of Condensation over a Flat Plate 

Nusselt (1916) first analyzed the film condensation of a quiescent vapor on an 
inclined plate. In his analysis, he assumed the following: the condensate film is 
laminar; all the properties of vapor and liquid are constant; the subcooling, 
momentum change, and interfacial shear stress are negligible; the temperature profile 
in the condensate is linear. 

Jakob (1949) was among the first to point out a high deviation of the Nusselt 
analysis with the experimental results involving a turbulent condensate film flow. 
Akers and Rosson (1960) also showed that the Nusselt analysis may not be applicable 
when the effects of the liquid subcooling or the vapor shear are significant or when 
the liquid film is turbulent. 

For condensation of vapor flow over a flat plate, Koh (1962) presented numerical 
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solution to a system of differential equations of mass, momentum, and energy 
balance. The equations included inertia and convective terms. The results of his 
analysis showed that for a thin condensate layer and for Prandtl number less than 10. 
the condensate velocity and temperature profiles are essentially linear. Koh showed 
that for higher Prandtl numbers, neglecting the convective terms in the energy 
equation considerably underestimated the heat transfer rates. 

Shekrilasze and Gomelauri (1966) expressed some skepticism about Koh’s result 
regarding the effect of the inertia forces and the convective terms on the condensation 
heat transfer rate. They developed a model of condensation on a flat plate, based on 
the assumptions that the temperature profile is linear and the inertia forces are 
negligible. They assumed that at high vapor velocities, the interfacial shear stress is 
equal to the momentum associated with the condensing vapor. Approximate solutions 
for predicting the heat transfer coefficients were given for condensation both with a 
uniform surface temperature and a constant heat flux. From the results of the model, 
they found that the assumption of negligible inertia forces is appropriate for the 
condensation of non-metallic liquids and that the subcooling of the liquid was 
insignificant for N < 1 (N=k,AT/h lg /x 1 ). They also reported a good agreement when 
comparing their model predictions with the experimental data of Jakob (1935). 

Sparrow et al. (1967) studied forced convection condensation on a horizontal flat 
plate. Their model included the effect of a non-condensable gas, but neglected the 
inertia and the convective terms. They used similarity analysis and numerical 
methods to solve the equations. They showed that in forced convection condensation, 
the effect of non-condensable gases was insignificant. They also indicated that the 
interfacial thermal resistance due to the effect of the transport of the molecules from 
the condensed vapor and the molecules evaporating from the surface was negligible 
for the forced convection condensation of steam. 

Uehara et al. (1984) analyzed the turbulent film condensation of a saturated vapor 
in forced flow over a flat plate. They assumed a turbulent liquid layer and a 
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turbulent vapor boundary layer above the liquid. The temperature and the velocity 

profiles in the liquid layer and the vapor boundary layer were assumed. In the 

analysis, they found that the thickness of the sublayer within the vapor boundary layer * 

influenced the heat transfer, and, therefore, was included in the equation for the heat 

9 

transfer coefficients. For their predictions to completely agree with experimentally 
determined correlations for heat transfer coefficients, the value of the sublayer 
thickness was adjusted. 

Narain and Kizilyalli (1991) developed a set of scaled mass, momentum, and 
energy differential equations to model the condensation of pure saturated vapor flow 
between two horizontal, parallel plates. In their study, the vapor flow was considered 
to be a pressure driven laminar flow. It was also assumed that condensation occurred 
only on the cooled bottom plate with the upper plate adiabatic, and that the 
temperature profile of the condensate was linear. Narain (1991) further simplified the 
model by dropping the gravity, inertia, and convection terms in the liquid phase. 

Comparison of the predicted heat transfer coefficients with the present experimental 
data showed agreement within +40 % . Using the data from the present experiments, 

Narain and Kamath (1991) refined the friction factor )del of the interface for both 
laminar and turbulent flows of vapor. The resulting agreement between the theory 
and some of the experiments with regard to film thickness was within ±10 %. 

2.2.2. Analytical Studies of Condensation inside Tubes 
Low Mass Flow Rate of Vapor 

For condensation of vapor inside a horizontal tube at low inlet mass flow rates, 
condensation occurs on the inside wall of the tube and the condensate drains to the 
bottom of the tube in the same vertical plane, similar to the external condensation of 
quiescent vapor on a horizontal cylinder. The condensate at the bottom of the tube is 
a thick liquid layer, while the condensate on the rest of the wall is a thin liquid film 
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(Fig. 2.2) . 

Chaddock (1957) studied 

* condensation of vapor in a 

S 

horizontal duct at low vapor 
t velocities. He postulated 

that the condensate was 
driven by the hydraulic 
gradient of the condensate 
(Fig. 2.2 a). The surface 
profile of the bottom 
condensate was estimated on 
this assumption, and defined 
the area covered by the 
condensate film. It was further assumed that condensation occurred only on the 
surface above the bottom condensate layer and that the condensation heat transfer rate 
was given by the Nusselt analysis on a horizontal cylinder. 

Chato (1957) adopted the Chaddock model to estimate the surface profile of the 
bottom condensate in a horizontal tube. The heat transfer to the bottom condensate 
was neglected, while the heat transfer to the condensate film on the upper part of the 
inside tube wall was derived from an analytical model involving the integral 
momentum and energy equations. The predicted heat transfer coefficient agreed with 
his experimental data when the vapor Reynolds number was less than 35,000. 
Comparing his analytical model with the Nusselt analysis, Chato concluded that the 
Nusselt analysis was applicable when the liquid Prandtl number is of order 1 or 

* 

greater, but when the liquid Prandtl number was less than 1 , using the Nusselt 
* analysis would underestimate the heat transfer coefficients. 

Rufer and Kezios (1966) regarded the condensate flow model assumed by 
Chaddock (1957) and Chato (1961) as unrealistic. They assumed that the condensate 



(b) Shear driven flow model (Rufer and Kezioa) 


Fig. 2.2 Condensate flow models 
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depth at the bottom of an inclined or a horizontal circular tube increased in the 
direction of vapor flow (Fig. 2.2 b). They also assumed that the flow of the bottom 
condensate was mainly by the pressure gradient. An equation predicting the depth of 
the bottom condensate was derived from the mass, momentum, and energy balance 
equations of the bottom condensate. Rufer and Kezios also included the heat transfer 
to the bottom condensate. The average heat transfer coefficient was evaluated from 
the Nusselt analysis. In the momentum equation, they neglected the interfacial shear 
stress between the vapor flow and the bottom condensate on the basis of the 
assumption that the bottom flow of the condensate was stratified. The predicted 
condensation depth at the bottom was compared with the experimental data in terms 
of the flow angle of the bottom condensate wetted on the inside tube surface. They 
found that the predicted value was in good agreement with the experimental data. 

Comparing the equation of Rufer and Kezios (1966) for the condensate depth with 
that of Chato (1962), Khabenskiy et al. (1981) concluded that both equations were 
basically the same. However, those two equations would give contradicting results 
due to the different initial conditions assumed for each equation (Chato assumed that 
the condensate depth increased along the condensing tube, while Rufer and Kezios 
assumed that the condensate depth decreased with increasing the length of the tube). 
Using an example for a condensing system including an horizontal or inclined 
condensing tube, Khabenskiy et al. indicated that the condensate flow patterns 
assumed by both Chato (1962) and Rufer and Kezios (1966) existed in practice. They 
considered that the flow model of Chato described the case of an incomplete 
condensation of vapor inside a horizontal or an inclined tube, while that of Rufer and 
Kezios simulated the situation of complete condensation inside the tube. 

Rosson and Myer (1965) found that their experimentally determined heat transfer 
rates across the bottom condensate inside a horizontal condensing tube (using 
methanol and acetone as condensing fluids) was about 23 % of the heat transfer rates 
across condensate film covering the upper part of the inside surface of the tube. They 
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therefore concluded that the heat transfer to the bottom condensate should not be 
neglected. The heat transfer to the bottom condensate was evaluated using the von 
• Karman analogy between heat and momentum transfers. The heat transfer to the 

condensate film was determined using the Nusselt analysis. To include the effect of 
| the interfacial shear stress in the Nusselt analysis, they used an experimentally 

determined function of vapor Reynolds number to replace an original constant in the 
Nusselt equation. A comparison of the model predictions with the experimental data 
showed that the predicted values were on the low side if the liquid Reynolds number 
was low and were on the high side if the liquid Reynolds number was high. They 
assumed that the deviation was caused by the use of the Lockhart-Martinelli (1949) 
correlation, which might not have been suitable for condensation. 

Moalem and Sideman (1976) analyzed the simultaneous process of condensation 
of a vapor flow inside a horizontal tube and the evaporation of a fluid on the outer 
surface of the tube. Neglecting the effect of the interfacial shear stress on the 
condensate film, which was considered to be appropriate for vapor velocities of less 
than 15.2 m/s, they solved the momentum and energy equations using integral 
methods. They predicted the variation of the local heat transfer coefficient with the 
periphery of the tube. Similar to the experimental results of Rosson and Myer 
(1965), they showed that the heat transfer coefficient at the bottom of the tube was 
about 20 % of the highest heat transfer coefficient. Comparison of the predicted 
overall average heat transfer coefficient of condensate and evaporation with one 
source of experimental data showed that the predicted values were generally 30 - 40 
% lower than the experimental data. The predicted low heat transfer coefficients 
might have been due to the neglect of the interfacial shear stress in the analysis. 

V 

' High Mass Flow Rate of Vapor 

High mass flow rate of vapor at the inlet to a horizontal tube leads to a thin 
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condensate film around the periphery of the tube. The condensate Film flows in the 
axial direction of the tube, resulting in the annular condensation regime. 

Soliman et al. (1968) developed an annular condensation model by modifying the t 

energy equation (for condensation inside a tube) given by Carpenter and Colburn 
(1951). The main modification was in the expressions for the friction force including f 

interfacial shear stress and the momentum change of condensate and vapor. In 
calculating the friction force, Soliman et al. neglected the effect of condensation of 
vapor (on the condensate) in the total interfacial shear stress and determined the 
friction force using the Lockhart and Martinelli method applied to an adiabatic two- 
phase flow. In evaluating the momentum change, they considered the condensate film 
as a turbulent liquid layer and assumed the ratio of the interfacial velocity to the 
average condensate velocity to be 1.25. By substituting the equations for the friction 
force and the momentum change into the Carpenter and Colburn expression, an 
equation predicting the heat transfer coefficient in an annular condensation was given. 

Two constants in the equation were determined based on a regression analysis with 
nine different sources of experimental data for condensation in both horizontal and 
vertical tubes. However, no comparison of the equation with other experimental data 
that were not used for determining the constants in the equation was reported. 

The model of Soliman et al. (1968) was modified by Bae et al. (1971). In the 
modification, the von Karman momentum-heat transfer analogy was used to determine 
the heat transfer coefficient, where the ratio of the eddy conductivity to the eddy 
viscosity was assumed to be 1 and the von Karman universal velocity distribution was 
used to determine the condensate velocity. Also, the ratio of the interfacial velocity 
to the average liquid velocity, which varied approximately between 2 and 1.1, was 
assumed to be a unique function of 5 + (5 + =5/i>|J g rjp). With increasing 5 + , the 
velocity ratio decreased. From the analysis, an expression for predicting the heat * 

transfer coefficient was derived, but the calculation of the heat transfer coefficient 
required considerable iteration. Subsequently, Traviss et al. (1973) simplified the 
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expression for heat transfer coefficients by neglecting a term representing the 
correction due to the difference between the wall and interfacial shear stresses. They 
claimed that neglecting the term was valid when the gradient of the quality of the 
mixture was not large. The predicted heat transfer coefficients were compared with 
their experimental data and the data of Bae et al. (1969). They found that the 
predicted values were generally lower than the experimental data. They considered 
that the deviation resulted from model not taking into account the effect of the 
entrainment of the condensing liquid; some of the experimental data was taken when 
entrainment occurred, which reduced the thickness of the condensate and increased 
the heat transfer coefficient. 

Analyzing a heat transfer model similar to the one given by Bae et al. (1971), 
Cavallini and Zecchin (1974) determined the major non-dimensional parameters 
affecting the average Nusselt number in an annular condensation. Based on the 
regression analysis of their experimental data, an equation for the Nusselt number is 
given. 

Jaster and Kosky (1976) studied condensation heat transfer in a mixed flow 
regime between an annular flow and a stratified flow (Fig. 2-1 a). First, they derived 
a criterion in terms of the ratio of the axial shear force and the gravitational body 
force to determine the annular and the stratified flow regimes. Based on this 
criterion, a heat transfer correlation for the mixed flow regime was given simply by 
using a linear interpolation from the known heat transfer correlation specified for an 
annular flow and from the available heat transfer correlation for a stratified flow. 
Comparison of the calculated heat transfer coefficient with the experimental data in 
the mixed regime showed an average deviation of ±34 %. 

Shah (1979) developed a correlation for condensation inside tubes. This 
correlation was obtained from a modification of the single phase heat transfer 
correlation of Dittus and Boelter (1930) for a turbulent flow of a single phase fluid. 
The Dittus and Boelter correlation was multiplied by a term involving the quality. 
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The modification was derived by an examination of correlations for film boiling based 
on the assumption of similarity between the mechanisms of heat transfer during film 
condensation and film boiling inside tubes. Shah compared his heat transfer 
correlation with 474 different sets of experimental data from 21 independent sources 
for condensation inside both horizontal tubes and vertical tubes. The mean deviation 
is reported to be within +15.4 %. 

Kutsuna et al (1985) derived a system of ordinary differential equations for the 
film condensation inside a horizontal rectangular duct. In the derivation, they 
assumed that the temperature profile in the condensate film was linear and the 
interfacial shear stress was calculated using the Shekriladze approximation (j x ~ 
rivVg). The predicted heat transfer coefficient was in good agreement with their own 
experimental data. 

2.2.3. Studies of Interfacial Shear Stress 

The interfacial sh_ar stress in condensation plays an important role in influencing 
the heat transfer rates. However, because of the difficulty in measuring the interfacial 
shear stress, different models have been proposed, but none has been validated. 

Rohsenow et al. (1956) used a force balance on an element of condensate flowing 
along a vertical condensing surface to derive the interfacial shear stress. In their 
derivation, the momentum changes through the element (including the momentum 
effect due to the condensation of the vapor were assumed to be small and negligible. 

A number of researchers [for example, Rosson and Myer (1965), Soliman et al. 
(1968), and Bae et al.(1971)] assumed that the interfacial shear stress in condensation 
is similar to that in an adiabatic two-phase flow. Therefore, some correlations 
normally applicable to adiabatic two-phase flows were used to determine the 
interfacial shear stress in condensation. In doing so, the effect of the momentum 
transfer caused by the mass transportation of the vapor condensed into liquid was 
neglected. 
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However, the aforementioned treatment on the interfacial shear stress was 
regarded as inappropriate by Shekriladze and Gomelauri (1966). They suggested that 
the interfacial shear stress depended mainly on the effect of the momentum transfer 
due to the condensation of vapor into condensate. For high condensate rates, the 
interfacial shear stress is presumably approximately equal to the value of this 
momentum transfer. 

Linehan et al (1969) suggested that the interfacial shear stress in the presence of 
condensation should be expressed as the sum of the adiabatic shear stress plus the 
product of the condensation rate and an average vapor velocity. This theory was 
tested in their experiments involving steam condensation inside a rectangular test 
section. The calculated local surface temperature and mean film thickness based on 
this theory was found to be close to the experimental results. 

Jensen and Yuen (1982) used an expression of the interfacial shear stress of a 
two-phase flow without phase change to represent the term of the adiabatic shear 
stress in the Linehan model. Expression for the interfacial shear stress was given for 
the case of a smooth interface as well as for a wavy interface. For high condensation 
rates, the value of the momentum due to the condensation of the vapor was found to 
be the dominant factor in the expression of the interfacial shear stress, supporting the 
assumption of Shekriladze and Gomelauri (1966). 

2.2.4. Experimental Studies of Condensation Inside a Duct 

Akers et al. (1958) measured heat transfer coefficients with R-12 and propane 
inside a horizontal tube. The results showed that increasing the vapor velocity 
enhanced the average heat transfer coefficient. When the vapor velocity increased 
beyond a certain value, the average heat transfer coefficient exhibited a nearly linear 
relationship with the velocity. The experimental data at high vapor velocities were 
plotted on the basis of the parameters derived from a single phase heat transfer 
correlation proposed by Eckert (1950). The data with high vapor velocities correlated 
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with Eckert equation. Akers et al. also investigated the effect of the differences 
between the saturated vapor temperature and the condensing surface temperature (AT) 
on the condensation. It was found that effect of AT on the average heat transfer 
coefficient was not clear. 

Chato (1962) conducted experiments with condensation of R-113 vapor flowing 
inside a 0.72 m long horizontal tube with 14.5 mm mean I.D. The condensate angle 
(Fig. 2.2) at the end of the test section was measured by visual observation and 
photography. The results of the experiment showed that increasing the slope of the 
test section to about 0.01 can lead to an increase in heat transfer rates. Further 
increasing the inclination did not increase the heat transfer rates appreciably. The 
experiment also showed that the orientation of the tube was not important at high 
vapor velocities. 

Rossen and Myers (1965) measured the peripheral variation of local condensation 
heat transfer coefficients with methanol and acetone inside a horizontal tube using a 
heat meter. The tube was divided into a number of small areas. The heat flux 
measured from each area was defined as the local heat flux. The local heat transfer 
coefficient was then found from the measured local heat flux. 

Goodykoontz et al (1966, 1967a, and 1967b ) described a method for measuring 
the local heat transfer coefficient with R-113 and steam condensing inside a vertical 
tube. The axial temperature gradient of the coolant in the cooling jacket was 
measured using thermocouples. Neglecting axial conduction in the tube, the local 
heat transfer rate was determined by computing the product of the temperature 
gradient, the coolant flow rate, and the specific heat of coolant. The local heat 
transfer coefficient was calculated from the local heat transfer rate and the 
temperature difference between the saturation temperature of the vapor and the inner 
wall temperature. 

Abis (1969) measured local heat transfer rate for forced convection condensation 
of R-12 inside a 2.44 m long, 12.7 mm I.D. horizontal tube. The vapor velocity was 
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considered sufficiently high (refrigerant mass velocity, 135.62-446.20 kg/m 2 s) to 
produce an annular flow pattern throughout the entire length of the tube. The local 
heat transfer coefficient was determined using a method similar to that of 
Goodykoontz et al (1966). The results showed that the local heat transfer coefficient 
decreased along the tube and increased when the vapor mass velocity increased. The 
effect of AT on the heat transfer coefficient was also studied, but no conclusion was 
drawn due to some conflicting results. 

Using R-12 and R-22, Bae et al. (1971) measured the condensation heat transfer 
coefficients in the annular regime. The test section of nickel tube was 5.5 m long and 
had a 12.5 mm I.D. They compared their experimental data with those of Akers and 
Rosson (1960), Altman et al. (1959), and Chen (1962) using the coordinates 
suggested by Akers and Rosson (1960). Bae et al. found that their heat transfer 
coefficients were generally higher than those of the other researchers. They 
concluded that this was because the coordinates used in the comparison did not 
include the parameter of quality; since the quality was related to the thickness of 
condensate film (a major resistance to the heat transfer), any difference in the quality 
would lead to different heat transfer coefficients even when the vapor velocity was the 
same. 

Sardesai et. al. (1981) reported measurements of the peripheral local heat transfer 
coefficients inside a 2.92 m long, 24.4 mm I.D. horizontal stainless steel tube, using 
R- 113, steam, propanol, methanol, and n-pentane as the condensing fluids. Using 
thermocouples embedded in the tube wall, the local heat transfer coefficient was 
measured. The measured peripheral local heat transfer coefficient was used to 
determine the flow regimes. If the distribution of the local heat transfer coefficient 
was about uniform around the periphery of the tube, the flow regime was considered 
to be an annular flow; if the local heat transfer coefficient at the top of the tube was 
much higher than that at the bottom of the tube, the flow regime was considered to be 
a stratified flow. From the analysis of the experimental data, they found that the 
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transition between the annular flow and the stratified flow regimes was determined by 
a parameter which was a function of the Martinelli number and the Froude number. 

Kutsuna et al. (1985) reported experimental results of condensation of vapor flow 
inside a horizontal rectangular duct. The test section was 100 mm wide, 30 mm high, 
and 500 mm long. The bottom plate of the duct was the condensing surface. The 
local heat flux was determined from a measurement of the temperature distribution in 
the condensing surface. The average heat flux was measured from the enthalpy of the 
condensate in the test section. The experimental results of the local heat transfer 
coefficient for pure steam showed a rapid decrease of the local heat transfer 
coefficient along the duct. 

Christodoulou (1987) experimentally investigated the film condensation of a vapor 
flow inside a horizontal rectangular duct with R-113 condensing on the bottom plate 
of the duct. The duct was 40 mm wide, 25 mm high, and 1 m long. The condensing 
surface was made of a copper block, which was hollowed out into three separate, 
equal sections. Cooling water was circulated in an channel formed by those three 
sections. In the experiment, Christodoulou measured the film thickness of condensate 
using an ultrasonic technique, the sectional heat transfer coefficients, and the total 
average heat transfer coefficient. The results showed that the thickness of the 
condensate film increased with increasing distance from the entrance of the duct and 
deceased with increasing vapor velocity; the sectional heat transfer coefficients 
decreased with increasing distance from the entrance. However, the effect of vapor 
velocity on the total average heat transfer coefficient was not clear. Since the 
experiment was conducted within a limited range of inlet vapor velocity (0.15—0.45 
m/s), no interfacial instabilities (interfacial waves) were observed in the experiment. 

Barry and Corradini (1988) experimentally studied film condensation of a steam- 
air mixture in the presence of interfacial waves. The test section was a horizontal 
square duct of 102 mm wide, 102 mm high, and 1.83 m long. The bottom plate of 
aluminum in the test section was used as the condensing plate. To develop a 
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technique for measuring the wave parameters, an isothermal experiment using air- 
water flow through the duct was carried out before the condensation experiment. The 
wave frequency, celerity, length, amplitude, and film thickness were measured using 
hot-wire anemometry and probabilistic analysis. The local heat transfer coefficient 
along the length of the condensing surface was determined at four locations by 
measuring the temperature gradient along the condensing wall. The results of the 
condensation experiment showed that the local heat transfer coefficient decreased 
along the condensing surface and increased with increasing vapor velocity. 

Comparing the results of the condensation experiment and the isothermal experiment, 
they concluded that condensation may have a damping effect on wave initiation. 
However, the effects of interfacial waves on condensation heat transfer are not shown 
in the study. 

2.3. Wave Initiation in Two Phase Flows 

It has been observed that condensation heat transfer is considerably influenced by 
the flow regimes. Bell et al. (1970) compared the results from various correlations 
for condensation inside horizontal tubes and found that the difference between the 
various correlations was fairly large, varying by a factor of 10 or more for a given set 
of conditions. For condensation in the annular flow regime, the condensate film may 
be smooth or wavy, depending on the magnitude of the vapor velocity. According to 
Dukler (1977), the interfacial waves in the annular wavy flow regime substantially 
increased the transport of energy and momentum in both the vapor and liquid phases. 
This indicates that the condensation heat transfer rate in the annular flow regime with 
a smooth condensate film may be different from that with interfacial waves. 

Therefore, a study of the instability causing the initiation of the interfacial waves is 
important. 

The instability of a liquid flowing horizontally adjacent to a gas flowing parallel 
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to the liquid but with a different velocity was recognized by Helmholtz over a 
hundred years ago. Later, Kelvin also studied this instability. In their studies, the 
domain of each of the fluids was semi-infinite and viscous effects were ignored. 
Chandrasekhar (1961) proposed a criterion to predict the instability. It is shown that 
both surface tension and gravity are stabilizing, but if either of them is neglected, the 
configuration is unconditionally unstable. 

When the fluids, instead of being semi-infinite, are flowing between two parallel 
plates, the condition of instability, including surface tension (but neglecting viscous 
effects), is given by Milne-Thomson (1969). It is shown that the surface tension is 
significant only for short wave lengths corresponding to high wave numbers. 

Feldman (1957) analyzed the stability of a liquid film with a very small thickness, 
in which the shear stress distribution was assumed to be linear, with a semi-infinite 
fluid on top of it; the two fluids flow with different velocities. He included the effect 
of the shear stress at the interface and gave the conditions for stability. From the 
results of his analysis, Feldman concluded that (in the following, subscripts 1 and 2 
represent the lower fluid and the upper fluid) 


a. If gravity and surface tension are ignored: 

1. for very small ratios of p 2 /p\, the flow is completely stable 

2. for a given value of pjp 2 , there is always a value of p 2 lp x for which the 
flow is unstable 

3. as p 2 /Pi 0, the ratio of the wave speed to interfacial velocity approaches 
approximately 0. 1 

4. the critical conditions depend only on the liquid flow rate and not on the gas 
flow rate except for the change in the interfacial shear stress 

5. all waves travel at speeds less than the velocity of the liquid— gas interface 


2.16 


2. Literature Review 


b. If gravity and surface tension are significant, they are destabilizing. 


Conclusions 2 and 3 seem to contradict conclusion 1. For air flowing over water 
at 100 °C, Feldman quoted a critical liquid Reynolds number (p,hiU,//X|) of 60,000 
which appears to be much greater than the observed values. The physical conditions 
of 100 °C water in contact with air at atmospheric pressure is not realistic as 
vaporization may become significant. Conclusion b is in direct contradiction with the 
results given by Chandrasekhar (1961) and Milne-Thomson (1969). 

The effect of viscosity stratification in two fluids was studied by Yih (1967). He 
pointed out that in both Couette and Poiselle flows, the interface was always unstable, 
no matter how small the Reynolds number was. He also indicated that because, in 
practice, the upper fluid is lighter than the lower fluid and due to the stabilizing effect 
of gravity, such instabilities caused by viscosity stratification alone may not be 
discernible. His main point is that, when considering instabilities in such flows, 
viscosity stratification should be considered and that one of the parameters defining 
instabilities should be the Reynolds number. 

The only study that takes into account phase changes in the flow appears to be 
that made by Kocamustafaogullari (1985). From a simplified equation for the 
stability, he derived a criterion for stability, which included the effects of gravity, 
surface tension, viscosity, and phase change. For condensation on a horizontal plate, 
the criterion predicted unconditional stability, (which contradicts experimental 
observations). 

On the experimental side, Hanratty (1982) summarized the results of experiments 
conducted in a horizontal, enclosed channel with air and water as the fluids. For very 
low velocities of air, the interface was smooth. As the air velocity increased, long 
crested two-dimensional waves with wave lengths of 2. 2-3.0 cm and with wave 
velocities greater than the liquid velocity appeared. Kao and Park (1972) reported 
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results of experiments directed towards determining the effect of the presence of a 

layer of fluid in laminar flow (at low Reynolds numbers) on top of another fluid at 

much higher Reynolds numbers. They observed that the shear instability in the lower * 

fluid occurred at a Reynolds number of 2300 and the disturbance at the interface 

appeared when the liquid became turbulent. From their studies on condensation 1 

inside tubes, Soliman and Berenson (1970) reported that film instabilities were 

observed to originate very close to the inlet of the condenser tube and the instabilities 

appeared as irregularly shaped wave trains. They speculated the waves to be a 

combination of Kelvin-Helmholtz and Tollmein-Schlichting types. 

2.4. Concluding Comments 

From the survey of the available literature, it is found that a detailed study of 
condensation in the annual flow regime, especially condensation in the annular wavy 
flow regime, has not been done. Nearly all the analytical studies are based on the 
assumption of smooth film condensation. Several models also include constants 
which have to be determined from experimental results. Some models are based upon 
the analogies between heat and the momentum transports, depending upon the 
available pressure drop correlations for adiabatic two-phase flow. None of the 
analytical equations is supported by broad experimental verifications. 

In addition to inadequate analytical work, the experimental work on condensation 
in a horizontal tube has also not been sufficient. Some researchers report only the 
average heat transfer coefficients, due to the relative ease of measuring them. Even 
though some researchers report experimental results of local heat transfer coefficient, 

* 

in reality, the local value is actually an average value over a finite area of the 
condensing surface. Also, the experimental results about the effects of the interfacial 
waves on the heat transfer of condensation have not been reported. 
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A survey of the literature reveals that the study on condensation in an annular 
flow regime is limited; more work is required in both the analytical modeling and the 
experimental investigation to understand and predict heat transfer in condensation. 
The objectives of the research are as follows: 

— to measure film thickness, local heat transfer coefficient, average heat 
transfer coefficient, wave length, and wave speed of condensation of a 
vapor flow inside horizontal rectangular duct 

— to identify the dimensionless parameters that are significant in the 
condensation heat transfer 

— to study the initiation of interfacial waves with phase change (condensation) 
and without phase change (air-liquid flow with liquid viscosity variation) 

— to develop a simple analytical model for condensation of a vapor flow 
inside a horizontal rectangular duct 
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i 4.1. Introduction 

The two main purposes of the experiment were to measure the heat transfer 
coefficients during condensation of a vapor flowing inside a horizontal rectangular 
duct (with the vapor condensing only on the bottom plate of the duct), and to study 
the effect of the interfacial waves on the heat transfer coefficient. Two separate 
experimental set-ups were used in the experiments. One was a condensation set-up 
for measuring the heat transfer with a smooth thin condensate film or with interfacial 
waves. To assist in the study of interfacial wave initiation, experiments were also 
conducted in an adiabatic air-liquid flow set-up. Since the geometric shapes and 
dimensions of the test sections in both experimental set-ups were identical, it is 
possible to compare the results of the interfacial wave initiation in the condensation 
experiment with those in the adiabatic air-liquid flow. 

The experimental set-ups and procedures are separately described for 
condensation, and wave initiation in an adiabatic air-liquid flow. 

4.2. Experimental Set-up of Condensation 

The experimental set-up of condensation was a modification of a previous work 
* by Christodoulou (1987) and is schematically shown in Figs. 4.1 and 4.2. The set-up 

included a condensation loop and a coolant loop. In the condensation loop, vapor was 
generated from a boiler and admitted to the rectangular duct of the test section 
through a converging unit. Part of the vapor was condensed on the bottom 
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Fig. 4.1 Schematic of experimental set-up of condensation 
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condensing plate inside the test section, while the rest of the vapor was discharged 
from the test section and was completely condensed in an auxiliary condenser. The 
condensates from the auxiliary condenser and the test section were gathered in a main 
liquid tank, from which the condensate was pumped back to the boiler. In the coolant 
loop, cooling water, which was provided by a temperature controller, flowed into the 
cooling channel of the test section via a set of rotameters. After passing through the 
cooling channel, the cooling water returned to the temperature controller. 

The experimental set-up was assembled and frequently checked for leakage. 

The following sections describe the apparatus used in the condensation set-up. 

4.2.1. Test Section 

The schematic of the test section is shown in Fig. 4.3. Details of the cross 
section of the test section are presented in Fig. 4.4. 

The test section consisted of a rectangular duct and a cooling channel. The inside 
dimensions of the rectangular duct were 0.04 m wide, 0.025 m high, and 1 m long. 
The bottom condensing surface was of 0.914 m long, 0.04 m wide, and 6.4 mm thick 
copper plate. The side and top plates of the duct were of transparent polycarbonate 
sheets, which permitted visual observation of the flow patterns of the condensate. 

The condensate metering tube at the exit of the test section, which was also used as a 
drain for all the condensate in the test section, had an inside diameter of 6.4 mm and 
was connected to a small piece of a stainless steel plate attached to the end of the 
copper condensing plate. The top of the stainless steel plate was flush with the 
surface of the copper plate. 

The cooling channel was underneath the condensing surface. The side plates of 
the channel were of brass; the bottom plate of the channel was of polycarbonate. The 
channel was divided into three equal sections by 3.2 mm thick plexiglass separators. 
Each section was 305 mm long and 17.5 mm deep. A mixing chamber (Fig. 4.5) was 
placed at the outlet of each section. The purpose of the mixing chamber was to 
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Fig. 4.3 Schematic of test section of condensation 
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Fig. 4.4 Cross section of test section of condensation 
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assure the measurement of the correct mean coolant temperature by mixing the 
coolant in the chamber. 


13.89 42.86 




Dimensions in mm 


Fig. 4.5 Details of mixing chamoer 


The test section was mounted horizontally. A liquid level was used to ensure that 
the test section was horizontal. To allow for the thermal expansion of the tubes 
connected to the test section, flexible tubes were installed at the inlet and outlet of the 
test section. The test section was insulated with fiber glass and styrofoam sheets. 
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4.2.2. Converging Unit 

To provide a uniform velocity distribution of vapor, a converging unit was 
installed at the inlet of the test section (Fig. 4.6). The converging unit consisted of a * 

rectangular section and a reducer section. The walls of the two sections were of 3. 18 
mm thick brass plates. The inside dimensions of the rectangular section were 101.6 
mm wide, 101.6 mm high, and 200.02 mm long. The length of the reducer section 
was 50.8 mm. The inside dimensions of the outlet of the reducer section were 40 
mm wide and 25 mm high, which were identical to those of the test section. To 
assist in a uniform flow, a flow deflector was installed in the converging unit near its 
inlet and thin glass tubes were packed in the rectangular section. 
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Fig. 4.6 Details of converging unit 
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The effectiveness of the converging unit was tested with air. The testing 
procedure and results are detailed in Appendix B, where it is shown that the velocity 
profile at the outlet of the converging unit was mostly uniform except near the side 
plates of the outlet where the velocities were lower than the average velocity due to 
the wall friction. In the test, the average deviation of the measured velocity at the 
outlet of the unit was ±9.10 %. 

4.2.3. Boiler 

Vapor 



* 

Details of the boiler are shown in Fig. 4.7. The 129 liter capacity boiler was 
fabricated from stainless steel. The maximum designed working pressure of the 
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boiler was 30 psi. The boiler usually contained about 23 liters of liquid. The liquid 
was heated by two 3 kw electrical heating elements mounted at the bottom of the 
boiler. A brass mesh screen covering the upper cross section of the boiler was used 
to reduce the liquid carryover. The mesh screen was located about 50 mm below the 
top of the boiler and about 80 mm above the normal level of liquid in the boiler. 

The temperature of the boiler was measured by a copper-constantan 
thermocouple. The pressure of the boiler was measured by a Bourdon Type Gauge 
through a pressure tap mounted on the top of the boiler. The level of the liquid in the 
boiler was observed through a view glass tube installed vertically on the side of the 
tank. This level was also automatically monitored by a level controller with a float 
on the surface of the liquid. If the liquid height was lower than 155 mm, the 
controller shut off the power supply to the heating elements until the normal level of 
fluid was restored. If the pressure of boiler was more than 70 kPa above the 
atmosphere pressure, a safety valve installed on the top of the boiler opened to relieve 
the pressure. Prior to its operation, the boiler was hydraulically tested at a pressure 
of 414 kPa. 

4 . 2 . 4 . Auxiliary Condenser 

Details of the auxiliary condenser are shown in Fig. 4.8. The shell of the 
condenser was of stainless steel; the tubes inside the condenser were of copper. The 
vapor flowing in the shell side was cooled by building water supply flowing inside the 
tubes. The total heat transfer surface of the condenser was 0.24 m 2 . 

* 

4 . 2 . 5 . Temperature Controller and Gravity Tank for Coolant , 

A temperature controlling unit was used to supply cooling water to the test 

I 

section. This unit included a heating circuit and cooling circuit. The controller had a 
capacity of delivering the cooling water to the test section at a rate of 0 to 7.6 liters 
per minute with an operating range from -30 °C to +70 °C. The temperature 
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Fig. 4.8 Details of auxiliary condenser 
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controller also served as a large coolant storage tank (=46 liters) during experiment. 

The gravity tank was a 0.5 m high and 0.36 m in I.D. container filled with water 
and connected to the temperature controller through a valve. When the experiments 
were completed for the day and the temperature controller was shut off, the valve 
connecting the gravity tank and the cooling loop was opened. Since the tank was 
placed approximately 1 m above the test section, the tank provided positive water 
pressure to the cooling channel in the test section to prevent any air leak into the 
channel from the environment. 

4.3. Experimental Set-up of Adiabatic Air-liquid Flow 

The shape and the dimensions of the test section in the experimental set-up of the 
adiabatic air-liquid flow were identical to those of the test section for condensation. 
Therefore, it is expected that a comparison of the result of the wave initiation with the 
air-liquid flow with the results of the wave initiation with condensation will lead to a 
better understanding of the effect of condensation on wave initiation. 

Fig. 4.9 is a schematic of the experimental set-up. A photograph of the set-up is 
provided in Fig. 4. 10. The set-up consisted of an air passage connecting to the 
building compressed air line, and a liquid loop. During experiments, the air from the 
building supply passed through a compressed air reservoir, a converging unit, and an 
inlet section into the test section. The air flow rate was controlled by adjusting the 
pressure regulators at the inlet and the outlet of the reservoir and was measured by a 
venturimeter. In the liquid loop, the liquid from a gravity tank entered the test 
section. The liquid from the test section drained to a collecting tank, from where the 
liquid was pumped back to the gravity tank. When the gravity tank was full, the 
liquid overflowed to the collecting tank. 

Fig. 4. 1 1 illustrates the inlet section. The inside dimensions of the inlet section 
were same as those of the test section (40 mm wide and 25 mm high). The inlet 
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Fig. 4.9 Schematic of experimental set-up of air-liquid flow 



Fig. 4.10 Photograph of experimental set-up of air-liquid flow 
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Dimensions in mm 


Fig. 4.1 1 Details of inlet section 

section was divided into two parts (an air channel and a liquid channel) by a H shaped 
horizontal separator inserted inside the section (Fig. 4-11). Based on the way of 
installing the separator, the height of the liquid channel was changed to either 1 mm 
or 3 mm. 

4.4. Measurements in Condensation Experiment 

The measured variables in the condensation experiment were 


volumetric flow rate of vapor 
volumetric flow rate of condensate 
volumetric flow rate of coolant 
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— temperatures of the condensing surface at different locations, vapor, and 
coolant at the inlet and the outlet of each section 

— vapor pressure in the test section 

— thickness of condensate film at five locations along the condensing surface 

— wave length and the wave speed in the condensate 

The various instruments used in the experiments are described in the following 
sections 

4 . 4 . 1 . Measurement of Vapor Flow Rates 

The vapor flow rates at the inlet of the test section was measured by a venturi 
meter and a rotameter. A turbine flow meter (turbine meter #1 in Fig. 4.1) was used 
to measure the condensate flow rate to the boiler. The rotameter and the venturi 
meter (venturi# 1) were installed in the vapor line between the boiler and the inlet of 
the test section; the turbine flow meter (turbine# 1) was installed in the condensate line 
between the inlet of the boiler and the outlet of the main liquid tank. 

The rotameter and the venturi meter were calibrated in the laboratory using a 
laminar flow meter with air as the fluid. The turbine flow meter was calibrated by 
the manufacturer. The procedures and results of the laboratory calibrations of the 
rotameter and the venturimeter are given in Appendix C. 

In the experiment, the vapor flow rate measured by the rotameter was used in the 
data analysis, while the measurements from the venturimeter and the turbine flow 
meter were only used as a check on the measurements with the rotameter. The 
differences were usually in the range of ±5-10 %. 

4 . 4 . 2 . Measurement of Condensate Flow Rate 

The condensate flow rate from the test section was measured by a gauge glass 
installed at the outlet of the test section (Fig. 4.1); details of the gauge glass are 
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shown in Fig. 4.12. The gauge glass was provided with a quick shutoff valve at its 
bottom. The glass tube was 455 mm long with 13 mm I.D. When the quick shutoff 
valve was closed, the time required for the condensate to fill the gauge glass from the 
lower marker to the upper maker (Fig. 4. 12) was measured; the temperature at the 
outlet of the gauge glass was measured by a thermocouple. From those 
measurements, the flow rate of the condensate was determined. 
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Fig. 4.12 Details of gauge glass 

The flow rate of the condensate was also determined as the difference in the 
vapor flow rate into and out of the test section as measured by the venturimeters. A 
turbine flow meter (turbine meter#2 in Fig. 4.1) was installed at the outlet of the 
liquid tank collecting the condensate from the test section. It also measured the 
condensate flow rate from the test section. 

The difference between the condensate measured by the gauge glass and that 
measured by the venturimeter varied from ±5 % to ±15 %; the difference between 
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the gauge glass and the turbine flow meter was ±5 % to ±20 %. The difference 
from the turbine flow meter was due to the variation of the liquid level in the liquid 
■* tank connecting the inlet of the turbine flow meter. Since the liquid level was 

influenced by the pressure in the test section, the vacuum in the auxiliary condenser. 

1 and the speed of the micro-pump in the condensate line, it was somewhat difficult to 

maintain a constant liquid level. Rate of condensation as measured by the gauge glass 
was used in the data analysis. 

4.4.3. Measurement of Cooling Water Flow Rate 

Cooling water flow rate was measured with a set of four rotameters in parallel. 
The rotameters were calibrated in the laboratory by weighing the water collected in a 
measured time interval. Results of the calibration are given in Appendix C. During 
experiments, the measurements from the rotameters was checked at regular intervals 
by weighing the coolant collected in a measured time. 

4.4.4. Measurement of Temperatures 

The following temperatures were measured in the experiment 

— condensing surface 

— vapor in the test section 

— condensate 

— cooling water in the cooling channel 

— vapor through the flow meters 

— surface of vapor tube, converging unit, and test section 

M — boiler and environment 

* 

I Temperatures of Condensing Surface 

The temperature of the condensing surface was measured by eight thermocouples 
embedded in the surface of the condensing plate. The layout of the thermocouples is 
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Fig. 4.13 Details of condensing surface 
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shown in Fig. 4.13. Each thermocouple junction was placed near the condensing 
surface through a 1.2 mm I.D. hole in the copper plate and soldered to the 
* condensing surface. The residues of the soldering material on the surface were 

cleaned by using fine sand paper and Brasso (a metal polish solution). The space 
i between the wire and the hole (3.2 mm I.D) in the plate was filled with electrical 

cement. 

Temperatures of Vapor in the Test Section 

The temperature of vapor in the test section was measured by four equally spaced 
thermocouples placed along the center line of the test section (Fig. 4. 3). The distance 
between two consecutive thermocouples was 280 mm. The distance from the entrance 
of the test section to its nearest thermocouple was 50 mm. The distance between the 
thermocouple hot junction and the side wall was 10 mm (Fig. 4.14). 



Thermocouple 


Thermocouple 
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Fig. 4-14 Arrangement of thermocouples for measuring vapor temperature 


Temperature of Condensate 

The temperature of the condensate in the gauge measuring glass was measured by 


4.19 


4. Description of Experiments 


a thermocouple. The location of the thermocouple is shown in Fig. 4.12. 
Temperatures of Cooling Water 

The temperatures of cooling water were measured using thermocouples installed 
inside the inlet and the outlet of the mixing chambers. The temperature difference of 
cooling water across each section of the cooling channel was measured using three 
pairs of differential thermocouples placed inside each mixing chamber. The total 
temperature change of the cooling water across the whole cooling channel was also 
measured by a pair of differential thermocouples placed at the inlet and the outlet of 
the channel. This total temperature change was compared with the value of the sum 
of the three temperature changes measured at each section of the cooling channel. 

The difference was usually within ±2 %. The arrangement of the differential 
thermocouples is shown in Fig. 4.15. 
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Fig. 4.15 Arrangement of differential thermocouples 
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Temperatures of Vapor at the Outlet of Rotameter and Venturi Meters 

For determining the vapor mass flow rate using the rotameter and venturi meters, 
the temperatures of vapor at the outlet of the flow meters were needed. The 
temperatures of vapor at the outlet of the flow meters were measured using 
thermocouples. 

Surface Temperatures of Vapor Tubes. Converging Unit, and Test Section 
To prevent condensation of the vapor on the walls of the vapor tubes, the 
converging unit, and the test section, the surface temperatures of these parts were 
maintained at a slightly higher temperature than the saturation temperature of the 
vapor using additional heating elements. The heating elements were either wrapped 
around the tubes and the converging unit or placed on the top of the test section. The 
power to the heating elements was adjusted, so that the vapor was superheated by 3-5 
°C. The surface temperatures of the vapor tubes, the converging unit, and the test 
section were measured using thermocouples. 

All thermocouples used in the experiment were made of 28 gauge copper- 
constantan wires and were calibrated. A description of the calibration is given in 
Appendix D. 

4.4.5. Measurement of Vapor Pressure 

The vapor pressure in the test section was measured using a bourdon tube 
pressure gauge. The position of the pressure tap is shown in Figs. 4.3 and 4.4. The 
inside diameter of the pressure tap was 1.6 mm. The tube (made of transparent 
teflon) connecting the pressure gauge with the pressure tap was installed vertically and 
was heated slightly by an auxiliary heating element. If any vapor condensed in the 
tube, the condensate flowed down into the test section without blocking the tube. 

The boiler pressure was measured by a bourdon tube pressure gauge. The 
atmospheric pressure was measured by a barometer in the laboratory. The vacuum in 
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the auxiliary condenser was observed through a vacuum meter. 

4.4.6. Measurement of Film Thickness 

An ultrasonic transducer was used to measure the condensate film thickness at 
five locations along the condensing surface. During the initial test of the transducer, 
it was found that with a copper plate, the attenuation of the reflected signal was 
significant. The signal strength was improved by using an aluminum plug. As a 
result, aluminum plugs were inserted in the copper plate at each location where the 
film thickness was to be measured. The aluminum plugs were made flush with the 
copper plate. The locations of the aluminum plugs are shown in Fig. 4.13. 

Fig. 4.16 illustrates the ultrasonic measurement equipment. The ultrasonic 
transducer sent a 20 MHZ signal to the condensate film through the aluminum plug. 
The signal was first reflected at the interface of the transducer and the aluminum 
plug, next at the interface of the condensate film and the condensing surface, and then 
at the interface of condensate and vapor. In Fig. 4.16, the peak numbered 1 




Ultrasonic Analyser 


Fig. 4.16 Illustration of film thickness measuring system 
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represented the reflection at the interface of the transducer with the bottom surface of 
the aluminum plug; peak number 2 was the signal reflected from the aluminum plug- 
condensate interface; peak number 3 was the reflection of the signal at the 
condensate-vapor interface. The interval (t f ’) between number 2 and 3 was the time 
for the signal to pass through the film to reach the vapor-liquid interface and then 
reflected from there. Having t f ’ and a, (sonic velocity of liquid), the film thickness 
was determined. The procedure of using t f ’ and a, to find the film thickness i> 
described later in Section 4.9.1. 

A typical trace of the ultrasonic signal displayed on the oscilloscope screen is 
shown on Fig. 4.17. A description of calibration of the ultrasonic transducer (for 
obtaining the sonic velocity of the condensing fluid) is given in Appendix E. 



Fig. 4.17 Typical pulse trace 


4.4.7. Measurements of Wave Length and Wave Speed 
Wave Length 

The wave length was determined by measuring the distance between two 
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consecutive crests. In the experiment, photographs of interfacial waves were taken 
with a ruler placed along the side of the test section. From the photographs, the wave 

length was measured (Fig. 4.18). ► 

% 


5 3 / 38 39 40 41 42 43 44 45 46 4 / 



Fig. 4.18 Photograph of wave length (The unit of the scale in the picture is mm.) 


Wave Speed 

The wave speed was measured using a mechanical device shown in Fig. 4. 19. 
The device was a moving belt driven by a variable speed motor. The top surface of 
the belt was marked with parallel white strips. During the measurement, the device 
was placed next to one of the sides of the test section, the speed of the motor was 
adjusted so that the strips on the belt moved at the speed of the waves. The speed of 
the belt was determined by measuring the rotational speed of the motor. 
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Fig. 4.19 Schematic of wave-speed measuring device 
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4.5. Measurements in Air-liquid Flow Experiment 

A The measured variables in the experiment of air-liquid flow were 

— volumetric flow rate of air; 

* — mass flow rate of liquid; 

— density of liquid; 

— viscosity of liquid; 

— thickness of liquid film. 

The procedure of measuring these variables is described in the following sections. 

4.5.1. Measurement of Mass Flow Rate of Air 

The mass flow rate of air was measured using a venturimeter. The venturimeter 
was calibrated in the laboratory. The calibration results are given in Appendix C. 

4.5.2. Measurement of Mass Flow Rate of Liquid 

The mass flow rate of liquid was measured by weighing the liquid collected in a 
known time. The liquid was collected at the outlet of the test section. 

4.5.3. Measurement of Density of Liquid 

The density of liquid was measured using a volumetric method. The volume of 
the measuring flask was 500 ml. 

* 4.5.4. Measurement of Viscosity of Liquid 

In the experiment of air-liquid flow, the liquid viscosity was varied by mixing 
( glycerin with water. The viscosity of liquid was measured using a Falling Ball 

Viscometer. As shown in Fig. 4.20, the main components of the instrument are a 
glass tube inclined at 10° with respect to the vertical and a set of calibrated balls for 
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Fig. 4.20 Schematic of viscosity meter 


different viscosity ranges. Two ring marks (A and B) on the tube are spaced 100 mm 
apart. During the measurement, a selected ball was placed inside the tube. The time 
for the ball to fall from marks A to B along the tube was recorded. Based on the 
falling time, the density of the liquid, and the ball and its drag coefficient (both 
provided by the manufacture), the viscosity of the liquid was determined. 

The Falling Ball Viscometer was calibrated in the laboratory using distilled water 
and different standard liquids provided by the manufacturer. A comparison was made 
between the measured viscosity of the liquid (mixture of glycerin and water) and the 
value reported by Leffmgwell (1945). The average difference of the comparison was 
within ±7 %. 

k 

4.6. Data Acquisition System 

1 

A data acquisition system was used in the condensation experiment to obtain the 
signals from the pressure transducers and most of the thermocouples. The system 
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consisted of a DT-100 data taker (made by Data Electronics Pty Ltd) and a P.C. The 
data taker had 23 differential channels or 46 single channels. The maximum sample 
% rate of the data taker was 30 samples/ second. The accuracy of the data taker was 

within ±0.15 %. 

1 

4.7. Liquids 

4.7.1. Liquids in Condensation Experiments 

R-113 (trichlorotrifluoroethane) and FC-72 (one of fluorinert™ electronic liquids 
developed by the 3M company) were used as the condensing liquids. The liquids 
were selected because of their low saturation temperature at atmospheric pressure (T sat 
= 47 °C), wetting properties, and moderate toxicity. The properties of R-113 were 
taken from the ASHRAE Handbook (1985); the properties of FC-72 were obtained 
from the 3M Product Manual for Fluorinert™ Electronic Liquids (1989). The 
thermodynamic properties of R-113 and FC-72 in the range of saturation temperatures 
employed in the experiments are given in Table 4.1. 

4.7.2. Liquid in Air-liquid Flow Experiments 

The liquid in the experiments of air-liquid flow was mixtures of distilled water 
and glycerin mixture. The glycerin was 99.9% chemically pure. The viscosity of the 
water- glycerin varied from 1 x 10' 6 kg/m s at 0 % weight concentration to 7.6 x 
10' 6 kg/m s at 60 % weight concentration. 

* 4.8. Experimental Procedures 

{ 4.8.1. Procedure for Condensation Experiments 

The procedure of the experiments had 3 major steps: (1) preparation of 
experiments; (2) initial stage of experiments; (3) measurements. 
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1 


Table 4.1 Selected thermodynamic properties of R-1 13 and FC-72 


R-1 13 

T Mt 

P 

P '« 

Pi 

h* 

ki 

c p> 

Pr, 

°C 

KPa 

kg/m 3 

kg/m 3 

kJ/kg 

mw/k m 

J/kg k 


40 

78.0 

5.8 

1526.1 

146.3 

72.0 

975 

7.7 

50 

109.3 

8.0 

1501.7 

142.8 

70.1 

986 

7.1 

60 

149.8 

10.8 

1476.8 

139.7 

68.0 

994 

6.6 

70 

218.1 

14.3 

1451.2 

136.2 

65.8 

1004 

6.2 

FC-72 

Tsat 

P 

P g 

Pi 

h* 

k, 

c p< 

Pr, 

°c 

kPa 

kg/m 3 

kg/m 3 

kJ/kg 

mw/k m 

J/kg k 

40 

54.8 

7.6 

1640.0 

89.7 

55.6 

1064 

6.5 

50 

79.5 

10.6 

1621.9 

86.9 

54.5 

1089 

6.0 

60 

112.4 

14.7 

1604.2 

84.1 

54.4 

1101 

5.4 

70 

154.9 

20.0 

1583.9 

81.1 

52.3 

1122 

4.9 
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Preparation of Experiments 

The preparation steps included removing air from the condensation loop, warming 

^ the test section and vapor tube, and heating the liquid in the boiler. 

» 

Some air would leak into the condensation loop during the period between 
1 experiments when the loop was under high vacuum. The air was removed from the 

condensation loop using a vacuum pump. To prevent vapor condensation on the wall 
of the test section and the vapor tubes, the test section and the vapor tubes from the 
boiler to the test section were heated to about 5 °C above the expected vapor 
saturation temperature in the test section using the heating elements placed on the top 
of the test section, and wrapped around the vapor tubes and the converging unit. The 
condensing liquid was heated in the boiler. When a slightly positive pressure was 
established in the boiler, a purge valve was slowly opened to discharge any air that 
might have leaked into the boiler. 

Initial Staee of Experiments 

The initial stage included adjusting the vapor flow rate to the inlet of the test 
section, the condensate flow rate to the boiler, the cooling water flow rate, and 
reaching a steady state. The vapor flow rate was gradually increased to a 
predetermined value by slowly opening a valve at the outlet of the boiler when the 
gauge pressure of the boiler reached about 14 kPa. The vapor from the boiler passed 
through the flow meters and the converging unit into the test section where a part of 
the vapor was condensed. In the test section, the temperature of the vapor was 
maintained approximately 3 °C above the vapor saturation temperature by adjusting 
. the power to the heating elements around the vapor tubes and the converging unit. 

The remaining vapor from the test section was condensed in the auxiliary condenser. 

{ The condensate from both the test section and the auxiliary condenser was 

continuously pumped back into the boiler by a variable speed micro-pump. The 
steady vapor flow rate was maintained by adjusting the power supply to the heating 
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elements in the boiler. 

A desired average condensing surface temperature was achieved by adjusting the 
flow rate and the temperature of cooling water to the cooling channel. After this 
adjustment, the surface temperature along the test section was with a ±1.25 °C of the 
average value and the total temperature difference of the coolant across the whole 
cooling channel was maintained at about 5 °C. The minimum temperature difference 
of coolant across the whole cooling channel was 1.5 °C to reduce the uncertainty in 
measuring the temperature difference of the coolant. 

The test section was visually inspected at regular intervals. Much attention was 
given to the inlet of the test section to ensure that there was no condensate entering 
the test section from the vapor inlet line. However, it was possible that some portion 
of vapor condensed in the inlet vapor line during the initial warming period of the 
system. If this occurred, such condensate was drained directly to the main liquid tank 
from the bottom of the converging unit through a valve in the tube connecting the 
tank and unit. 

Steady state of the experiments was assumed when the changes in the vapor flow 
rate, the condensate flow rate from the test section, the condensing surface 
temperature, and the total temperature difference of cooling water across the whole 
cooling channel were within ±1.5 %. Generally, it required 1 to \Vz hours to reach 
steady state. 

Measurements 

The experimental measurements were made after steady state was reached in a 
test run. The time for recording the data usually took about 15 minutes. During this 
time, the data from the venturi meters, the turbine meters, and the vapor 
thermocouples were taken repeatedly by the data acquisition system at 1 minute 
intervals and the measurements of the rotameter and the gauge glass for the 
condensate were repeated three times. The data was discarded if there was a large 
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discrepancy between the measurements collected at different times during the test run. 
Measurements were re-conducted one more time after 10—15 minutes. 

To check for repeatability of the results, experiment was repeated for the same 
values of the variable on different days. Due to the nature of the experiment, the 
inlet vapor velocities were within ±1 % and AT within ±1%; the heat transfer rates 
were within ±4 %. 

4.8.2. Procedure of Air-liquid Flow Experiments 

Water-glycerin mixture with a desired viscosity was prepared by mixing distilled 
water with glycerin in the liquid loop. The density and the viscosity of the liquid 
were measured. 

The liquid flow to the test section was started by opening a metering valve in the 
liquid tube. By adjusting the valve, a smooth liquid flow with an even distribution of 
a liquid film along the bottom plate of the test section was obtained. After a certain 
amount of liquid was accumulated in the collecting tank (Fig. 4.10), the liquid pump 
was turned on to pump the liquid back to the gravity tank. The flow rate of the liquid 
to the gravity tank was adjusted by regulating a valve at the outlet of the liquid pump. 
The adjustment of this liquid flow was to maintain a constant liquid level in the 
gravity tank, ensuring a constant liquid flow to the test section. The air flow into the 
test section was started by opening on the valve in the air line. The air flow rate was 
slowly increased until interfacial waves appeared on the liquid surface in the test 
section. Then the air flow rate, the liquid flow rate, and the thickness of the liquid 
film on the bottom plate in the test section were measured. Those measurements 
including the measurement of the viscosity and the density of the liquid were repeated 
twice to check the repeatability of the measurements. 
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4.9. Data Processing 

4.9.1. Condensate Film Thickness 

The condensate film thickness was determined by measuring the time for the 
ultrasonic signal to pass through the liquid film and relating it to the sonic velocity of 
the liquid. The sonic velocity of the liquid is density dependent and was determined 
from an equation suggested by Reid and Sherwood (1958) 


where c m is a constant involving the molecular weight of the fluid. The values of c m 
for R- 1 1 3 and FC-72 were found during calibration of the ultrasonic transducer 
(Appendix D). For R-113, c m = 1.97 x 10' 7 m‘°/kg 3 s and for FC-72, 
c m = 1.077 x lO' 7 m 10 /kg 3 s. 

Referring to Fig. 4.21, an infinitesimal 
film thickness (dy) of the condensate film, is 
expressed as 

dy = a f dr (4-2) 


Condensate Film 


f - - — 1 */: - _ 40 


Substituting Eq. 4-1 into Eq.4-2 yields 

dy = c m p, 3 dr 


Fig. 4.21 Condensate film 


(4-3) 


If the liquid density (p,), is assumed to be a linear function of the liquid temperature 

(T,), 

/>, = p 0 (l + aT,) ( 4 - 4 ) 


Substituting Eq. 4-4 into Eq. 4-3 and assuming a linear temperature profile in the 
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condensate film 


T = a + by 


dy = c m />o( 1 + aa + bay) 3 dr 


The condensate film thickness is found by integrating Eq. 4-5 

r s dy 


0 (1 + aa + bay) 


' = c Pq f 

+ bav) 3 m U 


which gives 


5 = 


3 ^o 


2c m (l + mi-«T s ) 3 p: 
2+0 


( 4 - 5 ) 


( 4 - 6 ) 


where 


0 = 


aAT 

1+aT 


An alternate way to determine the film thickness is to use a mean liquid density 
(pf in Eq. 4-3 to yield 

5 - VX l4 ' 71 

It was found that the difference between the film thickness evaluated from Eq. 4- 
7 and that from Eq. 4-6 was less than +0.1 % (See sample calculations in Appendix 
F). Eq. 4-7 was used in the data processing. 

4.9.2. Heat Transfer Rates and Heat Transfer Coefficients 
Total Heat Transfer Rate and Average Heat Transfer Coefficients 

The total heat transfer rate from the vapor to the whole condensing surface was 
measured using two different methods: (1) heat balance in the cooling channel; 
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(2) heat balance in the vapor duct of the test section. 

r 

(1) Heat Balance in the Cooling Channel 

The total heat transfer rate (q) of cooling water through the whole cooling 

i 

channel is given by 

q = m c AT (4-8) 

*■ w pw w 


(2) Heat Balance in the Vapor Duct of the Test Section 


m t , h g | 
Vapor in 


Heat Insulation 



(rn t - m, ), h go 


Vapor out 


I 

i 


Fig. 4.22 Illustration of heat balance in the vapor duct of the test section 

If the test section is taken as the control volume (Fig. 4.22), the total enthalpy of 
the vapor flowing into the control volume is m,!^ and the enthalpies of the vapor and 
the condensate flowing out of the control volume is [(iii, - m lt )h go + m, t h l0 ]. Applying 
the First Law of the thermodynamics to the control volume 

m,h gj = q + (m t -m lt )h go + m lt h lo (4-9) 

where 

h gl — specific enthalpy of the vapor at the inlet of the test section; 
h go — specific enthalpy of the vapor at the outlet of the test section; 
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h l0 — specific enthalpy of the condensate at the outlet of the test section; 
rh, — mass flow rate of the vapor at the inlet of the test section; 
rh„ — mass flow rate of the condensate at the outlet of the test section. 
The specific enthalpys of the vapor and the condensate were evaluated as 


t 

h . 

= \ 

+ c (T . 

pg v g' 

- T ) 

sat-' 


h go 

= \ 

+ 

O 

3 

3 

O 

- T ) 

sat'' 


K 

= h, 

+ Cp,( T ,o 

- T ) 

sat' 


where 

T gj — vapor temperature at the inlet of the test section; 

T g0 — vapor temperature at the outlet of the test section; 

T l0 — condensate temperature at the outlet of the test section. 

Rearranging Eq. 4-9 

q = m t (h gi - h go ) + m u (h go - h lo ) (4-10) 

The heat transfer rate computed by the two methods were usually within ±5 %. 
If the difference exceeded ±10 %, the data were discarded and the experiment 
repeated. 

As the sectional heat transfer rate could be computed only by Eq. 4-8, therefore 
for consistency, Eq. 4-8 was used to compute the total heat transfer also with Eq. 4- 
10 providing a check on the results given by Eq. 4-8. 

Based on the areas of the condensing surface, three average heat transfer 
coefficients were determined 
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(4-1!) 


(4-12) 


2A (T , - T ,) 

s v sat s2 7 


t 



AAT 


(4-13) 


where 

q, — total heat transfer rate on section 1 (q, = m w c pw AT wl ) 

q 2 — total heat transfer rate on sections 1 and 2 (q^ = m w c pw AT w2 ) 

q — total heat transfer rate on the entire test section (q = m w c pw AT w ) 

T sl — arithmetic mean condensing surface temperature of section 1; 

T s2 — arithmetic mean condensing surface temperature of sections 1 and 2; 

AT W , — temperature difference of the coolant through section 1; 

AT w2 — temperature difference of the coolant through sections 1 and 2. 

Sectional Heat Transfer Rates and Heat Transfer Coefficients 

The sectional heat transfer rate (q sj , i = 1,2,3) was calculated from 


where AT WI was the temperature difference of the coolant through section i 

Knowing q sj , the sectional heat transfer coefficient (h sj , i = 1,2,3) was evaluated 

fi = ^ (4-15) 

s ‘ A AT 

S SI 

Local Heat Transfer Rate and Heat Transfer Coefficient 

As in most analyses on film condensation, when the condensate film thickness is 
small (6 < 1 mm), the temperature distribution within the condensate film is assumed 


q . = m c AT 

*si w pw wi 


(4-14) 
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t 


to be linear. Based on this assumption, the local heat transfer rate is determined from 

k AT 

// = ( 4 - 16 ) 

4 5 

where AT, is the difference between the saturation temperature and the condensing 
surface temperature at x (distance from the leading edge). 

Assuming a power law relation between 5 and x 

5 = C,X Cl (4-17) 


where c, and c 2 are the constants determined by a regression analysis on the measured 
values of the film thickness. Substituting Eq. 4-17 into Eq. 4-16 


q 


// 


k.AT _ c 
_1 — -x 1 


(4-18) 


The local heat transfer coefficient is defined as 


AT x 


Substitution of Eq. 4-18 into the above equation 



(4-19) 


f 


Using the expression for the local heat transfer rate (Eq. 4-18), the total average 
heat transfer coefficient and the sectional heat transfer coefficients were deduced. 

For the total average eat transfer coefficient 

. . w J.Vdx (4 . 20 . 

' AAT 
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Substitution of Eq. 4-18 into Eq. 4-20 (assuming AT, = constant) 



( 4 - 21 ) 


I 


For the sectional heat transfer coefficients 






Substitution of Eq. 4-18 into the above equation (assuming AT, = AT si ) 



( 4 - 22 ) 


where 


L; —length from the leading edge of the condensing surface to the end of 
section i; 

Li.) —length from the leading edge of the condensing surface to the beginning of 
section i. 

The total average heat transfer coefficient and the sectional heat transfer 
coefficients calculated by Eqs. 4-20 and 4-21 were used to compare the values 
determined from (Eqs. 4-13 and 4-15) the method of heat balance in the cooling 
channel. From these comparisons, the technique using the film thickness 
measurement to obtain the local heat transfer coefficient (Eq. 4-19) was validated. 

4.9.3. Liquid Velocity 

In studying the interfacial wave initiation, a knowledge of the variation of the 
local, average liquid velocity along the test section was needed. The local, average 
liquid velocities are found from 
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v,(x) = 


5W 


( 4 - 23 ) 


In the adiabatic air-liquid flow, the liquid flow rate (rh,) was a constant, while in the 
condensation, m, was a function of x due to the condensation of vapor along the 
condensing surface 

*. = w dx 


where m„ was the condensate rate per unit area at x. From a heat balance, 


I 



dx = 



or 


m, = 


wj q"dx 


( 4 - 24 ) 


Substitution of Eq. 4-24 into Eq. 4-23 

v,(x) 



dx 


( 4 - 25 ) 


Substitution of Eqs. 4-17 for 5 and 4-18 for q" into Eq. 4-24 

v,(x) = -A_ f AT x x~ c ’ dx 

c.V'V 0 


After integration of the above equation, and assuming AT X * AT 
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v,(x) 


k.AT 


I -2c, 


CiO-c 2 )h fg 


( 4 - 26 ) 

r 


4.9.4. Vapor and Air Velocities 

Similar to the liquid velocity, the variation of the local, average velocities of the 
vapor and the air (in the air-liquid experiment) along the test section was accounted 
for in studying the interfacial wave initiation. The velocities of both the vapor and 
the air are defined as 


v g 00 


(H - 5)W 


( 4 - 27 ) 


» 


where m g was the mass flow rate of the vapor or the air at x. For the air, m g was 
constant, while for the vapor, m g varied along the condensing surface and it was 
evaluated from 

m g = m t - m, 


Substituting Eq. 4-24 into the above equation 


j„'q" dx 

h 


( 4 - 28 ) 


If H > 5, Eq. 4-27 simplifies to 

v g (x) 


HW 


V 

( 4 - 29 ) 


r 
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5.1. Introduction 

An analytical model for condensation of a vapor flowing parallel to a bottom 
condensing surface in a horizontal rectangular duct was formulated by employing the 
conservation laws. Mass and energy balance equations were obtained using the 
integral method. The resulting equations were numerically solved to predict 
condensate film thickness, condensate velocity, interfacial shear stress, and heat 
transfer coefficients. 

5.2. Physical Description of the Problem 


Adiabatic Wall 


Vapor 


V g (x). T gs 


"V ' ' V' 

y Condensate Rim 


Vapor Boundary Layer H 




Condensing Suface, T s - Constant 


Fig. 5.1 Physical model and coordinate system 


» 


The physical situation considered in the analysis is illustrated in Fig. 5.1. Pure 
vapor flows into a horizontal rectangular duct along the x direction with uniform inlet 
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velocity (v gj ) and saturation temperature (T sat ). The bottom plate of the duct is 
maintained at a uniform temperature (T s ) while the other three sides of the duct are 

» 

perfectly insulated. Since T sal > T s , vapor condenses on the bottom plate. The 
condensate film on the plate flows in the x-direction due to the action of the shear 
stress at the vapor-condensate interface and the momentum of the vapor. The flow of 
vapor is turbulent throughout the duct. 

5.3. Formulation 

5.3.1. Assumptions 

1 . Vapor enters the duct with a uniform velocity 

2. Vapor flows as a boundary layer flow over the condensate 

3. The interfacial thermal resistance is negligible, so that the condensate 
temperature at the interface is same as the saturation temperature 

4. The change in the static pressure of vapor along the duct is small, so that the 
saturation temperature along the duct is uniform 

5. The condensate film thickness is very thin and the temperature profile in the 
film layer is linear 

6. The height of the duct is much greater than the condensate film thickness 

7. The width of the duct is much greater than the condensate film thickness, so 
that the condensate flow is two-dimensional 

8. The properties are constant 

9. No non-condensable gases exist in the duct 

% 

5.3.2. Mass Conservation 

I 

Mass Conservation of Condensate Film 

Mass conservation to the control volume of condensate film (Fig. 5.2) yields 
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* 


* 


m x Ax 


m A* » I ! m, (x+Ax) 

m i' x ' 6 Condensate Film — ! — 




Fig. 5.2 Control volume of condensate 

m,(x+Ax) - m,(x) - m x WAx = 0 


Rearranging the equation and letting Ax -* 0, 

1 dm,(x) 


m. = 


W dx 


( 5 - 1 ) 


in, is expressed as 


d 1 ! = w /°i| o " v i( x iy)dy 


( 5 - 2 ) 


Substituting Eq. 5-2 into Eq. 5-1, an expression for the condensate rate is given 

d r f 5 


= p,^[| o v,( x -y) d y] 


( 5 - 3 ) 


Mass Conservation of Vapor 

Fig. 5.3 shows a control volume of vapor in the duct. Mass conservation to the 
control volume gives 
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m gW 


m,(x) 


y 


m g (x+Ax) 


rh ^x+Ax) 

a Condensing Surface : — - 

\/ //////////// // /////// /////////////, i 

AX 


Adiabatic Wall 

y//////s//£ Y//fr///s///z/s///////sz 

Vapor 

m.Ax 


Fig. 5.3 Control volume of vapor 

m g (x+Ax) + rh x WAx - m g (x) = 0 


Rearranging the equation and letting Ax 

dn^ 

dx 


-0, 

= -m x W 


The mass flow of vapor, rh g , at x is expressed as 

m g = p g (H-5)Wv g (x) 


Since H > 5 (assumption 6), the above equation is simplified as 

m g = p g HWv g (x) 


Substitution of Eq. 5-5 into Eq. 5-4 yields 


1 f * 


V (x) = v - f 

8 81 Hp g J o 


m dx 


9 


(5-4) 


(5-5) 


(5-6) 
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5.3.3. Energy Conservation 


4 


(m x Ax)h g 



Fig. 5.4 Control volume for energy conservation 


Neglecting changes in kinetic and potential energies of the condensate and viscous 
dissipation, the energy balance applied to the control volume shown in Fig. 5.4 is 

H,(x) + [m x WAx]h g =q // WAx + H,(x + Ax) ( 5 ‘ 7 > 


where 

H, — rate of enthalpy flow across the condensate film; 

q" — rate of heat transfer per unit area of the condensing surface. 

Rearranging Eq. 5-7 and letting Ax -» 0 yields 

q " (x) = — -H.+iii h (5-8) 

M v ' W Hv 1 * 8 


H, is expressed as 
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H 1 (x)=p 1 wJ o 5 v 1 (x,y)h,dy 


( 5 - 9 ) 


Substitution of Eqs. 5-9 and 5-3 into Eq. 5-8 yields 

d r * 


q" (x) = p,h fg _ j o V,(x,y)dy 


( 5 - 10 ) 


Based on the assumption of a linear temperature profile in the condensate, q" is 
given by 

k,(T a , - T s ) 


q (x) = 


( 5 - 11 ) 


Substitution of Eq. 5-11 into Eq. 5-10 yields 

d r r « . A, 


-f-[[ v(x,y)dy]=- 
dx Jo 


( 5 - 12 ) 


where 


k,AT 

P> h fg 


( 5 - 13 ) 


5.3.4. Condensate Velocity 

It has been well established that if the condensate film is thin the condensate 
velocity (Vl) in the film can be assumed to be linear 

v,(x,y) = c(x)y (5 ‘ ,4 > * 

The condition that t, = /u, (dv/dy) at y = 5 leads to c(x) = Thus ^ 
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( 5 - 15 ) 


where t, is the interfacia] shear stress varying in the x direction. 

Different hypotheses have been proposed for computing r,. Some researchers [for 
example, Rosson and Myer (1965), Soliman et al. (1968), and Bae et al. (1971)] 
assumed that r x in condensation was same as that in an adiabatic two-phase flow; 
therefore, the method applied in an adiabatic flow was used to determine the t x in 
condensation. Some other researchers [for example, Shekriladze and Gomelauri 
(1966), Linehan et al. (1969), Jensen and Yuen (1982)] assumed that t { in 
condensation resulted from a combination of an adiabatic shear stress exerted on the 
interface (without condensation) and a momentum shear stress due to the mass of 
vapor condensing on the film layer; when the condensate rate was high, the influence 
of the momentum was possibly the dominant contribution to the total value of r ; . It 
appears that the latter assumption is more appropriate to the situation of condensation 
and is thus used in this analysis. Accordingly, r x is expressed as 


where v is the average vapor velocity relative to the interfacial velocity (v^, i.e. 


and f c is a local adiabatic friction factor and assumed to be a function of Re x (vx/ v z ) 



2 


+ m x v 


( 5 - 16 ) 



( 5 - 17 ) 


f c = CfRe x n ' 


( 5 - 18 ) 


A discussion of the choice of the constants (q and n f ) is presented in Chapter 6. 
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5.3.5. Interfacial Velocity 

The mass rate of flow of the condensate at x is given by 

m, = p 1 v,(x)W5 = p,W JJv,(x,y)dy < 5 - |9 > 

_ 9 

where v, is the average liquid velocity varying only in the direction of x. After 

rearrangement of Eq. 5-19, 

v, = V,(x,y)dy (5-20) 


Employing Eq. 5-14 


v„ = 2v,(x) 


( 5 - 21 ) 


Using Eq. 5-1 and 5-19, it is found that 


v 


i 



Since 


m x dx = 


m f -idx 

h„ Jo 6 


( 5 - 22 ) 


( 5 - 23 ) 


Employing Eq. 5-23, substituting Eq. 5-22 into 5-21, the expression for the interfacial 
velocity is 


v 


li 


2 A, r *dx 

“ITJo T 


( 5 - 24 ) 


where A, is defined by Eq. 5-13. 


i 
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5.3.6. Differential Equation for Condensate Film Thickness 
Substitution of Eq. 5-15 into Eq. 5-12 leads to 

-Itrfl - 2 J^1 
dx ' 6 


or 


dr, + « , 2V, 

dx 'dx 5 2 


Substituting Eq. 5-18 into Eq. 5-16 (with considering Re x = vxh e ) 

r = f^iv 2 * n -x n ' + m v 
2 f' 

g 


Using Eqs. 5-3 and 5-12, rii x is expressed as 


. A iP ( 

m. = 

* 8 


Substituting Eq. 5-27 into Eq. 5-26 


T i 



i 

— _v 

8 


k 


Differentiating r, with respect to x 


dr, 

dx 


= [(2 + n f )A 2 xV* n ' + 


A,P| dv 

~5” J dx 


A ,P,v d5 
5 2 dx 




2 +n f 


where 


( 5 - 25 ) 


( 5 - 26 ) 


( 5 - 27 ) 


( 5 - 28 ) 


( 5 - 29 ) 
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A, = 


c _fz 


( 5 - 30 ) 


Substituting Eqs. 5-6 and 5-24 into Eq. 5-17 for v and using Eq. 5-27 


v = v 


AjPi r *dx _ 2 A, j- x<j x 

H^TJoT TjoT 


( 5 - 31 ) 


Differentiation of v with respect to x yields 

dv _ 2 Aj j- x(j x d6 _ A,p, ^ 2A j 
dx “ + ~b~’f 


Combining Eqs. 5-25, 5-26, 5-29, and 5-31, an ordinary differential equation for 
6 is obtained 

d5 x + [a(x)5+A,p 1 ](|5+T)) - 7 (x)5 3 (J _ 33) 

dx e[a(x)5+A,pJ + 25[0(x)5+/3(x)] - /?(x)5 

where 


a(x) = (2 + n f )A 2 x n 'v Un ', /3(x) = p,A,v 

7(x) n ( A 2 x n '- , v 2 *"', e(x) = 2A,p, f ^ 

J 0 0 


I = 



V = 2a, 


0(x) = A 2 v 2 * n 'x n ', x = 2A,p, 


4 
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Define dimensionless variables 


5 = 5*L, x = x *L 


( 5 - 34 ) 


Substituting Eq. 5-34 into Eq. 5-33, a non-dimensional differential equation for 
condensate film thickness is given by 


dx * e[a(x * )L5 * +A,pJ + 2L5 * [0(x * )LS * +0(x * )] - 0(x * )L6 * 

5.3.7. Heat Transfer Rate and Heat Transfer Coefficients 
Total Heat Transfer Rate and Average Heat Transfer Coefficients 

With AT constant, the total heat transfer rate from vapor to the entire condensing 
surface is given by 


The average heat transfer coefficient, h si , over a certain area of condensing 
surface, A s , is defined as 


d5* 


k + [ctfx^LS’+AjpJtfLS’+iO - y(x*)L 3 S*’ 


( 5 - 35 ) 


q = Wk,AT 


L dx 
o 6 


( 5 - 36 ) 


( 5 - 37 ) 


A AT 


where 


q,, = Wk,AT| 



and 


A s = W(L, - L ,) 


where the meanings of L s and Lj., is illustrated in Fig. 5.5. 
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Leading Edge Condensing Surface 



Fig. 5.5 Illustration of L, and L,., 


Rearranging Eq. 5-37 



L 


k, 

- Li-, 


l, dx 

k,T 


( 5 - 38 ) 


If L„, = 0 and = L, the expression for the total heat transfer coefficient (hj over 
the entire condensing surface is given by 



( 5 - 39 ) 


Local Heat Transfer Coefficient 

The local heat transfer coefficient at x is defined as 

h - q" 

* (T -T) 

^ sat s/ 

Substituting Eq. 5-1 1 into the above equation yields 


( 5 - 40 ) 
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5.4. Solution Procedure 

% Eq. 5-35 was solved numerically. The grid arrangement for the numerical 

* 

computation is shown in Fig. 5.6. The dimensionless length of the condensing 
^ surface is taken as one (L’=l). The total number of the grids was 600, of which 300 

grids were evenly spaced in the leading section of the condensing surface (from x‘ = 

0 to x* = 0.2) as 5 increases rapidly at the leading edge of the condensing surface; 
the rest of grids were evenly distributed along the remaining length of the condensing 
surface. The total number of grids (600) was determined from a pre-computation of 
Eq. 5-35. It was found that when the total grid number was greater than 600, the 
relative change in b was less than 0.05 %. 

In the calculations, the vapor properties were evaluated at T sat ; the liquid 
properties were determined at the arithmetic average of the saturation temperature and 
the surface temperature [(T sat + T s )/2], The effect of subcooling was considered by 
modifying h fg using the following equation (Rohsenow, 1970) 


f AJCl » 0.0007 
Ax2- 0.0027 


jv. Vapor 

, > 

Condensate film 


1 

1 
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Fig. 5.6 Arrangement of grids 
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Fig. 5.7 Flow chart of computing program 
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K = h f g + °.68c pl (T sat -T ) 

Eq. 5-35 was first solved to obtain 8. Then, v u , v„ v g , t x , h x , h sl , TT, were 
computed. In the computations, the values of c f and n f in Eq. 5-18 were determined 
based on a modified Schlichting friction equation (details are described in 
Chapter 6). The computational procedure (illustrated in Fig. 5.7) is as follows. 

1 . Input m„ T sa „ AT 

2. Find the necessary thermodynamic properties from the subroutine program 
for the properties 

3. Calculate A! (Eq. 5-13), A 2 (Eq. 5-30), and v gi 

4. In the initial calculation of iteration, assume that v H is 0 and v is constant 
(assuming v = v gi ) 

5. Solve Eq. 5-35 for 5* using the 4 th order Runge Kutta method 

6. Use Eq. 5-31 to calculate v including v g and v (j , employing the values of 8 
obtained in step 5 

7. Repeat step 5 until x* = 1 

8. Calculate the total heat transfer rate using Eq. 5-36 

9. Repeat steps 5 -* 8 once 

10. If the relative difference of the newly calculated heat transfer rate with the 
previously obtained value is less than 10' 3 (e = |q„ ew - q oW |/qnew < 0.001), 
halt the iteration; otherwise repeat steps 5 -* 8 

1 1 . Calculate m x , v u , V„ v g , r { using Eqs. 5-27, 5-24, 5-20, 5-6, and 5-28 

12. Calculate F t , h^, h x using Eqs. 5-39, 5-38, and 5-40 
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6.1. Introduction 

Condensation heat transfer with a vapor flow inside a horizontal rectangular duct, 
with condensation only on the bottom plate, was experimentally studied. 
Measurements were obtained both with a smooth condensate flow and a wavy 
condensate flow. The initiation of interfacial waves was investigated with 
condensation and with an adiabatic air- liquid flow. An analytical model was 
developed to simulate the condensation process. Typical results of the experimental 
measurements of the heat transfer coefficients, predictions of the analytical model, 
and the correlations relating the average heat transfer coefficient to the variables of 
significance, are discussed in this chapter. All the experimental results are given in 
Appendix G. 

6.2. Ranges of Experimental Variables 
6.2.1. Condensation Heat Transfer 

There were 118 valid runs of experimental data with R-113 and 75 runs with 
FC-72; interfacial waves were observed in 76 runs with R-113 and 56 runs with FC- 
72. The two main variables in the experiment were the inlet vapor velocity to the test 
section (v gi ) and the difference between the saturation vapor temperature and the 
condensing surface temperature (AT). In the test section, since the vapor pressure 
changed in the range approximately from 2 Pa to 10 Pa, the saturation temperature of 
vapor (T sat ), is assumed constant. Table 6.1 lists the ranges of the values of the inlet 
vapor velocity, (v gi ), the temperature difference (AT), and the inlet vapor Reynolds 
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number (Re gL = v gl L/»> g ). 

Table 6.1 Ranges of variables in the experiment of condensation heat transfer 


Condensing Fluids 

v^, m/s 

u 

o 

Re gL 

R-l 13 

0.31 - 4.39 

10-40 

210,000 - 2,700,000 

FC-72 

0.36 - 3.10 

10 - 50 

380.000 - 3,322,000 


6.2.2. Wave Initiation 

There were 46 runs of data with the adiabatic air-liquid flow and 8 runs of data 
with the condensation (R-l 13 and FC-72). For wave initiation with the air-liquid 
flow, the main variables were: air velocity (v gj ), mass flow rate of liquid (m,), and 
viscosity of liquid (/i^; for wave initiation with condensation, the main variables were 
the inlet vapor velocity (v gi ) and the temperature difference (AT). The ranges of the 
variables in the experiments are given in Table 6.2. 


Table 6.2 Ranges of variables in the experiment of interfacial wave initiation 


Fluids 

Ml X 10 6 
kg/m s 

m, X 10 2 
kg/s 

M. X 10 6 
kg/m s 

m/s 

R-l 13 (condensation) 

509- 517 

— 

11 

1.3-1. 6 

FC-72 (condensation) 

439 — 461 

— 

12 

0.9-1. 1 

Glycerine-water (adiabatic) 

1000 — 7550 

0.1 — 1.2 

18 

0.5-6. 0 


Since the experiments with the adiabatic air-liquid flow were performed at 
atmospheric pressure and room temperature ( —20 °C), the viscosity of air was 
assumed constant. Also, due to the insignificant change of the saturation temperature 
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of vapor in the experiment with condensation, the viscosity of vapor was 
approximately constant. However, the viscosity of the liquid (especially, the liquid 

v viscosity of glycerine-water) varied in the experiments (Table 6.2). 

* 

% 6.3. Observation of Flow Patterns 

The condensate flow patterns along the condensing surface were observed 
visually. The condensate Reynolds number (Re^ was less than 450 so that the flow of 
condensate is assumed to be laminar. It was observed that the vapor-condensate 
interface was either smooth or wavy, depending on the magnitude of the inlet vapor 
Reynolds number (Re gL ) and the temperature difference (AT). Table 6.3 lists the 
value of Re gL (critical Reynolds number) at which the transition of smooth interface to 
a wavy interface occurred (wave initiation). 


Table 6.3 Inlet vapor Reynolds number at which interfacial waves appear— condensation 


AT 

°C 

R-113 

Re gL 

Vp m/s 

FC-72 

Re gL 

Vp m/s 

10 

1,007,093 

1.59 

1,225,835 

1.13 

20 

886,762 

1.4 

1,108,890 

1.09 

30 

804,421 

1.27 

918,545 

0.90 

40 

802,312 

1.27 

883,375 

0.84 

50 

— 

— 

851,826 

0.81 


In the wavy flow regime, both two dimensional (2-D) and three dimensional (3- 
I D) waves can appear. Fig. 6. 1 shows photographs of typical 2-D waves and 3-D 

waves. It is seen that the 2-D waves have a clear wave length, while the 3-D waves 
are rather randomly shaped (a pebbled appearance). The 3-D waves appeared after 
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Fig. 6.1 Photographs of 2-D waves and 3-D waves 

a. AT = 11.1 °C; v gi = 1.85 m/s; Re gL = 1,166,472; R-113 

b. AT = 10.2 °C; v gj = 2.95 m/s; Re flL = 1,593,353; R-113 


4 


6.4 


6. Results and Discussion 


Leading edge of the condensing surface 



77 


73 79 

i . | 




31 » » 54 55 35~S /' 7 i« 5» 50 51 52 5: 





Trailing edge of the condensing surface 


Fig. 6.2 Photograph of condensate flow patterns along condensing surface 
FC-72; AT = 9.3 °C; v,, = 1.5 m/s; Re al = 1,593.353 
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the 2-D waves along the condensing surface. When Re gL was high, 3-D waves 
occurred all over the condensing surface. Table 6.4 gives the Re gL at which the 3-D 
waves covered all the whole condensing surface. Before Re gL reached the values 
listed in Table 6.4 at certain values of AT, all three flow patterns (smooth, wavy with 
2-D waves, and wavy with 3-D waves) usually appeared on the condensing surface. 

It is seen in Fig. 6.2 that the surface of condensate was smooth at the beginning of 
the condensing surface, and gradually became wavy with 2-D waves, which changed 
into 3-D waves towards the end of the condensing surface. 

Table 6.4 Conditions at which 3-D waves covering the whole condensing surface 


AT 

°C 

R-113 

Re gL 

v^ m/s 

FC-72 

R V 

Vp m/s 

10 

1,358,584 

2.20 

1,867,031 

1.80 

20 

1,121,217 

1.78 

1,700,628 

1.60 

30 

1,120,653 

1.73 

1,471,876 

1.38 

40 

976,646 

1.52 

1,471,163 

1.34 

50 

— 

— 

1,345,425 

1.30 


6.4. Temperature Distributions in the Test Section 

Fig. 6.3a gives the typical distributions of vapor temperature (T g ), condensing 
surface temperature (T s ), and coolant temperature (T c) in the test section when R- 113 
was used as the condensing fluid. Fig. 6.3b gives these measurements when FC-72 
was used. It is seen in both the figures that the vapor temperature is nearly uniform, 
except near the inlet of the test section where it is slightly higher due to the heating of 
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Fig. 6.3 Various temperature distributions in the test section 

a. R-113; AT = 30 AT; v* = 2.2 m/s 

b. FC-72; AT = 30 °C; v* = 1 .2 m/s 
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the vapor along the inlet tube. The maximum difference between the average vapor 
temperature and the saturation temperature was within 5 °C. The temperature of the 
condensing surface is seen to rise slightly along the condensing surface, but the 
overall variation of the surface temperature from its average value was less than ± 
1.25 °C. 

6.5. Condensate Film Thickness 

The condensate film thickness of condensate with a smooth interface and a wavy 
interface with small amplitudes and speeds was measured. When more waves 
appeared on the interface and the interfacial wave speed was high, the measurement 
of the film thickness was not possible because of the difficulty in distinguishing the 
interface due to its rapid oscillation. 

Fig. 6.4 illustrates the distributions of the film thickness (5), along the length of 
the condensing surface (x), using R-113 and FC-72 as the condensing fluids at two 
different values of AT. The film thickness increases rapidly in the leading section of 
the condensing surface (high condensate rate); however, the rate of the increase of the 
film thickness decreases with distance from the leading edge of the condensing 
surface; with increasing v gi , the film thickness decreases, but near the trailing edge of 
the condensing surface, especially for the cases with large AT (Fig. 6.4 b), the effect 
of v g , on the film thickness is not significant. 

6.6. Heat Transfer Rates and Heat Transfer Coefficients 

6.6.1. Total Heat Transfer Rates and Average Heat Transfer Coefficients 

The variations of total heat transfer rate (q) with the inlet vapor velocity (v gi ) for 
four different values of AT and the two condensing fluids are plotted in Fig. 6.5. 

This figure shows that the heat transfer rate increases with increasing inlet vapor 
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a. R-113; AT = 10 °C 


b. R-113; AT = 40 °C 
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c. FC-72; AT = 10 °C 


d. FC-72; AT = 40 °C 


Fig. 6.4 Distributions of condensate film thickness along the condensing surface 
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Fig. 6.5 Total heat transfer rates 
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velocity. The rate of the increase is higher when the inlet vapor velocity is greater 

than 1.6 m/s. (Note: when the inlet vapor velocity is higher than approximate 1.6 

v m/s, interfacial waves developed on the condensate surface). However, the influence 

♦ 

of the inlet vapor velocity on the total heat transfer rate at AT = 10 °C is not as 
* significant as it is on the total heat transfer rates at other values of ATs. 

The average heat transfer coefficients are calculated using Eq. 4-13. The 
variation of the average heat transfer coefficient with inlet vapor velocity is shown in 
Fig. 6.6. In the figure, h, denotes the value of the average heat transfer coefficient 
over the section 1; h 2 is the value over section 1 and section 2; h, is the overall 
average heat transfer coefficient over the entire condensing surface. The three 
different average heat transfer coefficients display similar trends: 

— the average heat transfer coefficient increases with an increase in the vapor 
velocity 

— beyond a certain value of the vapor velocity, the average heat transfer 
coefficient increases nearly linearly with the vapor velocity 

The total heat transfer coefficient can also be determined from the measurement 
of the measured film thickness using the equation (Eq. 4-21) 

h,.2£1 

c,(l-c 2 ) 

where c, and c 2 are the coefficients in the film thickness correlation (5 = c,x c2 ) and 
* are determined from a regression analysis of the experimental data. Fig. 6.7 shows a 

comparison between the total heat transfer coefficients obtained using the film 
( thickness distribution (h tf ) and those measured using the method of the heat balance of 

coolant through the test section (hj. Using the definition of the average deviation 
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Fig. 6.6 Variation of average heat transfer coefficients with v tf 
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Fig. 6.7 Comparison of heat transfer coefficients measued from the heat balance 
method with the values obtained from the film thickness data 
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the average deviation between the two different values is found to be within 10 %. 
This comparison validates the technique of measuring the local condensate film 
thickness by the ultrasonic transducer. As a check of the validity of the power law 
correlation for the film thickness, two typical correlations (8 = 0.48 lx 045 and 8 = 
0.483x 024 ) are shown in Fig. 6.8. The average deviation between the calculated 
thickness using the correlation and the measured values is less than ±6 % for both 
the correlations. As shown in Appendix G, the range of the values of the exponent 
(c 2 ) in all the correlations is 0.17 to 0.47 for R-113 and 0.15 to 0.50 for FC-72. 

Fig. 6.9 shows the effect of the temperature difference (AT) on the total average 
heat transfer coefficients (hj. In the figure, h, is plotted against the mass velocity 
(G) to compare the trends of the present experimental results with those of Akers and 
et al. (1958). It is seen that: 


— the influence of AT on ¥ t for R-113 is insignificant in the whole range of the 
mass velocity of vapor (G < 33 kg/s m 2 ) 

— the influence of AT on ¥, for FC-72 is also insignificant when G was less 
than or equal to 25 kg/s m 2 

— for FC-72, when G is greater than 25 kg/s m 2 , ¥, with AT = 10 °C is 
smaller than the others 

— for FC-72, when G is greater than 35 kg/s m 2 , h, with AT = 30 °C is less 
than ¥, with 20 °C, but greater than ¥, with 10 °C (The reason for this result 
is not clear and it needs further investigation) 

Akers et al. (1958) measured the total average heat transfer coefficient for R-12 
flowing inside a horizontal condensing tube. Their results also showed that the heat 
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Fig. 6.8 Typical examples of film thickness correlations 

a. R-113; AT = 22.8 °C; v* = 1.8 m/s 

b. FC-72; AT = 30.8 °C; v* = 0.8 m/s 
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c. Akers, Dean, and Crosser (1958) 


Fig. 6.9 Total average heat transfer coefficients as a function of vapor mass velocity to the inlet 

of the test section 
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transfer coefficient with a lower AT is smaller than that with a higher AT when the 
total vapor mass velocity was higher than 25 kg/s m 2 . Abis (1969) also studied the 
effect of AT on h, in a horizontal condensing tube. However, no conclusion about the 
relationship between AT and h t was drawn from his experiments. 

6.6.2. Sectional Heat Transfer Rates and Sectional Heat Transfer Coefficients 

The condensing surface was divided into three equal sections; the heat transfer 
from the vapor to each section is the sectional heat transfer rate (q 3 ). The sectional 
heat transfer rate is obtained from an energy balance on the coolant in each section. 
Fig. 6.10 illustrates the variation of the three sectional heat transfer rates with the 
condensing surface. It is seen that the heat transfer rate at section 1 is higher than the 
heat transfer rates at sections 2 and 3, while the change of heat transfer rate from 
section 2 to section 3 is quite small. 

The sectional heat transfer coefficient is an average heat transfer coefficient in the 
section and is evaluated using Eq. 4-15. Two typical distributions of the sectional 
heat transfer coefficients along the length of the condensing surface are shown in Fig. 
6.11. Similar to the trend of the sectional heat transfer rates, the sectional heat 
transfer coefficients decrease from section to section along the length of the 
condensing surface. 

Fig. 6.12 compares the sectional heat transfer coefficients of R-113 with those of 
FC-72 at approximately the same values of AT and m, (one plot for lower m, and the 
other for higher m,). It is seen that when m, is lower, the sectional heat transfer 
coefficient of R-113 is larger than that of FC-72 (figure a). However, when m t is 
higher, the reverse happens except in the first section (figure b). 

The phenomenon shown in Fig. 6.12a may be explained by comparing the film 
thickness of R-l 13 with that of FC-72 at a lower m t . At approximately the same m t 
and AT, the value of the film thickness of R-113 is smaller than that of FC-72 
(Fig. 6.13a) and this is more obvious when AT is higher (Fig. 6.13b). Studying the 
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Fig. 6.10 Sectional heat transfer rates 

a. R-113; AT = 21.7 °C; v, = 1.3 m/s 

b. FC-72; AT = 19.4 °C; v* = 2.5 m/s 
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Each point represents the average heat transfer coefficient at a section whose 
midpoint corresponding to the x-value. 
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Fig. 6.11 Sectional heat transfer coefficients 

a. R-113; AT = 10.8 °C; v* = 1.6 m/s 

b. FC-72; AT = 39.1 °C; v* = 2.5 m/s 


Note: Each point represents the average heat transfer coefficient at a section whose 

midpoint corresponding to the x-value. 
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Fig. 6.12 Comparison of sectional heat transfer coefficients of R-1 13 with that of FC-72 


Note: Each point represents the average heat transfer coefficient at a section whose 

midpoint corresponding to the x-value. 
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thermal properties of the condensing fluids (Table 4.1), it is found that the thermal 
conductivity of liquid R-l 13 is about 1.3 times higher that of FC-72. The thinner 
film thickness and higher thermal conductivity of R-l 13 results in a lower thermal 
resistance across the condensate film, and therefore leads to a higher heat transfer 
coefficient. 

At higher m, (implying higher Re gL ), interfacial waves appear on the condensate 
film. The interfacial waves influence the heat transfer coefficients and the extent of 
the effect of the interfacial waves on the heat transfer coefficients may vary for 
different condensing fluids. 

Based on the measured distribution of the film thickness along the condensing 
surface, the sectional heat transfer coefficients can also be calculated using the 
measured film thickness (Eq. 4-22). The sectional heat transfer coefficients obtained 
using the measured film thickness (TT, f ) are compared with those measured using the 
method of the heat balance of the coolant flow (h s ) and the results are shown in 
Fig. 6.14. The average deviation between F sf and ¥ s in section 1 is within ±10 %; in 
section 2 and 3 the average deviation is between ±10% and ±20 %. 

6.6.3. Local Heat Transfer Coefficients 

The local heat transfer coefficients (hj are obtained from the measurement of the 
condensate film thickness along the condensing surface 



where c t and c 2 are the coefficients in the film thickness power correlation (Eq. 4-17) 
Fig. 6.15 shows some typical trends of local heat transfer coefficients (hj varying 
with the axial distance of the condensing surface (x). It is seen that with increasing 
inlet vapor velocity (v gi ), the local heat transfer coefficient increases; however, the 
effect due to the increase in the vapor velocity on the local heat transfer coefficient 


6.21 


6. Results and Discussion 



X, 171 


a 



X, m 


b 

Fig. 6.13 Effect of condensing fluids on film thickness of condensate 

a. AT ■» 20 °C; m t »• 0.008 kg/s 

b. AT - 40 °C; rii, - 0.008 kg/s 
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Fig. 6.14 Comparisons of sectional heat transfer coefficients measured from heat 
balance method with values obtained from film thickness data 
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Fig. 6.15 Variation of local heat transfer coefficients along the condensing surface 
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decreases with increasing axial distance from the leading edge; near the trailing edge 
of the surface, this effect is almost negligible. 

, In experiments of film condensation with steam-air mixture inside a horizontal 

square duct (0.1 m on each side and 1.8 m long), Barry and Corradini (1988) also 
« found that the local heat transfer coefficients decrease along the condensing surface 

and increase when the velocity of the steam-air mixture increases. In their 
experiments, the local heat transfer coefficients were obtained from the local heat flux 
obtained by measuring the temperature gradient across the thickness of the condensing 
wall. 

Azer et al. (1971) also report that the local heat transfer coefficients decrease 
along the condensing surface and increase with an increase in the vapor velocity. 
However, their results show that the local heat transfer coefficients vary nearly 
linearly with the axial distance. Their experiment was conducted inside a 12.7 mm 
I.D. and 2.4 m long horizontal tube using R-12 as the condensing fluid. The local 
heat transfer cot icients were obtained from the local heat transfer flux determined 
from the product of the measured temperature gradient of the coolant, the coolant 
flow rate, and the coolant specific heat. 

Traviss et al. (1973) conducted experiments of condensation inside a horizontal 
tube using R-12 as the condensing fluid. The trends of their local heat transfer 
coefficients with the axial distance and the inlet vapor velocity are very similar to the 
trends of the present experiments. 

6.7. Dimensionless Parameters 

6.7.1. Derivation of Dimensionless Parameters 

j To reduce the numbers of variables in any correlation, dimensionless variables 

are employed. Dimensionless variables can be obtained from differential equations or 
from empirical methods. 
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A differential equation for the film thickness of condensation of vapor flow inside 
a horizontal rectangular duct was derived in Chapter 5 and is reproduced below 

d5* * + [a(x*)LS* + A l/0l ](|L5*^) - 7 (x*)L 3 S*’ 

dx * f [a(x * )L5 * +A,p,] + 2 U * [0(x * )L6 * +0(x * )] - 0(x * )U * 


where 

a(x) = (2 + n f )A 2 x n, v Un ', /3(x) = p,A,v 
7 (x) = n f A 2 x n, ' 1 v 2 *" , > e(x) = 2A iPiJ 0 *^ 


£ = 


Vi 



0(x) = A 2 v 2tn, x n ', 


V = 2A, 
k = 2A,^, 


Rearranging the above equation 


d5* _ *(x*) 
dx* «(x*) 


where 


*(x‘) = 2 + [k,(^)(^)Reix* n '5* + (^-)][(i)(-)5* + l] 

P K Pr, H Pg 


+ k3Re g k L ‘(-^)(-)(^) 2 x* k ’6* 3 
Ja p g m, 


and 
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0(x‘) = kj x '^l[k 7 Re^(^) + (A)]-(^)(^)Re gL 5* 

J° 6 P, Pr , P g Pi 

+ 2Re gL (^)(— )(— )^ * [k^RegL 0 (— )x * n '5* +A] 

Ja p g ft, /x, Pr, 

where lq (i = 1 to 10) represent some constants. Eq. 6-1 implies that the film 
thickness is a function of Pr,/Ja, p,/p g , p,/p g , and Re L . Therefore, the following 
relation is assumed 

5* = f s [(5i, Re gL , (^i), (— ),x*] (6-2) 

Ja g p, p g 

Using Eq. 5-39, the average heat transfer coefficient (TT;) over the condensing 
surface with a length (1^) measured from the leading edge can be calculated 

H. - Jil f *1 , i = 1 to 3 (6-3) 

’ L, Jo 6* 


Defining the Nusselt number (Nuj) 


Nu, = 


k ! 


Eq. 6-3 becomes 


Nu, = 


l,- dx * 

o ~y~ 


( 6 - 4 ) 


* If Eq. 6-4 is integrated by substituting Eq. 6-2 for the dimensionless film 

thickness, the parameter (x*) in Eq. 6-2 will be replaced by L* (L,/L), while the 

l 

parameters [Pr,/Ja, Re gL , (p g /pi), p g /pj will remain in the relation for the Nusselt 
number 
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These parameters in Eq. 6-5, in addition to quality, can also be found in the heat 
transfer correlations of some other researchers [for example: Akers et al. (1958), 

Azer et al. (1971), Shah (1979)]. The quality included in those correlations accounts 
for the change in the local vapor velocity along the duct. 

In Eq. 6-5, both Nu ( and Re gL involve a characteristic length. There are two 
possible choices for the characteristic length used for condensation of vapor flow 
inside a duct: one is the hydraulic diameter (Dh), which is suitable for fully developed 
velocity and temperature profiles; the other is the length of the condensing surface (x 
or L), which is appropriate for a boundary layer problem. In the present study, the 
vapor enters the test section with a uniform velocity with a boundary layer on the 
condensate. It is assumed that the boundary layer thickness is much less than the 
height (25 mm) of the test section. (It has been estimated that the boundary thickness 
is less than 5 mm.) With such a boundary layer, the distance from the leading edge 
is more appropriate to be used as the characteristic length. 

6.7.2. Average Nusselt Numbers and Stanton Numbers 


Average Nusselt Numbers 

The variation of the average Nusselt numbers with the Reynolds number (Re gL ) 
based on the inlet vapor velocity for AT = 10 °C and AT = 40 °C is shown in 
Fig. 6.16, where Nu,, Nu 2 , and Nu t are defined as the following 


NUj 


Hi L , 




The Reynolds number at which the interfacial wave initiation was observed in the 
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experiments is denoted by the dashed lines. Before the initiation of the interfacial 
waves, the Nusselt numbers increase slightly with increasing the Reynolds number. 
However, after the interfacial waves appear, the Nusselt numbers increase rapidly. It 
can, therefore, be concluded that interfacial waves have a very significant effect on 
the condensation heat transfer rate. 

Stanton numbers 

The Stanton number [St, = NUj/(Re gl Pr,)] can be expressed as h/(p|C p ,v gi ). The 
variation of the Stanton numbers with the inlet vapor Reynolds number is shown in 
Fig. 6.17, where St, is based on the average heat transfer coefficient (h,) in section 1; 
St 2 on the average heat transfer coefficient (h 2 ) of section 1 and section 2 taken 
together; and St, on the average heat transfer coefficients (hj over the entire 
condensing surface. All the plots in the figure show that the Stanton numbers have 
the highest value at the leading edge of the condensing surface; they then decrease 
rapidly when the Reynolds number increases; beyond a certain value of the Reynolds 
number (about Re gL > 800.000), the Stanton numbers tend to be a constant. The 
Stanton numbers of both R- 113 and FC-72 closely follow the same trend, varying 
with the Reynolds number and overlapping each other. The collapsing of all the 
experimental measurements opens the possibility of a correlation between St and Re gL . 

6.8. Comparisons with Correlations and Data in Literature 

6.8.1. Comparison of Heat Transfer Correlations 

There are a number of correlations predicting the heat transfer coefficient for 
condensation inside tubes. Some correlations are developed to predict the average 
heat transfer coefficients, while others can be used to predict both the local heat 
transfer coefficient and the average value. The latter correlations are more 
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Fig. 6.17 Variation of Stanton numbers with Re gL 
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appropriate to compare the predictions of local, sectional, and average heat transfer 
coefficients with the present experimental data. As a result, four correlations of 
Soliman et al. (1968), Azer et al. (1971), Traviss et al. (1973), and Shah (1979), 
predicting both the local and the average heat transfer coefficients of condensation 
inside a tube are chosen. 

The backgrounds and techniques for developing those selected correlations are 
described in Chapter 1 . A brief introduction for each correlation follows. 

Correlation of Soliman. Schuster, and Berenson (1968) 

Their correlation is derived based on a theory of an annular condensation and 
experimental data of condensation inside both horizontal and vertical tubes. The 
condensing fluids in the experiments were water, ethanol, methanol, toluene, 
trichloroethylene, R-113, and R-22. The correlation is 

h x = 0.036— Pr, 06 y (F f +F m )p, (6-6) 

P. 

where 

h = 0.045Re~ L o 2 [X' 8 + 5.7(^i) 0 0523 (1 - X) 0 47 X ,M 

* M g 

x(^£) 0261 + 8.1 1(—) 0 105 (1 -X) 0M X a86 (^£) 0 - 522 ] 

Pi P g Pi 

— , = 0.5D H (— )[2(1 - X)(^£) 2/3 + (I - 3 + 2X)(— ) 4/3 

8rh t 2 /(ir 2 p g D H ) dx P. X P. 

+ (2X - 1 - /SX)(^*) l/3 + /}(— ) 5 ' 3 (2 - I - X) + 2(1 - X -0 + 0X)(^i)] 

P, P, x Pi 
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0 = 1.25 for turbulent film; /3 = 2.0 for laminar condensate film. Since in the 
present experiment, the flow of the condensate film is considered to be laminar, a 
value of 2.0 is chosen for /?. 

Correlation of Azer. Abis. and Swearingen (1971) 

The correlation is derived based on an analytical model of condensation and their 
own experimental data of condensation inside a horizontal tube of 12.7 mm I.D. and 
2.4 m long. In the experiment, R-12 was used as the condensing fluid; the range of 
the total vapor mass velocity varied from 137 to 448 kg/m 2 s; the quality (X) at the 
exit of the tube varied from 0.18 to 0.77. The correlation is 

h = 0. 153— Pr.Re„ 0 L 9 (— )(— ) 0J X 09 — (6-7) 

d h «, P, C 


where 

= 1 + 1.0986x2' 039 

C = 3.88Pr, 0663 (4.67 - X) Xtt = A 01 (^)°- 5 (i-^) 09 

M g P, X 


Correlation of Traviss. Rohsenow. and Baron (19731 

The correlation is based on an annular flow model using the von Karman 
universal velocity distribution for the condensate film. The correlation is 

h x = 0.15— Pr ' Rei [_L + 2.85xn° 476 ] ( 6 ' 8) 

D h f 2 l Xn 

where 

F 2 = 0.707Pr,Re, 05 , Re, ^ 50 

F 2 = 5Pr, + 51n[l + Pr,(0.09636Re l ° 535 - 1)], 50 < Re, < 1125 
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F 2 = 5Pr, + 51n(l +5Pr,) +2.51n(0.00313Re° 812 ), Re, > 1125 



t 


Re, ■ 0(1 - X)D 


Correlation of Shah (1979) 

Shah developed his correlation from a broad base of experimental data of 
condensation inside horizontal, vertical, and inclined tubes with diameters ranging 
from 7 to 40 mm. The condensing fluids involved in the experiments were water, 
R-ll, R-12, R-22, R-113, methanol, ethanol, benzene, toluene, and trichloroethylene. 
The total vapor mass velocity varied from 10.83 to 211 kg/m 2 s. His correlation is 


where h, is the heat transfer coefficient for a single phase fluid in fully developed 
flow in a circular tube. From the Dittus-Boelter correlation 


and p r represents the reduced pressure (actual pressure/critical pressure). Eq. 6-6 was 
further integrated to derive an expression for the average heat transfer coefficient over 
a certain length of a condensing surface. 



( 6 - 9 ) 


Pr 




r (1-X) 1 * 3.8.X 1 76 0.04X 276 .,x, 

, O + —tSt^TZ - — ^-=tz — )Jx, 


( 6 - 10 ) 
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where X, is the quality at the beginning of the condensing surface and X 2 is the 
quality at the end of the surface. 

Shah’s correlation appears to have the simplest form among the four selected 
correlations. However, his correlation predicts a zero value for the local heat transfer 
coefficient at the leading edge of the condensing surface where the quality (X) is 1, 
which is not acceptable. As a result, this correlation is used only to predict the 
average heat transfer coefficients, while the other three correlations are used to 
calculate both the local and the average heat transfer coefficients. 

In using the above the correlations, the thermodynamic and physical properties of 
the vapor are evaluated at the saturation temperature, the properties of the condensate 
are determined at the arithmetic average of the saturation temperature and the 
condensing surface temperature. The quality (X) is calculated as follows. 

The quality is defined as 

X = 1 - 5 (6-11) 

m, 


The expression for the condensate flow rate (m,) is found by substituting Eq. 4-16 into 
Eq. 4-24 


Wk, , * AT 


m, = 


-dx 


Using the film thickness correlation (Eq. 4-17) and assuming AT X to be constant, the 
above equation becomes 


X = 1 - 


W K AT } 
c,(l-c 2 ) m,h f 


l-c. 


( 6 - 12 ) 


Fig. 6.18 illustrates a typical distribution of the quality along the condensing 
duct. The maximum deviation between the calculated exit quality and the measured 
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Fig. 6.18 Typical variation of quality along the condensing surface 
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Fig. 6.19 Comparison of predictions for local heat transfer coefficients from correlations 

of Azer, Soliman. and Tranviss with experimental data 
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value (which is determined using the expression: X c = 1 - m lt / rh,) is less than 3.5 % 
for all the 51 sets of experimental data with the film thickness measurement. 

Fig. 6.19 shows a comparison of the predicted local heat transfer coefficient from 
correlations of Azer, Soliman, and Traviss with the present experimental data. The 
prediction from Traviss’s correlation agrees well with the experimental data. Both 
Azer and Soliman ’s correlations give a nearly constant local heat transfer coefficients 
along the condensing surface. All the three correlations predict lower local heat 
transfer coefficients than the present experimental data. From a comparison of the 
results with predictions from the three correlations, it was found that the predictions 
are not good when the exit quality is less than about 0.5. 

The four correlations of Soliman et al. (1968), Azer et al. (1971), Traviss et al. 
(1973), and Shah (1979) are used to calculated sectional heat transfer coefficients. 
Assuming constant AT, the following equation is employed to estimate the sectional 
heat transfer coefficients using the Azer, Soliman, and Traviss correlations 


H = 


L-L 


-Y h Ax 

/ ' XI 1 


1 > = I 


( 6 - 13 ) 


while Shah’s correlation (Eq. 6-10) for the average heat transfer coefficients can be 
used directly. 

The comparison between the predicted sectional heat transfer coefficients and the 
present experimental data are shown in Fig. 6.20. The experimental data are about 
within ±50 % of the predicted values. The average deviations are shown in Table 
6.5. It is seen that Soliman’ s correlation has a fairly low deviation in section 1 but 
was higher deviations (> ±30 %) in sections 2 and 3; Azer’s correlation leads to 
higher deviations (> ±30 %) in all the three sections; Traviss’s correlation gives the 
highest deviation in section 1 , but also produces the lowest deviations in sections 2 
and 3; the deviations in all the three sections from Shah’s correlation are nearly the 
same (around ±29 %). 


6.37 


6. Results and Discussion 


O 

O 

C\J 

e 


t 

^Q- 
, tn 

ISZ 


-ion |~ 
300 


□ 


^2 

*s3 


+50 % 





l 00 


20C 


iOO 




500 


CexpJ , w/m 2 °C 

a. Comparison with the correlation of Soliman et. al (1968) 


o 

o 

c\j 

E 

\ 

$ 



i oo 


400 


500 


h s CexpJ , w/m 2 °c 

b. Comparison with the correlation of Azer et. al (1971) 

Fig. 6.20 Predictions of sectional heat transfer coefficients 
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c. Comparison with the correlation of Traviss (1973) 



d. Comparison with the correlation of Shah (1979) 

Fig. 6.20 Predictions of sectional heat transfer coefficients 
from other researchers' correlations 
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Table 6.5 Average deviations in comparison of sectional heat transfer coefficients 


Correlations 

Deviations (%) 

Section 1 

Section 2 

Section 3 

Soliman 

22.2 

37.0 

31.9 

Azer 

30.3 

39.9 

37.8 

Traviss 

50.8 

27.2 

24.6 

Shah 

27.9 

30.0 

29.5 


The four correlations are also used to predict the total average heat transfer 
coefficients (TTj over the whole range of the present experimental conditions. Since 
there is no film thickness measurement available when the vapor velocity is high, the 
quality can not be evaluated using Eq. 6-12. Instead, it is estimated by the following 
equation assuming a linear relation with the length (x) of the condensing surface 

v rh 

X = 1 - 2L_!I (6-14) 

L m ( 

Reviewing Fig. 6.18, it is seen that the quality changes nearly linearly along the 
length of the condensing surface except close to the leading edge. The distribution of 
the quality along the condensing tube was also studied by Shah (1979) based on a 
large amount of the experimental data. He showed that in estimating the heat transfer 
coefficient the error due to the assumption of a linear variation of quality is quite 
small if the change of the quality less than 40 % . Since the change of quality in the 
present experiment through the whole test section is less than 40 % for most of test 
runs, the use of Eq. 6-14 to estimate the quality distribution is appropriate. 

Fig. 6.21 is a comparison of the predicted total average heat transfer coefficients 
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from the four correlations with the present experimental data. It is seen that the 
correlations of Azer and Traviss over-predict the heat transfer coefficients when the 
vapor velocity is high, while the correlations of Shah and Soliman generally give 
good predictions of the heat transfer coefficients. The over-prediction for the heat 
transfer coefficients from the correlations of Azer and Traviss were also observed by 
Luu (1980) after he compared the correlations with his experimental data. Also, 
those trends were found in the original works of Azer and Abis (1971) and Traviss et 
al (1973). 

Table 6.6 lists the deviations as computed from Eq. 6-0 (where let h tf = TT t(pred) ) 
between the predicted total average heat transfer coefficients from each of the four 
correlations and the present experimental data. Correlations of Shah and Soliman are 
the best of the four, with prediction that is within ±22 % of the experimental value. 
The correlation of Traviss has the largest deviation. However, its prediction at lower 
vapor velocity is much better than the predictions from the other correlations (Fig. 
6.21c). 


Table 6.6 Average deviations in comparison of the 
total average heat transfer coefficients 


Correlations 

Deviations (%) 

Soliman 

22.2 

Azer 

29.9 

Traviss 

36.7 

Shah 

20.7 
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c. Comparison with the correlation of Traviss at. al (1973) 
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6.8.2. Comparison with Experimental Data in the Literature 

A comparison of the present experimental data with those of Akers and Rosson 
(1960) is shown in Fig. 6.22. The experimental data of Akers and Rosson was 
obtained with a 0.31 m horizontal condensing tube with a 9.5 mm I.D. Methanol, R- 
12, and propane were used as the condensing fluids. It is seen in the figure that the 
present experimental data display nearly the same trend as those of Akers and Rosson. 
There is an apparent difference between the two group of data when the parameter 
[D H G//X|(p,/p g )' A ] is low. However, this difference reduces when D H G/p,(p 1 /p g )' /5 
becomes larger. It is probably because the length of their condensing tube is similar 
to that of section 1 (L, = 0.31 m) of the present condensing surface, the value of the 
present Nu! appears closer to those of their Nusselt numbers (Fig. 6.22a) than Nuj 
and Nu t . 

6.9. Predictions by Analysis 

The analytical model described in Chapter 5 predicts the film thickness, the 
interfacial shear stress, and the local and the average heat transfer coefficients. The 
predicted values are compared with the present experimental data. 

A key step in solving the analytical model is the determination of the interfacial 
shear stress. In the model, this shear stress is considered as the sum of an adiabatic 
shear stress and an equivalent shear stress due to the condensation of vapor on the 
vapor-condensate interface (See Eq. 5-16). The friction factor of the adiabatic shear 
stress is calculated in Eq. 5-18 

f c = c,Re"' 

where the most appropriate values of the constants (Cf and n f ) are found from the 
available correlations for the local adiabatic friction factor and a trial procedure. 
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Fig. 6.22 Comparison of the present experimental data with that of Akers and Rosson (1960) 
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For a single phase turbulent boundary layer over a smooth flat plate, Schlichting 
(1979) proposed the following correlation for the local adiabatic friction factor 

f = 0.0592Re;° 2 (6-15) 

In the model, it is assumed that a vapor boundary layer exists over a flat condensate 
film. Although the vapor boundary layer may not be a turbulent one all along the 
condensate film, considering the difference between single-phase flows and two-phase 
flows, the uncertainty of the inlet conditions, and the idea of using a trial method to 
model the shear force, Eq. 6-15 is first chosen to evaluate the local adiabatic friction 
factor in Eq. 5-18. 

Comparing Eq. 6-15 with Eq. 5-18, the constants, c f and n f , in Eq. 6-15 are 
found to be 0.0592 and -0.2 respectively. The model (Eq. 5-35) was solved with 
these values. The results from the initial calculation of the model are compared with 
the experimental values to evaluate the validity of Schlichting ’s friction equation in 
the model. 

The film thickness predicted from the analytical model is compared with that 
measured in the experiment. This comparison is for Re gL < 1,120,000 for R-113 and 
Re gL < 1,550,000 for FC-72 in the experiment since for higher values of Re gL , 
interfacial waves set in and the film thickness measurement is not available. Fig. 

6.23 shows some typical plots of the predicted film thickness and the experimental 
data for AT* 20 °C. The predicted film thickness values are always larger than the 
experimental values. When the inlet vapor Reynolds number (Re gL ) increases the 
difference between the predicted film thickness and the experimental value reduces. 

Predicted sectional heat transfer coefficients are compared with the measured 
sectional heat transfer coefficients for the whole range of the experiment. Fig. 6.24 
shows some typical comparisons. The experimental data used in this figure are the 
same as those used in Fig. 6.23. The deviation between the predicted value and the 


6.46 


mm 6 , mm 


6. Results and Discussion 


l 




a. R-113; AT = 21.4 °C; v„ = 0.5 m/s 
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Fig. 6.23 Comparisons of predicted film thickness with measured values (c f = 0.0592) 
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Fig. 6.23 Comparisons of predicted film thickness with measured values (c f = 0.0592) 
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experimental data for each section is shown in the figure. When the inlet vapor 
Reynolds number is less than 600,000, the agreement of the predicted value with the 
experimental data is within +10 % in the first section; however, in the remaining two 
sections, the difference between the predicted value and the experimental data is large 
(Fig. 6.24 a, b, e, and f); with further increase in the inlet vapor Reynolds number, 
the difference reduces substantially (Fig. 6.24 c, d, g, and h); when the vapor 
Reynolds number increases beyond about 2,000,000, the difference again becomes 
large (Fig. 6.25). 

Fig. 6.26 shows a comparison of the total average heat transfer coefficients 
predicted from the analytical model using Schlichting’s correlation for the adiabatic 
friction factor with those measured in the present experiment. The predicted total 
heat transfer coefficient agrees quite well with the experimental data when the inlet 
vapor Reynolds number (Re gL ) ranges from 500,000 to 1,500,000 for both R- 113 and 
FC-72. However, at very low and high Re gL , the model under-predicts the heat 
transfer coefficients in most cases. The predicted heat transfer coefficients for AT = 
10 °C are close to the experimental values. 

The prediction from the model with Schlichting’s friction factor agrees with the 
experimental data within +30 % when 300,000 < Re gL > 1,7500,000. Beyond this 
range of Re gL , the model under-estimates the heat transfer coefficients. The reason 
for the underestimation of the heat transfer coefficients at very low Re gL is not yet 
clear. 

At high Reynolds number, interfacial waves set in. The waves affect the heat 
transfer rate in two ways. It has been shown by Suryanarayana (1972) that the rate of 
the heat transfer rate with a mean condensate film thickness (5) with interfacial waves 
of amplitude (a) to the heat transfer rate without waves is [l-fa/TT) 2 ]' 1 . The waves are 
likely to increase the turbulent shear stress as they act as surface roughness on the 
vapor boundary layer, leading to a thinner condensate film thickness. Thus the 
increase in heat transfer rate due to the interfacial waves may be the result of both the 
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with measured values (c, = 0.0592) 
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Fig. 6.26 Predicted total average heat transfer coefficients (c f = 0.0592) 
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effects. As the amplitude of the waves is not available, only the effect of 
condensation on the shear stress can be accounted for. In Eq. 5-16, the interfacial 
shear stress is considered to be composed of a momentum shear stress (r^ due to 
condensation and an adiabatic shear stress (r,). A comparison of the values of the 
two components at higher inlet vapor Reynolds number is given in Fig. 6.27, where 
7 a is estimated using Schlichting’s expression for the friction factor. As seen from 
this figure, Near the leading edge of the condensing surface, the value of the 
equivalent shear stress due to momentum is much higher than that of the adiabatic 
shear stress; the difference between the two values decreases with increasing distance 
from the leading edge; beyond a certain distance, the adiabatic shear stress tends to be 
slightly higher than that of the equivalent shear stress. The figure shows that the 
contribution of both the adiabatic shear stress and the momentum shear stress should 
be considered. To determine an expression for r, that will better predict the 
condensation heat transfer, the shear stress (Eq. 5-16) needs to be modified. 

Based on some available experimental data, Jensen and Yuen (1982) derived a 
correlation for adiabatic shear stress for two phase flows in a horizontal channel. For 
a smooth interface, the correlation is expressed as 

T a = 5 X I0' 3 p g (v g - V,) 2 (6- 16a) 

A different form of the correlation is given for a wavy interface with a developing 
flow [The definition for the fully developing flow is the distance ^ 37/(H-S)] 

t, = 6.02xl0- 7 p g [A]'*(v g - Vl ) 3 * < 6 - 16b ) 

Wavy interface is assumed when the dimensionless parameter (v g -v,)/(4<rg/pi) I/4 , is 
greater than 17. 

Using Jensen and Yuen’s correlation for the adiabatic shear stress, the analytical 
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Fig. 6.27 Comparisons of adiabatic shear stress (r.) with momentum shear stress (r m ) 

(c, = 0.0592) 
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Fig. 6.28 Predicted total average heat transfer coefficients using Jensen 
and Yuen's friction equation 
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a. R-113; AT = 21.7 °C; v„ = 1.3 m/s 


b. R-113; AT = 20.2 °C; v, = 1.8 m/s 
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c. FC-72; AT = 20.1 °C; 1.0 m/s 


d. FC-72; AT = 20.0 °C; v„ = 1.3 m/s 


Fig. 6.29 Predicted film thickness using Jensen and Yuen's friction equation 
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Fig. 6.30 Predicted sectional heat transfer coefficients using Jensen and Yuen's friction equation 
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model was re-computed. Fig. 6.28 shows a comparison of the predicted total average 
Nusselt numbers from the model using Jensen and Yuen’s friction correlation with the 
experimental data. The average deviation between the re-computed value and the 
measured value is 44.0 % for R-113 and 45.3 % for FC-72 in comparison with 27.8 
% for R-113 and 31.3 for FC-72 (See Fig. 6.26) using Schlichting’s friction equation. 
Comparing the newly computed film thickness (Fig. 6.29) and sectional heat transfer 
coefficients (Fig. 6.30) with the experimental data, it is further found the predictions 
of the model using Schlichting’s friction equation are better than those using Jensen 
and Yuen’s correlation. 

Schlichting’s friction equation applies to turbulent single phase flow over a 
smooth flat plate. But, the vapor-condensate interface may not be always smooth. 
When interfacial waves set in, the surface of the condensate become very rough, 
which increases the adiabatic shear stress. Once this occurs, the use of Schlichting’s 
friction correlation may result in an underestimate of the interfacial shear stress and 
the heat transfer coefficients. One possible way to account for the effect of the 
interfacial waves on the shear stress is to modify the Schlichting friction correlation. 

Schlichting friction correlation 


was modified by adjusting only the 
constant (c f ) in Eq. 5-18. The 
model was computed repeatedly, 
each time with a new value of c f 
using 100 typical experimental runs. 
The variation of average deviation 
(S) between the predicted heat 
transfer coefficients and the 
experimental values with c f is 



Cf 

Fig. 6.31 Effect c f on S 


shown in Fig. 6.31. The average deviation decreases when c f is increased from its 
original value (c f = 0.0592); the deviation reaches a minimum value at c f = 0.245; 
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however, when c f is further increased, the deviation increases again. Based on this 
comparison, c f was modified as 0.245, while n f remains at its original value of -0.2. 

Substituting c f = 0.245 and n f = -0.2 into Eq. 5-18 for the friction factor, the 
analytical model was re-computed. Fig. 6.32 shows a comparison of the total average 
Nusselt numbers predicted from the model with those measured in the present 
experiment. Examining Fig. 6.32 and Fig. 6.26, it is seen that the average deviation 
between the predicted Nu, and the experimental data is reduced from 27.8 % to 16.3 
% for R-113 and from 31.3 % to 21.2 % for FC-72. Also, the newly predicted film 
thickness is closer to the experimental value (Fig. 6.33). Fig. 6.34 shows 
comparisons between the predicted sectional heat transfer coefficients with the 
experimental data. It is seen in this figure that the predicted sectional heat transfer 
coefficients approach the experimental data in section 2 and section 3; however, in 
section 1, the model with c f = 0.245 predicts a lager sectional heat transfer 
coefficients than the experimental values. At higher Re gL , the model still under- 
predicts the sectional heat transfer coefficients (Fig. 6.35), but the difference between 
the predicted value from the model using c f = 0.245 and the experimental data is 
reduced in comparison with the prediction with c f = 0.0592. 

With c f = 0.245, the performance of the model improves over the whole 
experimental range. For 500,000 ^ Re gL $ 1,500,000, the prediction of the model 
agrees well with the experimental data. For this range of Re gL , the average deviation 
between the predicted total average Nusselt number and the experimental value is less 
than ±15 % for R-113 and less than ±17 % for FC-72. 

Narain and Kamath (1991) developed a model to simulate the process of 
condensation of a pure saturated vapor flow between two horizontal parallel plates. 
They identified the friction factor equation for the interface (with the help of the data 
collected in the present experiment) for both laminar and turbulent flows of vapor. 
Narain and Kamath further used their model to predict the film thickness of 
condensate in the present experiment. A comparison between the film thickness 
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Fig. 6.32 Predicted total average heat transfer coefficients (c, = 0.245) 
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a. R-113: AT = 21.7°C;v ai = 1.3 m/s b. R-113;AT = 20.2 °C; v,, = 1.8 m/s 
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Fig. 6.33 Predicted film thickness (c, = 0.245) 
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Fig. 6.34 Predicted sectional heat transfer coefficients (c, = 0.245) 
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Fig. 6.35 Predicted sectional heat transfer coefficients at higher Re gL (c f = 0.245) 
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Fig. 6.36 Comparison between the present model and Narain's model 
R-113; AT = 39.7 °C; v* = 1.28 m/s 

prediction of Narain’s model and that of the present model (using Cf = 0.245) is 
given 


in Fig. 6-36. It is seen that the predictions from both the models are in good 
agreement with the experimental data. However, it should be noted that the 
experimental measurements of heat transfer were used to derive the friction factor. 

Fig. 6.37 presents typical distributions of the predicted vapor velocity and the 
liquid velocity along the condensing surface. The vapor velocity and the liquid 
velocity are calculated using Eq. 5-6 and Eq. 5-22. It is seen from this figure that the 
vapor velocity decreases along the condensing surface from 1.28 m/s to 1.16 m/s, 
while the liquid velocity increases slightly from 0.05 m/s to 0.07 m/s. 

Fig. 6.38 shows a typical distribution of the predicted interfacial shear stress 
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Fig. 6.38 Distribution of the interfacial shear stress along the test section 
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along the test section. The interfacial shear stress decreases shaiply near the leading 
edge of the condensing surface, and then reduces slowly along the test section. 

6.10. Heat Transfer Correlations 

Two different correlations for calculating the average heat transfer coefficient in 
condensation inside a horizontal rectangular duct are derived. One is a semi- 
empirical correlation based on the condensation model and the present experimental 
data. The other is an empirical correlation based merely on the present experimental 
results and dimensionless parameters. 

6.10.1. Semi-empirical Heat Transfer Correlation 

Based on the analysis on the condensation model, the average Nusselt number 
(NUi) is found to be possibly a function of Re gL , Pr/Ja, ujn 1 > p/p g > and L/L (Section 
6.7.1). Hence, it is assumed that 

N U| = a 0 (^l) 4 -[Re (6-17) 

Ja gL L /*, Pg 

where ao, a,, a 2 , a 3 , and a 4 are constants, which can be determined from the 
experimental results. 

Since 

Re ,h = Re , (6-18) 

gL L gL, 

Eq. 6-17 is rearranged as 



Based on a regression analysis on Eq. 6-19 using the present experimental data, a 
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correlation for the average Nusselt number is derived 

NUj = 0.00292Re g ° L 7ll (— ) 0 0853 (— ) 2 537 (— )' 1103 ( 6 - 20 ) 

' Ja ^ P t 

The range of the experimental data used in the regression analysis for Eq. 6.20 is 
listed in Table 6.7. 


Table 6.7 Ranges of experimental data used in correlating Eq. 6-20 


Fluids 

Re gL x 10- 3 

Pr,/Ja 

/Vm g 

P\'P * 

R- 113 

210-2,700 

24.0 - 104.7 

44.6 - 48.46 

177.7 - 212.6 

FC-72 

380 - 3,322 

13.7 - 79.4 

36.1 - 39.8 

112.7- 129.7 


Fig. 6.39a shows a comparison of the calculated average Nusselt numbers with 
the experimental values. The calculated Nusselt numbers agree fairly well with most 
of the experimental data. The average deviation between the calculated Nusselt 
numbers and the measured values is ±17.5 % for the entire condensing surface. 
However, when the inlet vapor Reynolds number is very small (approximately when 
Re gL < 445,000), the deviation is high. At the lowest Reynolds number (Re gL = 
202,338), the difference is the highest (~ ±50 %). 

Eq. 6-20 can be used to predict the average Nusselt number over any specific 
length of the condensing surface from its leading edge by replacing Lj in Eq. 6-20 
with that specific length. As a result, the average Nusselt number (Nu 2 ) over sections 
1 and 2 m = L 2 ) and the average Nusselt number (Nu,) over section 1 (L, = L) are 
calculated. Comparisons of Nu 2 and Nu, with the experimental data are presented in 
Fig. 6.39b and Fig. 6.39c. The average deviation between the predicted Nusselt 
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Fig. 6.39 Comparisons of the heat transfer correlation (Eq. 6-20) with experiment 
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numbers and the experimental values increases with decreasing length of the 
condensing surface. Based on the whole length of the condensing surface, the 
average deviation is only ±17.6 %, while for the section 1, the deviation increases to 
±24 %. The reason for increasing the deviation with decreasing the length of the 
condensing surface is not clear. Eq. 6-26 is probably more relevant for predicting the 
average Nusselt numbers with an appropriate longer length of condensing surface. 

6.10.2. Empirical Heat Transfer Correlation 

It has been shown previously in Fig. 6-17 that the Stanton number with an 
average heat transfer coefficient can be related to the inlet vapor Reynolds number, 
vapor. Therefore, it may be assumed that 

St, = foRe,L, (6 - 2l > 

where the values of the constants, f 0 and f,, are determined using the present 
experimental data involving only the average heat transfer coefficient (TTd and the 
regression analysis. The correlation so obtained is 

St, = 0.00573Re gL °- 274 (6-22) 

The range of Re gL used for the regression analysis of Eq. 6-22 is listed in Table 6.7. 

The Stanton number with the average heat transfer coefficient over the entire 
condensing surface (L ( = L) is calculated from Eq. 6-22. Also, the Stanton number 
with the average heat transfer coefficient over section 1 and section 2 and the Stanton 
number with the average heat transfer coefficient over section I are calculated by 
replacing L ( with Lj and L, respectively. Predictions from Eq. 6-22 are compared 
with the experimental data (Fig. 6.40). The average deviation (Eq. 6-0) between the 
predicted and experimental values is ±20 %. However, as with Eq. 6-20, when Re gL 
is approximately less than 445,000, the difference between the calculated Stanton 
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Fig. 6.40 Comparison of the heat transfer correlation (Eq. 6-22) with experiment 
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number and the measured values is high. The maximum deviation at the lowest Re gL 
(Re gL = 202,338) is about 67.18 %. 

6.10.3. Comparison between Heat Transfer Correlations 

In Section 6.8, the heat transfer correlations of Soliman, Azer, Traviss, and Shah ’* 
are discussed. Among the four correlations, it is found that the correlations of Shah 
(Eq. 6-10) and Soliman (Eq. 6-6) are better than the other two in predicting the total 
average heat transfer coefficients of the present experiment. The correlations of Shah 
and Soliman are compared with Eq.s 6-20 and 6-22 on the basis of the present 
experiment. It appears that all the predictions for Ti, from all the four correlations are 
in good agreement with the experimental data except for low values of Re gL (Fig. 

6.41). At low values of Re gL , Shah’s correlation greatly underestimates the heat 
transfer coefficients, while Soliman’s correlation overestimates the heat transfer 
coefficients. Eqs. 6-20 and 6-22 also underestimate the heat transfer coefficients 
when Re gL is low, but the difference between the experimental data and the predicted 
values from Eqs. 6-20 and 6-22 is small than that from the correlations of Shah and 
Soliman. 

6.11. Interfacial Waves 

6.11.1 Observation of the Interfacial Wave Initiation 

The interfacial wave initiation was observed in two separated experiments: (1) 
condensation experiment; (2) adiabatic air-liquid flow experiment. 

Interfacial Wave initiation in Condensation • 

In the condensation experiment, the initiation of the interfacial waves on the 
condensate film (Fig. 6.42) depends on the inlet vapor Reynolds number (Re gL ) and - 

AT. As AT increases, the value of Re gL initiating the waves decreases (Fig. 6.43). 

This may be because higher AT leads to a larger condensate film thickness (5) and a 
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Fig. 6.41 Comparison of different heat transfer correlations 
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Fig. 6.42 Photograph of interfacial wave initiation in condensation 
FC-72; AT = 9.7 °C; v„ = 1.1 m/s; Re gL = 1,225,835 
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Fig. 6.43 Relationship between Re gL causing the interfacial wave initiation and AT 
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a. Outlet view 



b. Side view 


Fig. 6.44 Photographs of interfacial wave initiation in air-liquid flow 

a. p, = 5.2 CP; v g = 3.2 m/s; rh, = 1.8 x 10' 3 kg/s 

b. p, = 4.5 CP; v 9 = 6.0 m/s; m, = 6.1 x 10"* kg/s 
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larger 8 is more likely to cause the instability of the condensate film than a smaller 8. 
Interfacial Wave Initiation in Adiabatic Air-liquid Flow 

Typical photographs of the initiation of the interfacial waves in the experiment of 
the air-liquid flow are shown in Fig. 6.44. One of the photographs in the figure (Fig. 
6.44 a) was taken at the outlet opening of the test section and the other (Fig. 6.44 b) 
was taken from the side of the test section. From a comparison between Fig. 6.44 
and Fig. 6.42, it is seen that the structures of the interfacial waves in the 
condensation and those in the air-liquid flow are similar (both are two dimensional 
waves); the waves appear somewhere on the liquid film and then propagate towards 
the outlet of the test section. As illustrated in Fig. 6.45, the distance between the 
inlet of the test section and the first wave is defined as the onset length (LoJ. 
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Tesct Section 
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First Wave 

Lon 




Wavy Flow 

Liquid 

Smooth Row 


Fig. 6.45 Illustration of onset length 


For both the condensation and the air-liquid experiments, the relationship between 
the onset-length and the liquid flow rate is show in Fig. 6.46, where m, in figure a is 
the condensate rate at the point where the first wave (Fig. 6.45) appears and m, in 
figure b is the liquid flow rate. In the condensation experiment, the onset length does 
not change much with increasing m, (L^ = 0.4 m). In the air-liquid experiment, the 
onset-length varies from 0 to 0.75 m and the relationship between the onset length 
and the liquid flow rate displays two different trends: (1) increases when jx, = 

0.91 CP; 2) L ot decreases when /x, = 4.46 CP (figure b). Those conflicting trends of 
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Fig. 6.46 Variation of onset length with liquid flow rate 
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L ot with m, may indicate that L OT is influenced by other parameters besides the liquid 
flow rate. 

Fig. 6.47 shows some typical variations of the liquid film thickness in the air- 
liquid flow experiment. It is seen that the liquid film thickness decreases along the 
test section. This trend is different from that of the condensate film thickness in the 
condensation experiment (Fig. 6.4). This difference may be attributed to the different 
driving forces in the two separate experiments. In the condensation experiment, the 
condensate is mainly driven by the momentum of the condensate and the shear stress 
of the vapor flow; in the air-liquid flow experiment, the flow of the liquid is caused 
by the hydraulic gradient of the liquid and the shear stress at the interface. 
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Fig. 6.47 Variation of liquid film thickness along the test section, air-liquid flow 


Because the distributions of the liquid film thickness in condensation and the air- 
liquid flow experiments are different, the variations of the liquid velocity along the 
test section may differ. However, those differences do not affect the following 
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discussion since only the liquid velocity at the point when the first wave appears is 
considered. 

Fig. 6.48 shows the relationship between the air velocity (v g ), at which wave 
appears, and the liquid velocity (v,) calculated using Eq. 4-23. The air velocity at 
which waves appear decreases with increasing liquid velocity. Fig. 6.49 shows the 
effect of the liquid viscosity (/*,) on the air velocity for wave initiation. It is found 
that this air velocity increases with increasing liquid viscosity; when the liquid 
viscosity and the liquid flow rate are greater than 2.37 CP and 0.002 kg/s 
respectively, the increase of the liquid viscosity is insignificant for the wave initiation. 
However, when the liquid flow rate is less than 0.002 kg/s, the effect of the liquid 
viscosity on the initiation of interfacial waves is significant (Fig. 6.50) and the air 
velocity for wave initiation increases with increasing liquid viscosity. The liquid 
viscosity has a stabilizing effect on the initiation of the interfacial waves. 

The trends of the vapor velocity for wave initiation in condensate film seem to be 
similar to those of the air velocity for wave initiation. It is seen from Fig. 6.51 that 
the vapor velocity (at which waves appear) decreases with increasing condensate rate. 
The vapor velocity of R-133 for initiating the waves is slightly higher than that of FC- 
72. This may be because the viscosity of liquid R-l 13 is greater than that of FC-72 
(about 1.2 times). 

The experimental results indicate that the initiation of the interfacial waves is 
influenced by the gas (vapor or air) velocity, the liquid velocity, and the liquid 
viscosity. The effect of liquid viscosity decreases gradually with increasing liquid 
velocity. 

Fig. 6-52 shows a relationship between Re g and Re, at which waves appear, 
where Re g and Re, are defined as v t Uv g and v gj 5///,. For condensation v g and m, are 
the vapor velocity and the condensate flow rate at the point where the first wave (Fig. 
6.44) appears; for air-liquid flow, v g and m, are the air velocity and the liquid flow 
rate. At similar values of Re,, the vapor Reynolds number is apparently larger than 
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Fig. 6.48 Variation of v g for wave initiation with v, in air-liquid flow 


co 


E 


cn 

> 


10 
8 
6 
4 
2 
0 

0.000 0.004 0 008 0.012 0.016 



— 1 1 — 

1 

— i r 

j 

1 





• 

2.37 CP 

_ 




O 

3.26 CP - 





■ 

3.95 CP 

- 




□ 

4.46 CP 





A 

5.24 CP 

V* 




▼ 

7.55 CP - 

- 1 

1 ▼ 1 




- 

- \ 

Pa M 1 

^ * 



□ 

1 

\ 


A 

ft 



• 





_ 

• • 

• 



_ 


• i 


_J t 

• 

1 — 



mi, kg/s 


Fig. 6.49 Variation of v g for wave initiation with ih, in air-liquid flow 
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Fig. 6.50 Effect of liquid viscosity on wave initiation 
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Fig. 6.51 Variation of v^ for wave initiation with m, in condensation 
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Fig. 6.52 Damping effect of condensation on wave initiation 

the air Reynolds number. This may indicate that the condensate flow is more stable 
than the air- liquid flow. This damping effect of condensation on the initiation of the 
interfacial waves was also reported by Barry and Corradini (1988) 


6.11.2 Criteria Predicting the Instability of Liquid Flow 

For gas and liquid flowing parallel and horizontally, the instability of the liquid 
flow is discussed by Chandrasekhar (1961) assuming that the domain of each of the 
fluids was semi-infinite. The condition for the onset of waves is given 


( v g " v i) 2 > 


«i«2. 


<rg(ai - a 2 ) 
Pi + Pg 


( 6 - 23 ) 


where 




P. 

P. + P g ’ 


a. 


Pg 

P. + Pg 


Considering gas and liquid flowing horizontally between two horizontal plates, 
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r 


Considering gas and liquid flowing horizontally between two horizontal plates, 
Milne-Thomson (1960) derived the instability condition of the liquid 


(v g - v,) 2 > 


P, cothkH, + p g cothkH g ^g(p,-p g ) ^ ^ 
Pfi g cothkH, cothkH g k 


(6-24) 


where 

H, - depth of the liquid; 

Hg - depth of the gas; 
k - wave number (2x/X); 

X - wave length. 

For large values of the depths of the fluids making both kH,, and kH 2 much greater 
than 1 , Eq. 6-24 becomes 


( v g-v ,) 2 > 


Pg+P^g (PrP,) 
PgP, k 


+ <rk] 


(6-25) 


To find a value of k to make the right hand side minimum 


_d_ r g (PrPg) 
dx ^ k 


+ <rk] = 0 


o 


Substitution of the above equation into Eq. 6-25 yields 

(v - v,) 2 = 2^-^[<rg(p, - p g )] l/2 (6-26) 

P,P 2 

It is seen that Eq. 6-26 is the same as Eq. 6-23. 

For gas and liquid flowing inside a horizontal channel, Jensen and Yuen (1982) 
proposed a criterion to determine the transition of liquid flow from a smooth surface 
to a wavy surface 
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( y g - v .) 

w j 


> 17 


(6-27) 


where 


W, = 

p\ 


The only criterion applied for predicting the stability of condensate flow appears 
to be the one given by Kocamustafaogullari (1985) 

Ka5*1c* 4 - [g*5” + 2p * (— ) 2 + i-5*‘(l + — )]k*‘ - J_[^] ^ 0 (6-28) 
^ Pr 1 £ * S • > Pr J 


5*’ Pr 


where 


Ka = (4^-r 
v- Ap% x 


CAT 

Ku = -_il — 


5 * = 5[ P|Ap !*] 1/3 

2 

Pi 


k * = k( 


f* 1 y/3 


pAps 


T- 


'5* a* = °y 


g - = iz p* = — — 

g* p . 


ox 

For condensate flowing horizontally, the component of the gravity force in the x 
direction (gj is zero. In this case, Eq. 6-28 can be simplified as 

—(—) s; 0 (6-29) 

6*’ Pr 


Eq. 6-29 predicts that for condensation inside a horizontal duct, the condensate 
flow is always stable. This prediction is contradicted with the present and other 
experimental results. Hence, the criterion of Kocamustafaogullari may not be 
appropriate for condensation on a horizontal surface. 
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Eqs. 6-23, 6-24, and 6-27 are used to predict the instability of the present 
experiments. Some typical comparisons of the predictions from those equations with 
the experimental results of the air- liquid flow are shown in Table 6.8, where v gl 2 
represents (v g - v,) 2 and the values in a column with equation numbers are the right 
hand side values of the equations. It is seen that the right hand side values of both 
Eq. 6-23 and 6-24 are larger than the values of (v g - v,) 2 in the same row, especially 
when the air velocity (v g ) is small. The left hand side values of Eq. 6-27 is close to 
17 (the criterion value in Eq. 6-23), but it is in a lower side. When Eqs. 6-23, 6-24, 


and 6-27 are used to predict the onset of the waves in condensation, it is found that 
none of these equations give satisfactory results (Table 6.9). 

The present experimental results show that the instability of the liquid flow may 
also be a function of the liquid viscosity (Figs. 6.50 and 6.51). However, this 
viscosity influence is not taken into account in any of the three equations (Eqs. 6-23, 
6-24, and 6-27). To include the viscosity effect in Eq. 6-27, it is assumed that 


where a, and bj are constants and are determined only from the experimental data of 
the air-liquid flow and the regression analysis. As a result, Eq. 6-30 becomes 


Using Eq. 6-31 to evaluate the instability in the condensation experiment, it is 
found that the predictions from the equation are in good agreement with the 
experimental data (Fig. 6.52). However, some predictions from Eq. 6-31 for the 
instability of the air-liquid flow fall inside the smooth flow region, which is in 
contradiction with the experimental observation. Hence, Eq. 6-31 still needs more 
improvement and experimental support. 



(6-30) 


C± ~ > 2.0 

W, /i g 


(6-31) 
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Table 6.8 Comparison of the predictions from Eq.s 6-23, 6-24, and 6-27 with the 
experimental results of sir-liquid flow 


Run 

Mi 

v * 

V| 

V 

Eq.6-23* 

Eq.6-24 b 

Eq.6-27 c 

26 

1.00 

1.78 

0.01 

3.13 

43.75 

50.22 

7.66 

69 

1.12 

3.08 

0.02 

9.36 

43.48 

50.00 

13.23 

27 

1.71 

1.94 

0.00 

3.76 

44.43 

50.39 

8.51 

31 

2.37 

1.78 

0.01 

3.13 

44.55 

49.80 

7.88 

29 

2.37 

2.73 

0.01 

7.40 

44.55 

49.81 

12.11 

40 

2.67 

1.44 

0.05 

1.93 

44.38 

49.70 

6.18 

36 

2.67 

5.94 

0.04 

34.81 

44.38 

49.67 

26.23 

41 

3.26 

2.38 

0.09 

5.24 

44.51 

49.62 

10.23 

43 

3.26 

3.92 

0.03 

15.13 

44.51 

49.62 

17.38 

46 

3.26 

5.99 

0.01 

35.76 

44.51 

49.60 

26.72 

47 

3.95 

2.26 

0.05 

4.88 

44.42 

49.41 

9.90 

50 

3.95 

3.85 

0.01 

14.75 

44.42 

49.41 

17.20 

51 

3.95 

5.97 

0.02 

35.40 

44.42 

49.39 

26.66 

52 

4.46 

2.91 

0.05 

8.18 

44.93 

49.73 

12.89 

55 

4.46 

3.48 

0.01 

12.04 

44.93 

49.73 

15.63 

56 

4.46 

6.00 

0.01 

35.88 

44.93 

49.72 

26.99 

62 

5.24 

2.41 

0.04 

5.62 

44.97 

49.75 

10.68 

60 

5.24 

3.21 

0.01 

10.24 

44.97 

49.75 

14.42 

64 

7.55 

3.53 

0.03 

12.25 

45.36 

49.94 

15.87 

67 

7.55 

5.90 

0.00 

34.81 

45.36 

49.92 

26.75 

59 

5.24 

5.98 

0.02 

35.52 

44.97 

49.73 

26.86 


a: right hand side values of Eq. 6-23 
b: right hand side values of Eq. 6-24 
c: left hand side values of Eq. 6-27 


> 
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Table 6.9 Comparison of the predictions from Eqs. 6-23, 6-24, and 6-27 with 

the experimental results of condensation 


Run 

AT 

v « 

v, 

V 

Eq . 6-23 

Eq . 6-24 

Eq . 6 - 

27 

273 

39.6 

1.24 

0.013 

1.51 

4.09 

4.10 

8.80 

274 

31.48 

1.25 

0.016 

1.52 

3.95 

3.96 

8.83 

278 

22.42 

1.38 

0.016 

1.86 

3.91 

3.91 

9.74 

279 

10.75 

1.56 

0.017 

2.38 

4.01 

4.02 

11.00 

345 

40.58 

0.81 

0.013 

0.64 

1.79 

1.92 

6.63 

346 

28.21 

0.87 

0.015 

0.74 

1.85 

1.98 

7.10 

388 

19.53 

1.07 

0.021 

1.10 

1.86 

2.00 

8.66 

389 

9.69 

1.11 

0.012 

1.21 

1.72 

1.83 

9.05 
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Fig. 6.53 Validity of the instability criterion (Eq. 6-31) 
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6.11.3 Interfacial Wave Length and Speed 

To study the interfacial waves in condensation, the wave length and the wave 
speed were measured. The experimental data of the wave length include the local 
wave length measured along the condensing surface and the average wave length, 
which is the arithmetic mean of all the values of the local wave length. Since the 
measured wave speed is an average value over the length of the wave-speed 
measuring device (Fig. 4.19) and also no obvious change of the wave speed was 
found when moving the device along the test section, the data of the wave speed are 
considered as the average wave speed over the whole condensing surface. 

Fig. 6.54 presents two typical distributions of the two dimensional local wave 
length (L**) along the condensing surface. For both R-l 13 and FC-72 used as the 
condensing fluids, the largest wave length occurs near the leading edge of the 
condensing surface; the wave length decreases along the surface; at some distance 
down stream of the condensing surface, the two dimensional waves become three- 
dimensional waves. 

For the average wave length (LJ, it is seen in Fig. 6.55 that the effect of AT on 
L* is not quite significant and the influence of the inlet vapor velocity (v gi ) on L* is 
not clear (for R-l 13, it shows that L* decreases with increasing v gi ; for FC-72 it 
shows that L* does not change much with v gi ). 

Fig. 6-56 shows the relationship between the average wave length and the liquid 
Reynolds number. The liquid Reynolds number varies approximately from 40 to 130 
for R-l 13 and from 85 to 400 for FC-72. Within those ranges, the average wave 
length is approximately constant. This trend is similar to the experimental results 
summarized by Hishburg and Florschuetz (1982). 

The variation of the wave speed with the inlet vapor velocity using R-l 13 and 
FC-72 as the condensing fluids is shown in Fig. 6.57. It is seen that with increasing 
vapor velocity, initially the wave speed is mainly constant; when the vapor velocity 
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Fig. 6.54 Distributions of local wave length along the condensing surface 
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a. R-113 
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Fig. 6.55 Variations of average wave length with inlet vapor velocity 
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Fig. 6.56 Variations of average wave length with Reynolds number of condensate 
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further increases, the wave speed increases; the temperature difference (AT) generally 
has no noticeable effect on the wave speed except for AT * 10 °C, for which, the 

* wave speed is lower. 

Fig. 6.58 shows comparisons of the wave length and wave speed of R-113 with 

♦ 

those of FC-72. It is seen that the wave length of R-113 is always larger than that of 
FC-72 at the same vapor velocity; the wave speeds of both fluids are about the same. 
From the above results, two main conclusions are drawn 

— the wave length reduces along the condensing surface and finally the two 
dimensional wave becomes the three dimensional wave; 

— the wave speed increases with increasing the vapor velocity. 

One unexpected feature that was observed was that the interfacial waves with 
condensation were not perpendicular to the side plates but made an angle of 
approximately 30 ° to 45 ° with the side plates, as shown in Fig. 6.2. A slight 
rotation of the duct about its axis had no influence on the orientation of the waves. 
However, with air-liquid adiabatic flows, the waves were perpendicular to the side 
plates even when the bottom surface was tilted by rotating the duct about its axis. 
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Fig. 6.58 Comparison of wave length and wave speed of R-1 13 with those of FC-72 


7. CONCLUSIONS AND RECOMMENDATIONS 


• 7.1. Conclusions 

Condensation heat transfer in an annular flow regime with and without interfacial 
waves was experimentally investigated. The study included measurements of heat 
transfer rate with condensation of vapor flowing inside a horizontal rectangular duct 
and experiments on the initiation of interfacial waves in condensation, and adiabatic 
air-liquid flow. An analytical model for the condensation was developed to predict 
condensate film thickness and heat transfer coefficients. The conclusions drawn from 
the study are as follows: 

— The condensate film thickness was very thin (<0.6 mm). The film thickness, 
growing rapidly at the leading edge of the condensing surface, increased 
gradually along the surface. With increasing the inlet vapor velocity, the 
condensate film thickness decreased. 

— The average heat transfer coefficient increased with increasing the inlet vapor 
velocity. After the appearance of interfacial waves, the increasing rate of the 
heat transfer coefficient with the vapor velocity enhanced greatly (it appeared that 
the effect of interfacial waves became significant when Re gL > 1 x 10 6 ). 

£ — The local heat transfer coefficient decreased with the axial distance of the 

condensing surface, with the largest change at the leading edge of the test section. 

c 

— The prediction of the average heat transfer coefficients from the present analytical 
model agreed within ±17 % with the experimental data for 5X10 3 ^Re gL ^ 

7.1 
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1.5 xlO 6 . However, outside this range, the predicted average heat transfer 

coefficients were lower than the measured values with an average deviation of 

approximate -40 %. * 

The interfacial shear stress, which consisted of the momentum shear stress and • 

the adiabatic shear stress, appeared to have a significant effect on the heat 
transfer coefficients. In the analytical model, a modified Schlichting’s friction 
equation (Eq. 6-15) was used in calculating the adiabatic shear stress. 

The heat transfer correlation based on the present analytical model and the 
experimental data (Eq. 6-20) and the correlation developed from the present 
experimental results (Eq. 6-22) were in good agreement with the experimental 
data. The average deviation between the calculated heat transfer coefficients and 
the experimental values was ±18 %. 

In the experiment, the condensate flow along the condensing surface experienced 
a smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy 
flow. The change of the flow patterns depended on the vapor velocity and the 
viscosity of liquid. 

In the condensation experiment, the local wave length decreased with the axial 
distance of the condensing surface and the average wave length decreased with 
increasing inlet vapor velocity, while the wave speed increased with increasing 
vapor velocity. 

The heat transfer measurements are reliable. The difference between the heat % 

transfer rate obtained from a heat balance in the cooling channel and that obtained 

from a heat balance in the vapor duct was usually within ±5 %. * 

The ultrasonic technique was effective for measuring the condensate film 
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thickness when the surface was smooth or had waves of small amplitude. 

7.2. Recommendations for Further Study 

Although a large amount of data and results were obtained in the present study, 

further work is still required for a clear understanding of the condensation process. 

Based on the present study, the following recommendations are proposed: 

— Only two different condensing fluids were used in the present experiment. For 
investigating the effects of the parameters (Pr„ Ja, p,/p g , p,/p g ) on the 
condensation heat transfer, more fluids are needed. Also, for studying the 
geometric effect, experiment should be conducted inside different test sections 
with various height, width, and length dimensions. 

— The present experiment showed that for R-l 13 the effect of AT on the heat 
transfer coefficients was not significant, but for FC-72 this effect was not quite 
clear. More experiments are needed to study the effect of changing AT on the 
heat transfer coefficients. 

— A better analytical model may be achieved by obtaining experimental values of 
interfacial shear stress and including a vapor velocity profile inside the model. 
Also, since the heat transfer is influenced by the flow patterns (smooth flow or 
wavy flow), it is necessary to derive the different models according to the 
different flow patterns. 

— More work is need to study interfacial waves, including the wave length and 
wave speed. To support the criterion (Eq. 6-31) predicting the initiation of 
interfacial waves, future work should provide more data on the instability in 
condensation, as well as in the air-liquid flow with a wider range of liquid 
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viscosity and surface tension and, possibly, other parameters. 

— Besides the ultrasonic technique, some other methods (such as, using capacitance 
probe, conductance, and optical glass) may be considered to measure the 
condensate film thickness. 

— The average heat transfer coefficients were found to increase rapidly with 
increasing inlet vapor velocity after the appearance of the waves. However, 
when the vapor velocity was more than approximately 4 m/s, it seemed that the 
rate of the increase of the heat transfer coefficients started to decrease (Fig. 6.6). 
To verify this result, more experiments with higher vapor velocity are needed. 

To do this, two major modifications on the present experimental set-up of 
condensation are needed: (1) to add more heating power to the boiler, (2) to 
increase the cooling capacity of the auxiliary condenser. 


* 
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APPENDIX A 


UNCERTAINTY ANALYSES 


A.1. Total Average Heat Transfer Coefficient 

The total average heat transfer coefficient (TT*) is expressed as 

E = 3 (A-l 

‘ A(T Ml - T s ) 

The value of the total average heat transfer coefficient (Fj depends on the 
measurements of the heat transfer rate (q) the total area of the condensing surface 
(A), and the difference between the saturation temperature of vapor (T ut ) and the 
condensing surface temperature (TJ. The uncertainties in measuring q, A, T„, - T, 
will propagate through the calculations to produce an uncertainty in determining F t . 
Based on the method given by Kline and McClinrock (1953), the uncertainty of F t is 
estimated as flows. 

The basic equation for the uncertainty of F, is derived according to Eq. A-l 



(A-2) 


where 


dE t l 

dq ' A(T ut -T.) 


(A-3) 


5 

dA 


A 2 (T Mt -T.) 


q 


(A-4) 
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and 


3H t _ q 
aT .« A(T stl -T s ) 2 


_ q 

dT . ~ A(T Mt -T,) 2 


(A-5) 


(A-6) 


Substituting Eq.s A-3, A-4, A-5, and A-6 into Eq. A-2 and using Eq. A-l yield 



( C TMt +C T «) ll/2 

(T„ t -T .) 2 


(A-7) 


where 6 q , e A , e T „ t , and 6 X , represent the uncertainties of q, A, T ut , and T, and are 
estimated as follows. 


(i) 


The equation for the heat transfer rate (q) is expressed as 

q = m w C pw AT w 


(A-8) 


Based on Eq. A-8, the equation for e q is derived 

e = U ^ f f )i+( ) 2 V /2 +a 

q ^ dC Cm ' AAT fi V * 


pw 


3AT 


(A-9) 


where ri\, is the mass flow rate of coolant; AT W is the temperature difference of 
coolant across the entire cooling channel; A q is the possible heat transfer between the 
test section and the environment. 

Since the value of the specific feat of coolant (C pw ) has been well established, it is 
assumed that 


€ 


C w 


« 0 


S 


A.2 


Appendix A: Uncertainty Analyses 


With this assumption, Eq. A-9 simplifies to 


e = [(.d-S-g. ) 2 +( dq ) 2 ] l/2 +A 
q 3m w 3AT„, AT " q 


1 w 


From Eq. A-8 


= C AT 

am.. pw w 


(A- 10) 


(A- 11) 


aq 

aAT w 


=c_. m„ 


pw 


(A- 12) 


Substitution of Eq.s A-ll and A-12 into Eq. A-10 yields 

The uncertainty of e Aw is mainly due to the accuracy of the rotameter measuring 
the coolant flow rate. It is assumed that 

= 0.02m w (A- 14) 


The uncertainty of the temperature difference of coolant (e AT J is caused by the two 
uncertainties: (1) due to the inaccuracy of the thermocouple wire and junction (e te ) and 
(2) due to the uncertainty of the thermocouple readout (e me ). It is estimated that 

e te = 0.2 °C and € me - 0.1 °C 


therefore 

e 4Tw = 0.2 + 0.1 = 0.3 °C (A- 15) 

Substitution of Eq.s (A-14) and (A-15) into Eq. (A-13) and using Eq. (A-8) yields 
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«, = 0.02q[l +(_^L) 3 ]^+A j) (A-16) 

w 

The heat transfer between the test section and the environment (AJ includes the 
heat transfer between the vapor duct of the test section and the environment (q,) and 
the heat transfer between the cooling channel and the environment (%), i.e. 

A q = q, + q 2 (A-17) 

Here, q, and qj are calculated using one-dimensional heat transfer analysis. The 
parameters and their values used in the calculation are listed below: 

A, area of the side wall of the vapor duct in the test section, 0.05 m 2 ; 

A 21 area of the side wall of the cooling duct, 0.032 m 2 ; 

A 22 area of the bottom wall of the cooling duct, 0.036 m 2 ; 

A 23 area of the end wall of the cooling duct, 5.4 X 10 4 m 2 ; 

lq, thermal conductivity of brass, 110 W/m-K; 

kg thermal conductivity of fiber glass, 0.036 W/m-K; 

kp thermal conductivity of polycarbonate, 0.19 W/m-K; 

k, thermal conductivity of rubber gasket, 0.36 W/m-K; 

k, thermal conductivity of styrofoam, 0.033 W/m-K; 

h* convective heat transfer coefficient of coolant, 600 W/m 2 -K; 

h a convective heat transfer coefficient of environment (air), 10 W/m 2 -K; 

T w temperature of coolant, 5 °C; 

T 0 surface temperature of the flange of converging section connecting to the the 
front end of the cooling channel, 40 °C; 

T a temperature of environment (air), 20 °C. 

Calculation of q . 
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Due to the additional heating applied to 


the top wall of the vapor duct, it is 
assumed that the top wall is adiabatic and 
the heat transfer (q,) occurs through the 
side wall of the duct. Fig. A. 1 illustrates 
the side wall and the attached insulation 
involving the heat transfer. This figure 
shows that the temperature of the inner 
surface of the side wall is assumed to be 
the saturation temperature of vapor, T Mt . 
Referring to Fig. A.l, q, is calculated as 



Fig. A.1 Illustration of heat transfer 
through the vapor duct 


q. 


= A 


T_ - T 


k k, h_ 

p J » 


= 0.05 


20 - T. 


(A- 18) 


, 11.9 37 

(: + - 


0.19 0.033 


)x,°-=.-L 


= 0.04(20 - TJ 


Calculation of q^ 

The heat transfer between the environment and the cooling channel (q^ consists 
of heat transfer through the sides of the channel (qj,), heat transfer through the 
bottom of the channel (q^), heat transfer through the front end side of the channel 
(q 2 3 ), and heat transfer through the rear end side of the channel (q^). Therefore 

q 2 = q2. + q« + q23 + q 2 4 (A_19) 

Referring Figs. A.2 through A.4, qj,, ^ q 23 , and q^ are estimated in the following 
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** 32 

Dlm«ntions In mm 


Fig. A. 2 Heat transfer through the sides and the 
bottom of the cooling channel 



Converging 

Secdon 


Dimensions in mm 




Condensing surface 


10 


<_ Cooing Channel 
T* h w 



h. 


4 

* 


Fig. A. 3 Heat transfer through the front 
side of the cooling channel 


Fig. A. 4 Heat transfer through the rear 
side of the cooling channel 
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P 




f 

* 


Estimate of 

^21 = \\ 


s p o 


oo w 


= 0.032. 


20-5 


_L + _L +(_£L +1L1 +_L) X 10- 3 

10 600 0.033 0.19 110 


= 0.37 W 


Estimate of 


Q 22 ~ 


T_ - T. 


* '+± H n + ni )x io* 


= 0.027- 


20-5 


h . h 

= 0.92 W 


± + J-*J°- + ii£)xio* 

10 600 v 0.036 0.19' 


^23 - ^23 


T - T 


= 5.4X10 4 - 


40-5 


( 


1 


. + o^i2 )xlo . 


0.39 110 


600 


= 4.33 W 


Estimate of q^ 


T - T 

q 24 = a, 3 : * 

24 1 1 7 9 10 

Tr -+ ir + ^'T~ + T -) x10 ’ 3 

h » K k„ k 


= 5.4 xlO' 4 - 


20-5 


— +— +(— + 1 - ° ) x 10' 3 
10 600 M10 0.036 


= 0.02 W 


Substitution of the values of q^, q 2 2 , q 23 , and q 2 4 into Eq. A-19 yields 

q 2 = 5.64 W 


(A-20) 


A.7 


Appendix A: Uncertainty Analyses 


Substitution of Eqs. A- 18 and A-20 into Eq. A- 17, the heat transfer from the 
environment to the test section is found to be 

A q = 6.44 - 0.04T jat W (A-21) 


Substitution of Eq. A-21 into Eq. A- 16 yields 


f, = 0.02q[lt(Jl) ! ] 5 -0.04T„ + 6.44 


(ii) e A 

The total area of the condensing surface (A) is calculated by 

A = WL 


(A-22) 


(A-23) 


where W and L are the width and the entire length of the condensing surface. 
From Eq. A-23 


= rt-^e ) 2 +f— f ) 2 l l/2 

K aw W L ' J 


(A-24) 


where 




(A-25) 


Substitution of Eq. A-25 into Eq. A-24 yields 

6 a = [(L6 w ) 2 +(W €l ) 2 ] ,/2 (A-26) 

Both £ W and e L are related to the accuracy of machining the condensing surface. 

By assuming e w = 0.001 W and e L = 0.001L and using Eq. A-23, Eq. A-26 becomes 
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e = 1.41 x 10' 3 A (A- 27 ) 

A 

» 

* 

(iii) *r»t 

* 

The uncertainty of the measured saturation temperature of vapor (T ul ) results 
from the accuracies of the pressure gage measuring the vapor pressure in the test 
section and the barometer measuring the atmospheric pressure. It is estimated that the 
uncertainty of e Tsit is 

€ T .. t = 0.5°C (A-28) 


(iv) f T . 

The uncertainty of the measured surface temperature (e Tf ) is considered to be a 
summation of the following factors: 


— uncertainty due to thermocouple wire and junction (e^ ~ 0.2 °C); 

— uncertainty due to the data acquisition system (e, c = 0.01 °C); 

— uncertainty in the location of the thermocouple junction inside the 
condensing surface (e em ) . 

Hence 


C T. = f ,c 


+ €„ + €. 


(A-29) 


The value of is estimated as follows 

It is shown in Fig. A. 5 that the junction of the thermocouple is 7 mm deep inside 
the condensing surface. The junction of thermocouple may not coincide with the 
condensing surface, so that the temperature obtained by the thermocouple (T’ s ) is not 
same as the surface temperature (TJ. The difference between T’ s and T, is the 
uncertainty (e^J, which is estimated using one-dimensional heat conduction analysis 
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where A is the distance between the condensing 
surface and the junction of the thermocouple; k, 
is the thermal conductivity of tin (See Fig. 

A. 5); q t " is the heat transfer rate from the 
condensing surface to the junction of the 
thermocouple. Substituting e = 7 mm and k, 

= 66.6 W/m k into Eq. A-30 yields 



Fig. A. 5 Illustration of the surface 
thermocouple 


e = 1.05 x10V' °C 

cm 


Substitution of the expression for , the values of € te (0.2 °C) and e, c (0.01 °C) into 
Eq. A -29 yields 

e T =0.21 + 1.05 x10V" °C ( A ' 31 > 


By substituting Eq.s A-22, A-27, A-28, and A-31 into Eq. A-7, an equation for 
estimating the uncertainty of Ti, is derived 




[ 0 . 02(1 + 


15 2 _)i/2. 


6.44-0.04T 0.5 2 +(0.51 +1 .05 x 10'V f 


tat ^2 


at: 


r + - 


AT 2 


+ 1.99x 10' 6 
(A-32) 


The uncertainty increases with decreasing temperature difference of coolant (AT W ), 
and decreasing temperature difference between temperatures of saturated vapor and 
the condensing surface (AT). 


•* 


« 

¥ 


Sample Calculations 
(1) Run 180 (R- 113) 
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h t = 188.91 W/m 2 k, q = 229.83 W, T Sil = 51 °C, AT = 33.83 °C, and 
AT W = 2.68 °C. Substituting these values into Eq. (A-33) and assuming 
q t " = q/A = 6384.17 W/m 2 yields e ¥l = 25.48 W/m 2 k. Hence 

= 18891 ± 25 ‘ 48 W / m 2 -°C 

The percentage uncertainty of the total average heat transfer coefficient in Run 180 is 
13.49 %. 

(2) Run 360 (FC-72) 

¥, = 585.27 W/m 2 k, q = 407.75 W, T„ t = 56.66 °C, AT = 19.37 °C, and 
AT W = 2.44 °C. Substituting these values into Eq. A-32 and assuming 
q," = q/A = 11020.2 W/m 2 yields e ¥l = 82.54 W/m 2 k. Hence 

lu = 585 - 27 ± 82.54 W/m 2 - °C 

The percentage uncertainty of the total average heat transfer coefficient in Run 320 is 
14-1 Q %, 


A. 2. Inlet Vapor Reynolds Number 


The inlet vapor Reynolds number (Re gL ) is defined as 


Re tL - 

", 


(A-33) 


The uncertainty equation associated with Eq. A-33 is expressed as 

dRe , , 5Re . 

v dv „ - dL dv„ » 


Expansion of the partial derivative terms in the above equation yields 
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e Rv = Re gL [(^) 2 + + (^)T 2 (A-34) 

gi g 

where e n , e L , and € Vti are the uncertainties of v v L, and v gl . 

The value of v t is considered to be well established (v % is taken from the 
thermodynamics property table in the ASHRAE Handbook of Fundamentals). It is 
assumed that ~ 0. Also, similar to the consideration in deriving of Eq. A-27, e L 
is assumed to 0.001L. Hence Eq. A-34 becomes 

- ReJIO' 6 * < A ’ 35 > 


The vapor velocity (v gi ) at the inlet to the test section is calculated by 


rh t 

Pg WH 


(A-36) 


The uncertainty equation for v gi is expressed as 



(A-37) 


where is the uncertainty in m g due to the inaccuracy of the vapor rotameter 
measuring iiy, e H and e w are the uncertainties of H and W due to the inaccuracy of 
machining and assembling the vapor duct of the test section; is the uncertainty in 
determining the value of p g . It is assumed that 

= 0.01m t , = 0.01W, e H = 0.01H (A-38) 

Substitution of Eq. A-38 into Eq. A-37 and considering = 0 (p g is considered to be 
well established) yields 


V 
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\ = 0.0173v fi 

Substitution Eq. A-39 into Eq. A-35 leads to 

f Re . = 0.0173Re gL 


(A-39) 


(A-40) 


The percentage uncertainty of Re c . is 1.73 %. 


A.3. Total Average Nusselt Number 


The total average Nusselt Number (NuJ is defined as 

EX 

Nu. = -L. 

* k, 


(A-41) 


Based on the above equation, the uncertainty equation of Nut, is derived 


dNu. , 3Nu, , dNu. , ... 

V - H-jfV * * <TfV" 


3L 


dk, 


Expansion of the partial derivative terms in the above equation using Eq. A-41 yields 

<NU = Nu t [(^) 2 + A 2 + A 2 ] 1 ' 2 (A-42) 

E, L K i 

Assuming e L = 0.001L (uncertainty of L) and € b = 0 (uncertainty of k,), Eq. (A-42) 
becomes 


* Nu - nu,[A 2 + 10-r 2 

E. 


(A-43) 


Sample Calculations 
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(1) Run 180 (R- 113) 

Nu t = 2428.78; ¥ t = 188.91 W/m 2 -°C; e Kl = 25.48 W/m 2 o C (See page A-ll). 
Substituting these values into Eq. A-43 yields e Nui = 327.64. Hence 

Nu t = 2428.78 ±327.64 


The percentage uncertainty of Nu. in run 180 is 13.49 %, 

(2) Run 360 (FC-72) 

Nu t = 9606.34; ¥, = 585.27 W/m 2o C; e F , = 85.54 W/m 2o C (See page A-ll). 
Substituting these values into Eq. A-43 yields e Nui = 1354.81. Hence 

Nu, = 9606.34 ± 1354.81 


The percentage uncertainty of Nu. in run 360 is 14.10 %, 


A.4. Total Average Stanton Number 


The total average Stanton Number (StJ is defined as 


E t 

St, = ' 




Nu. 

or St, = 


‘ Re iL Pr, 


(A-44) 


Based on the above equation, the uncertainty equation of St, is derived 

e -[( dSt ' e ) 2 + ( 3St< f ) 2 + (^ e ) 2 1> /2 

€st » 3Nu Nu ' ^ Re,t ^ aPr Pr ' 


aRe gL 


apr, 


Expanding the partial derivative terms in the above equation using Eq. A-44 




(A-45) 
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where e P „ is the uncertainty of Pr ( . Assuming e Pn = 0 and using Eq. A-40, Eq. A-45 
becomes 

c_ = St t [(-!^) 2 + 2.99 x 1 O' 4 ] 1/2 (A' 46 ) 

s Nu t 

Sample Calculations 

(1) Run 180 (R-113) 

St, = 2.53 x 10- 4 ; Nu, = 2428.78; e Nu = 327.64 (See page A. 14). Substituting 
these values into Eq. A-45 yields e Stt = 0.344 x 10 -4 . Hence 

St t = 2.54X10' 4 ± 0.344X10" 4 

The percentage uncertainty of St in run 180 is 13.60 %. 

(2) Run 360 (FC-72) 

St, = 1.32 x 10"*; Nu, = 9606.34; e Nu = 1354.81 (See page A. 14). Substituting 
these values into Eq. A-45 yields e Sn = 1.32 x 10" 4 . Hence 

St t = 1.34X10' 4 ± 0.188X10' 4 
The percentage uncertainty of St, in run 320 is 14.21 %, 
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TEST OF CONVERGING UNIT 


The purpose of the converging unit was to obtain a uniform velocity at the inlet 
of the test section. A deflector was installed inside the unit (Fig. 4.6). The deflector 
was selected among eight deflectors of different shapes. 

The converging unit was tested with air as the working fluid. The air passed 
through the converging unit assembled with a deflector and packed with thin glass 
tubes. The air velocity was measured at 40 equally spaced locations at the outlet of 
the unit using an air velocity meter. The average deviation of the local velocities was 
evaluated for each testing deflector. The deflector corresponding to the lowest 
average deviation (±9.1 %) was selected to be used in the experiment. 

A schematic of this selected deflector is shown in Fig. B. 1 . The profile of the 
measured air velocity at the outlet of the converging unit is given in Fig. B.2, where 
each grid joint indicates each measured local velocity. It is seen that the velocity 




Fig. B.1 Schematic of the selected deflector 
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profile is basically uniform except that at the two sides of the unit, the velocity is a 
little lower. The maximum difference between the local velocity and the average 
velocity is -23.0 %. The average difference is ±9.1 %. 



Fig. B.2 Air velocity profile at the outlet of the converging unit 


APPENDIX C 


CALIBRATION OF FLOW METERS 


Vapor rotameter, venturimeters, and coolant flow meters were calibrated. The 
following describes the calibration and the calibration results. 

Vapor Rotameter 

The vapor rotameter was used to 
measure the inlet vapor flow rate to 
the test section of condensation. It 
was specially calibrated for use with 
R-l 13 by the manufacturer. It was 
also calibrated in the laboratory. 

Fig. C.l illustrates the system 
for the calibration of the rotameter. 

This system included the rotameter, a 
calibrated laminar flow meter, and two regulating valves, with air as the calibrating 
fluid. The rotameter was calibrated using air flow. The results of the present 
calibration (V^) were compared with the calibration data provided by the manufacture 
(V ftc ) and are shown in Fig. C.2. The values presented in the figure are based on one 
of the experimental conditions with T„ t = 55 °C, P„ t =1.14 bar, and using R-l 13 
as the condensing fluid. The maximum difference between the two different 
calibrations was ±2.6 %. The results of the calibration is given in Table C.l. In the 
table, the scale indicates the reading on the gage glass of the rotameter; the symbol 
SCMM means cubic meter of air per minute at the standard condition (the air 
pressure and temperature at 1 atm and 20 °C respectively). Having the flow rate of 


Air 



Fig. C.l Rotameter calibration 
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air at the standard condition, the flow rate of vapor in experiment was determined by 


V 


vap 



(C-l) 


where 

V* — volumetric flow rate of air at the standard condition, m 3 /s 
V vap — volumetric flow rate of vapor in experiment, m 3 /s 
p air — density of air at the standard conditions, kg/m 3 
p vap — density of vapor in experiment, kg/m 3 



Fig. C.2 Comparison between the laboratory calibration results of the vapor rotameter and 

the manufacturer data 
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Table C.1 Calibration data of vapor rotameter 


Vapor Rotameter 


Scale 

SCMNFxlO 3 

Scale 

SCMNPxlO 3 

3.54 

3.36 

15.22 

64.60 

4.80 

6.00 

16.50 

75.09 

5.55 

8.79 

17.58 

85.44 

6.40 

11.51 

20.25 

114.50 

7.20 

14.65 

22.75 

142.31 

8.10 

18.22 

25.45 

175.57 

9.10 

22.36 

27.60 

209.19 

10.93 

32.89 

29.85 

245.55 

12.45 

43.23 

32.90 

297.76 

13.85 

53.68 

35.10 

344.01 


Air 

I 


Venturimeters 

Venturimeters were used to measure 
the inlet and the outlet vapor flow rates 
of the test section in the condensation 
experiment and the air flow rate in the 
air-liquid flow experiment. All the 
venturimeters were made of brass and 
were calibrated using the laminar flow 
meter and air in the same calibration 
system as used for the rotameter 
(Fig. C.3). During the calibrations, the pressure drop across the venturimeter was 
measured with a diaphragm type pressure transducer. The pressure drop was then 
transformed into a voltage output via a carrier demodulator. The pressure transducers 
for the venturimeters were calibrated by Christodoulou (1987). It was shown in his 



Fig. C.3 Venturimeter calibration 


C.3 



Appendix C: Calibration of Flow Meters 


results that the measured pressure drop was linear with the output voltage. 

The calibration of the venturimeters was repeated several times. It was found 
that the results were quite repeatable (Fig. C.4). The results of the calibration of the 
venturimeters are listed in Table C.2, where venturi# 1 and venturi#2 were those used 
for measuring the inlet and outlet vapor flow rate respectively, and venturi (air) was 
used for measuring the air flow rate in the air-liquid experiment. 



Fig. C.4 Repeatability of the calibration of venturimeter 


Coolant Flow Meters 

Four rotameters were used to measure the coolant flow rate. All the rotameters 
were calibrated using water and measuring the time taken for a certain amount of 
water to be collected. The results of the calibration are given in Table C.3. 
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Table C.2 Calibration data of venturimeters 


Venturi# 1 

Vol. SCMM 2 x 10 3 

Venturi#2 

Vol. SCMM 2 x 10 3 

Venturi(air) 

Vol. SCMM 2 x 10 3 

0.29 

28.06 

0.33 

14.78 

0.10 

1.26 

0.48 

46.65 

0.96 

42.55 

0.47 

4.67 

1.03 

97.95 

1.42 

64.49 

1.27 

12.12 

1.49 

146.44 

2.06 

95.87 

1.75 

15.68 

2.02 

195.66 

2.76 

125.66 

2.46 

21.93 

2.50 

244.29 

3.48 

161.87 

3.78 

34.21 

3.05 

300.87 

4.20 

194.14 

5.87 

49.54 

3.48 

344.96 

4.88 

226.90 

7.27 

60.91 

4.04 

402.28 

5.18 

243.46 

9.10 

76.51 

5.06 

497.93 

5.70 

264.78 

— 

— 

6.34 

634.32 

6.05 

285.23 

— 

— 
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Table C.3 Calibration data of coolant rotameters (unit in kg/s) * 


Scale 

No. 1 

No. 2 

No. 3 

No. 4 

10 

0.0013 

0.0013 

0.0011 

— 

20 

0.0034 

0.0034 

0.0032 

0.0006 

30 

0.0059 

0.0062 

0.0061 

0.0015 

40 

0.0080 

0.0082 

0.0074 

0.0020 

50 

0.0102 

0.0109 

0.0100 

0.0026 

60 

0.0122 

0.0123 

0.0115 

0.0033 

70 

0.0143 

0.0145 

0.0140 

0.0039 

80 

0.0169 

0.0168 

0.0163 

0.0046 

90 

0.0188 

0.0190 

0.0185 

0.0052 

100 

0.0212 

0.0212 

0.0209 

0.0058 

110 

0.0234 

0.0237 

0.0228 

0.0063 

120 

0.0258 

0.0261 

0.0251 

0.0070 

130 

0.0279 

0.0282 

0.0274 

0.0076 

140 

0.0304 

0.0302 

0.0301 

0.0083 

150 

0.0318 

0.0321 

0.0316 

0.0088 
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CALIBRATION OF THERMOCOUPLES 


The thermocouples for measuring the vapor (condensate) temperature and the 
differential thermocouples for measuring the coolant temperature change in the 
cooling sections were calibrated in the laboratory. 

Vapor and Condensate Thermocouple 

Fig. D.l shows the calibration 
system. The thermocouples were 
connected to the data acquisition system 
(also used in the condensation 
experiment); the measuring junction of 
the thermocouple was placed inside an 
insulated bath along with a high 
accuracy thermometer. The 
thermocouple was calibrated at three 
different temperatures: 1) ice-water; 2) room; 3) boiling water. The temperatures 
measured by the thermometer and the thermocouple (through the data acquisition 
system) were compared with each other. Table D. 1 lists two typical calibration 
results. 

Differential Thermocouples 

The calibration system for the differential thermocouples is shown in Fig. D.2. 
The junctions of one differential thermocouple were placed inside two baths with 
different temperatures. The voltage generated by the differential thermocouple (Vol^ 


Thermometer 

/ 

Data Aqulstion 


System 



Thermercouple 

^Bath 


Fig. D.l Thermocouple calibration 
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Table D.1 Sample results of the thermocouple calibrations 


Condensate (°C) 

Boiler (°C) 

Thermometer 

Thermocouple 

Thermometer 

Thermocouple 

0.19 

0.20 

0.12 

0.13 

22.62 

22.48 

23.63 

23.89 

94.87 

94.91 

93.67 

93.60 


was measured by a micro-voltmeter. 

The difference (AT,) of the temperatures 
in the two baths was measured by two 
thermometers was compared with the 
temperature difference (AT 2 ) determined 
by 


AT 2 


Vol* 

0.039 



„ . . .. ..... Fig. D.2 Differential thermocouple calibration 

where 0.039 is the gradient of the linear 

curve of temperature-voltage for type T thermocouple with temperature ranging from 
0 - 20 °C (See the general thermocouple table). The results of the comparison 
between AT, and AT 2 for three differential thermocouples are given in Table D.2, 
where S denotes the average deviation and the term — chamber means the coolant 
mixing chamber (See Fig. 4.17). As seen in Table D.2, the maximum average 
deviation is less than ±4 %. 


D.2 


t » 
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Table D.2 Results of calibration of differential thermocouples 


Chamber 1-2 

(°C) 

Chamber 2-3 (°C) 

Chamber 3-4 (°C) 

AT, 

at 2 

AT, 

at 2 

AT, 

at 2 

0.83 

0.82 

0.35 

0.37 

0.90 

0.92 

1.17 

1.15 

1.10 

1.16 

1.29 

1.24 

2.00 

2.17 

1.69 

1.72 

1.42 

1.35 

3.06 

3.14 

2.60 

2.65 

2.23 

2.27 

4.04 

4.08 

2.71 

2.75 

2.74 

2.78 

4.95 

4.99 

3.08 

3.07 

3.67 

3.79 

6.01 

6.30 

3.43 

3.54 

6.70 

6.77 

6.37 

6.25 

3.63 

3.65 

8.75 

8.95 

8.34 

8.40 

4.19 

4.30 

9.42 

9.70 

21.30 

22.50 

4.84 

5.03 

21.19 

21.80 

S: 

3.8 % 

S: 3.2 % 

S: 

2.9% 
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CALIBRATION OF THE ULTRASONIC TRANSDUCER 

The ultrasonic transducer was calibrated prior to its use in the experiment. The 
calibration device consisted of a pointer attached to a micrometer traversing 
arrangement mounted on a small 
container (Fig. E.l). The transducer 
was placed underneath an aluminum 
plug inserted into the bottom copper 
plate of the container. The depth of the 
fluid in the container was determined by 
moving the pointer connected to the 
micrometer to the surface of the liquid. 

The time for the ultrasonic signal to 

pass through the liquid film was 

... ... Fig. E.l Schomatic of calibration device of 

measured with an oscilloscope. ultrasonic transducer 

Knowing the time and the thickness of 

the liquid film, the sonic velocity in the fluid was determined. The arithmetic average 
(which was within ±0.8 % of the measured values for R-113 and within ±1.3 % for 
FC-72) was taken as the sonic velocity at the measured temperature. The sonic 
velocity is function of the density of the fluid, which is temperature dependent. Reid 
and Sherwood (1958) suggested that the sonic velocity in a liquid be expressed as 

a = cj>] <E-1> 

where the constant c m was obtained from the measured velocities. Table E. 1 shows 
the sonic velocities computed by Eq. E-l and the comparison with the values reported 
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by Meyer (1969) for R-113 and by the 3M company (1989) for FC-72. The values of 
c m in Eq. E-l determined from the calibration is 1.97 x 10‘ 7 m‘°/kg 3 s for R-113 and 
1.08 x 10' 7 m 10 /kg 3 s for FC-72. 




a. R-113 


b. FC-72 


Fig. E.2 Sonic velocities of condensing liquids (R-113 and FC-72) 


Table E.1 Comparison of measured and reported sonic velocities 


R-113 

FC-72 

Temp. 

Meyer 

Calib. 

Diff. 

Temp. 

3M 

Calib. 

Diff. 

°C 

m/s 

m/s 

% 

°C 

m/s 

m/s 

% 

20.1 

714.5 

767.1 

6.9 

10 

555.0 

548.02 

1.3 

14.9 

732.4 

784.8 

6.7 

20 

527.5 

521.5 

1.2 

10.2 

748.2 

801.0 

6.6 

30 

498.0 

501.0 

0.6 
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APPENDIX F 


COMPARISON BETWEEN EQ. 4-6 AND EQ.4-7 

* 

% 


In Chapter 4, two equations (Eqs. 4-6 and 4-7) are introduced to calculated the 
condensate film thickness. Eq. 4-6 is expressed as 


= 


2c m (l +|8) 2 (1 +aTJ 3 p^ 
2+j3 


where 

r a AT 

0 ' Twf; 

Eq. 4-7 is expressed as 

^2 = C aPl *f 


The ratio of the values of 8 given from these two different equations is 

6, _ 2(l+ft) 2 (l+aT t ) 3 Po 
5 2 ' (2+18)?; 

Assuming that the equation for the mean liquid density (p\) is 

^ = Pod + « T ref) 

where 

x T. + T„ t 

ref 2 

Substitution of Eq. F-2 into F-l yields 


(F-l) 


(F-2) 


F.l 
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5, _ 2(1 +_ 0 ? 1+aT 

5 2 2+0 ' 1 +ocT re / 


(F-3) 


Sample Calculations 

(1) condensing liquid: R-113; T, = 7 °C; T ut = 47 °C; (AT = 40 °C and T ref 
°C); p, = 1619.48(1-1.45 x 10 3 (p 0 = 1619.48 and e*=-1.45 x 10' 3 ) 

0 = -^-^xlQ- 3 x(47-7) _ _ 0 05 g 
1 -1.45 x 10' 3 x7 


= 27 


ft = 2(l-0.058) 2 1-1.45X10- 3 X7 3 _ Q 
5 2 2-0.058 1-1.45X10- 3 X27 ; 




-I = |1 - --- 

6, 1 1 0.999 


= ± 0.1 % 


(2) condensing liquid: FC-72; T, = 27 °C; T„, = 57 °C; (AT = 30 °C and T ref = 
42 °C); p, = 1737.86(1-1.26 x 10' 3 (p 0 = 1737.86 and a=-1.26 x 10’ 3 ) 

o _ -1.26xl0- 3 x(57-27) = _ Q Q39 
1-I.26xl0' 3 x27 


«i _ 2(l-0.039) 2 , 1-1.26 x10* 3 x 27 )3 _ n 00Q 
5 2 2-0.039 1-1.26 x10- 3 x42 


h = |1 - = |1 - _L_ 

S 2 5, 1 1 0.999 


= ± 0.1 % 


The above two examples show that the difference between the value of film thickness 
evaluated from Eq. 4-6 and that from Eq. 4-7 is less than ±0.1 %. 
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Table G.l Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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Table G.l Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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Table G.1 Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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Note: 


Based on the data of the average heat transfer coefficients provided in 
Table G.l, the sectional heat transfer coefficients can be 
evaluated using the following equations. 


H ti = H, 


(G-l) 


n t2 = 2H, - n, (G-2) 

Ti# = 3 n 3 - 2 (G-3) 
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Table G.2 Film thickness 
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* x is the distance between the leading edge of the condensing surface to the point where 8 is measured, x is in mm. 


Table G.2 Film thickness (continued) 
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FC-72 1.09 9.38 0.160 0.202 0.216 0.220 0.226 0.244 0.128 


BIBLIOGRAPHY 


3M, (1989), Fluorinert™ Electronic Liquids. 

Abis, L. V., (1969), Forced Convection Condensation inside Horizontal Tubes , 

Ph.D. Dissertation, Kansas State University. 

Altman, M., Norris, R. H., and Staub, F. W., (1960), Local Heat Transfer and 
Pressure Drop for Refrigerants Evaporating in Horizontal Tubes, Journal of Heat 
Transfer, August, 189-198 

Altman, M., Staub, F. W., and Norris, R. H., (1960), Local Heat Transfer and 
Pressure Drop for Refrigerant-22 Condensing inside Horizontal Tubes, Chem. Eng. 
Prog. Symp. Series, 56, 151-159. 

Akers, W. W., Akers, Dean, H. A., and Crasser, O. K., (1958), Condensing Heat 
Transfer within Horizontal Tubes, Chem. Engg. Prog., 54, 89-90. 

Akers, W. W. and Rosson, H. E., (1960), Condensation inside a Horizontal Tube, 
Chem. Eng. Prog. Symp. Series, Heat Transfer, Storrs, 56, 145-149. 

Andritsos, N. and Hanratty, T. J., (1987), Interfacial Instabilities for Horizontal Gas- 
Liquid Flows in Pipelines, Int. J. Multiphase Flow, 13, 583-603 

ASHRAE Handbook of Fundamentals, (1985), SI Edition, American Society of 
Heating Refrigerating and Air-Conditioning Engineers, Inc., Atlanta, GA. 

Azer, N. Z., Abis, L. V., and Swearingen, T. B., (1971), Local Heat Transfer 
Coefficients during Forced Convection Condensation inside Horizontal Tubes, 

Trans. ASHRAE, 77, 182-201. 

Bae, S., (1969), Refrigerant Forced-condensation inside Horizontal Tubes, Ph.D. 
Thesis, MIT. 

Bae, S., Maulbetsh, J. S. and Rohsenow, W. M., (1971), Refrigerant Forced- 
convection Condensation inside Horizontal Tubes, ASHRAE Trans, Part two, 104- 
116. 

Bell, K. J., Taborek, J., and Fenoglio, F., (1970), Interpretation of Horizontal In- 
tube Condensation Heat Transfer Correlations with a Two-Phase Flow Regime 
Map, Chem. Eng. Prog. Symp. Series, 66, 150-163. 


H.l 


Bibliography 


Barry, J. J. and Corradini, M. L., (1988), Film Condensation in the Presence of 
Interfacial Waves, ASME Proceedings, Vol. 2, 1988 National Heat Transfer 
Conference, Houston, TX. 

Carpenter, F. G. and Colburn, A. C., (1951), The Effect of Vapor Velocity on 
Condensation inside Tubes, Proceedings of the General Discussion of Heat 
Transfer, The Institute of Mechanical Engineers and the ASME, July, 20-26. 

Cavallini, A. and Zecchin, R., (1974), A Dimensionless Correlation for Heat Transfer 
in Forced Convection Condensation, 5 th Int. Heat Transfer Conference, Tokyo, 3, 
309-313. 

Chaddock, J. B., (1957), Film Condensation of Vapor in a Horizontal Tube, 
Refrigerating Engineering, 65, 36-41 and 90-95. 

Chandrasekhar, S., (1961), Hydrodynamic and Hydromagnetic Stability, Oxford 
University Press. 

Chang, J. S., Ichikawa, Y., and Irons, G. A., (1982), Flow Regime Characterization 
and Liquid Film Thickness Measurement in Horizontal Gas-Liquid Two-Phase 
Flow by an Ultrasonic Method, AIAA/ASME, June, 7-12. 

Chato, J. C., (1962), Laminar Condensation inside Horizontal and Inclined Tubes , 
Ashrae Journal, 4, 52-60. 

Chen, C. J., (1962), Condensing Heat Transfer in a Horizontal Tube, M.S. Thesis, 
Kansas State University. 

Christodoulou, C., (1987), Condensation in a Horizontal Rectangular Duct, An 
Experimental Study, M.S. Thesis. 

Dittus, F. W. and Boelter, L. M. K., (1930), University of California, Berkeley, 
Publications on Engineering, 2, 443. 

Dukler, A. E., (1977), The Role of Waves in Two Phase Flow: Some New 
Understandings, Chem. Eng. Education, Summer, 108-117. 

Eckert, E. R. G., (1950) Introduction to the Transfer of Heat and Mass, p. 115, 
McGraw-Hill. 

Feldman, S., (1957), On the Hydrodynamic Stability of Two Viscous Incompressible 


H.2 


Bibliography 


Fluids in Parallel Uniform Shearing Motion, J. Fluid Mech., 2, 343-370. 

Jacob, M., (1949), Heat Transfer, Vol. 1, John Wiley, New York. 

Goodykoontz, J. H. and Dorsch., R. G., (1966), Local Heat Transfer Coefficients for 
Condensation of Steam in Vertical Downflow within a 5/8 Inch-diameter Tube, 
NASA-TN-D-3326. 

Goodykoontz, J. H. and Brown, W. F., (1967a), Local Heat-Transfer and Pressure 
Distributions for Freon-113 Condensing in Downward Flow in a Vertical Tube, 
NASA TN D-3952. 

Goodykoontz, J. H. and Dorsch, R. G. (1967b), Local Heat-transfer Coefficients and 
Static Pressures for Condensation of High-velocity Steam within a Tube, NASA 
TND-3953. 

Hagiwara, Y., Esmaeilzadeh, E., Tsutsui, H., and Suzuki, K., (1989), Simultaneous 
Measurement of Liquid Film Thickness, Wall Shear Stress and Gas Flow 
Turbulence of Horizontal Wavy Two- phase Flow, Int. J. Multiphase Flow, 15, 
421-431. 

Hanratty, T. J., (1982), Interfacial Instabilities Caused by Air Flow over a Thin 
Liquid Layer, in Waves on Fluid Interface, Proc. Symp. Math. Rec. Center, The 
University of Wisconsin, Madison, Ed. Richard E. Meyer, Academic Press. 

Hirshburg, R. I. and Florschuetz, L. W., (1982), Laminar Wavy-film Flow: Part II, 
Condensation and Evaporation, Journal of Heat Transfer, 104, 459-464. 

Jensen, R. J. and Yuen, M. C., (1982), Local Heat and Mass Transfer Correlation in 
Horizontal Stratified Concurrent Flow, Proc. 7 th Int. Heat Transfer Conference, 
Munichen, 95-100. 

Kao, T. W. and Park, C., (1972), Experimental Investigation of the Stability of 
Channel Flows, Part 2, Two-layered Co-current Flow in a Rectangular Channel, J. 
Fluid Mech., 52, Part 3, 505. 

Khabenskiy, V. B., Granovskiy, V. S., and Morozov, P. A., (1981), Condensation in 
Stratified Two- Phase Flow inside Horizontal or Sloping Tubes, Heat Transfer- 
Soviet Research, 13, 85-97. 

Kline, S. J. and McClinrock, F. A., (1953), The Description of Uncertainties in 


H.3 


Bibliography 


Single Sample Experiments, Mech. Eng., 75, 3-8. 

Kocamustafaogullari, G., (1985), Two-fluid Modeling in Analyzing the Interfacial 
Stability of Liquid Film Flows, Int. J. Multiphase Flow, 11, 63. 

Koh, J. C. Y., (1962), Film condensation in a Forced-Convection Boundary-Layer 
Flow, Int. J. Heat Mass Transfer, 5, 941-954. 

Kosky, P. G. and Jaster, H., (1976), Condensation Heat Transfer in a Mixed Flow 
Regime, Int. J. Heat Mass Transfer, 19, 95-99. 

Kutsuna, J., I., K. and Nakanishi, S., (1987), Filmwise Condensation of Vapor 
Containing Non-condensable Gas in a Horizontal Duct, Heat Transfer Science and 
Technology, Ed., Wang Bu-Xuan, Hemisphere Publishing Corporation, 
International Symposium on Heat Transfer, Bejing. 

Linehan, J. H., Petrick, M., and El-Wakil, M. M., (1969), On the Interface Shear 
Stress in Annular Flow Condensation, Journal of Heat Transfer, August, 450-452. 

Lockhart, R. W. and Martinelli, R. C. (1949), Proposed Correlation of Data for 
Isothermal Two-phase Two-component Flow in Pipes, Chemical Engineering 
Progress, 45, 39. 

Luu, M. , (1980), Augmentation of Horizontal In-Tube condensation ofR-113, Ph.D. 
Dissertation, Department of Mechanical Engineering, Iowa State University, Ames, 
Iowa. 

Milne-Thomson, L. M., Theoretical Hydrodynamics, MacMillan, (1960). 

Miner, C. S. and Dalton, N. N., (1953), Glycerol, Reinhold Publishing Corporation. 

Moalem, D. and Sideman, S., (1976), Theoretical Analysis of A Horizontal 
Condenser-Evaporator Tube, Int. J. Heat Mass Transfer, 19, 259-270. 

Narain, A.(1991), Condensing Vapor Flows and the Case of Flow between Parallel 
Plates, Recent Advances in Mechanics of Structured Continua, ASME, AMD-Vol. 
117. 

Narain, A. and Kizilyalli, Y., (1991), Pressure Driven Flow of Pure Vapor 
Undergoing Laminar Film Condensation between Parallel Plates, International 
Journal of Non-linear Mechanics, Vol. 26, 501-520. 


H.4 


Bibliography 


Narain, A. and Kamath, R., (1991), Internal Flows of Vapor Undergoing Film 
Condensation, Proceedings of the Twenty-second Midwestern Mechanical 
Conference, Development in Mechanics . Vol. 16, University of Missouri-Rolla, 
411-412. 

Nusselt, W., (1916), Die Obetfldchenkondesation des Wasserdampfes, Z. Ver. Dt. 
Ing. 60, 541-546. 

Reid, R. C. and Sherwood, T. K., (1958), The Properties of Gases and Liquids, 
McGraw-Hill. 

Rohsenow, W. M., Film Condensation, (1970), Applied Mechanics Review, 23, 487- 
495. 

Rohsenow, W. M., Webber, J. H., and Ling, A. T., (1956), Effect of Vapor Velocity 
on Laminar and Turbulent-film Condensation, Transactions of the ASME, 
November, 1637-1643. 

Rosson, H. E. and Myer, J. A., (1965), Point Values of Condensing Film Coefficient 
inside a Horizontal Pipe , Chem Engg. Prog. Symp. Series, 61, 190-199. 

Rufer, C. E. and Kezios, S. P., (1966), Analysis of Two-phase, One-component 
Stratified Flow with Condensation, Journal of Heat Transfer, 88, 265-275. 

Sardesai, R. G., Owen, R. G., and Pulling, D. J., (1981), Flow Regimes for 

Condensation of a Vapor inside a Horizontal Tube, Chemical Engineering Science, 
36, 1173-1180. 

Schlichting, H. (1979), Boundary Layer Theory, Seventh Edition, McGraw-Hill. 

Shah, M. M., (1979), A General Correlation for Heat Transfer during Film 
Condensation inside Pipes, Int. J. Heat Mass Transfer, 22, 547-556. 

Shekriladze, I. G. and Gomelauri, V. I., (1966), Theoretical Study of Laminar Film 
Condensation of Flowing Vapor, Int. J. Heat Mass Transfer, 9, 581-591. 

Soliman, M. and Berenson, P. J., (1970), Flow Stability and Gravitational Effects in 
Condenser Tubes, Proc. 4 th Int. Heat Transfer Conf., Vol. 6, Paris. 

Soliman, M., Schuster, J. R., and Berenson, P. J., (1968), A General Heat Transfer 
Correlation for Annular Flow Condensation, Journal of Heat Transfer, 90, 267- 


Bibliography 


276. 

Sparrow, E. M. and Gregg, J. L., (1959), A boundary Layer Treatment of Laminar 
Film Condensation, J. Heat Transfer, 81, 13-18. 

Sparrow, E. M., Minkowycz, W. J., and Saddy, M., (1967), Forced Convection 
Condensation in the Presence of Non-condensables and Interfacial Resistance, Int. 
J. Heat Mass Transfer, 10, 1829-1845. 

Suryanarayana, N. V. and Merte, H., Jr., (1972), Film Boiling on Vertical Surface, 

J. of Heat Transfer, 94, 377-384. 

Traviss, D. P., Rohsenow, and Baron, A. B., (1973), Forced-convection 
Condensation inside Tubes: A heat Transfer Equation for Condenser Design, 
ASHRAE Trans. Vol 79, Part 1, 157-165. 

Uehara, H., Nakaoka, T., and Kurosaki, M., (1984), Turbulent Film Condensation in 
a Forced-convection Boundary Layer Flow, Heat Transfer-Japanese Research, 13, 
89-104. 

Yih., C. S., (1967), Instability due to Viscosity Stratification, J. Fluid Mech. 27, Part 
2, 337. 


H.6 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing Instructions, searching existing data sources 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of thte 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports. 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington. DC 20503. 

1. AGENCY USE ONLY (Leave blank) 

2. REPORT DATE 

January 1993 

3. REPORT TYPE AND DATES COVERED 

Final Contractor Report 

4. TITLE AND SUBTITLE 

Film Condensation in a Horizontal Rectangular Duct 

5. FUNDING NUMBERS 

WU-506-41-51 
NAG3-71 1 

6. AUTHOR(S) 

Qing Lu and N.V. Suryanarayana 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Michigan Technological University 
College of Engineering 
Houghton, Michigan 49931 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E-7550 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA CR-191061 


11. SUPPLEMENTARY NOTES 

Project Manager, John B. McQuillen, Microgravity Fluids Branch, NASA Lewis Research Center, (216) 433-2876. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Condensation heat transfer in a horizontal rectangular duct was experimentally and analytically investigated. To prevent the dripping of condensate on the Film, 
the experiment was conducted inside a horizontal rectangular duct with vapor condensing only on the bottom cooled plate of the duct. R- 1 13 and FC-72 
(Fluorinert Electronic Fluid developed by the 3M Company) were used as the condensing fluids. The experimental program included measurements of film 
thickness, local and average heat transfer coefficients, wave length, wave speed, and a study of wave initiation. The measured film thickness was used lo obtain 
the local heat transfer coefficient. The wave initiation was studied both with condensation and with an adiabatic air-liquid flow. The test sections used in both 
experiments were identical. Experimental results showed that the average heal transfer coefficient increased with increasing inlet vapor velocity. There was a 
significant increase in the heat transfer after the appearance of interfacial waves (when the inlet Reynolds number was approximately greater than 1 ,000, (XX)). 
The local heat transfer coefficient decreased with axial distance of the condensing surface. The local heat transfer coefficient decreased rapidly with axial 
distance near the leading edge of the condensing surface but was nearly constant towards the trailing edge. It was observed that the condensate flow along the 
condensing surface experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy flow. The change in the flow pattern depended on 
the vapor velocity and the difference between the saturation temperature of the vapor and the condensing surface temperature. The wave length decreased with 
axial distance and the inlet vapor velocity, while the wave speed increased with vapor velocity. An analytical model simulating the condensation process was 
formulated by employing the universal conservation laws. The resulting equations were solved numerically. The heat transfer coefficients predicted from the 
model arc within ±20% of the measured values. Using the present experimental data and analysis, correlations for the average heat transfer coefficients in the 
annular flow regime were developed. The average deviation between the predictions and the experimental values is within ±18%. Also, predictions from some 
correlations selected from the literature for the heat transfer coefficients of condensation in the annular flow regime were compared with the experimental data, it 
is found that the correlations of Shah (1979) and Soliman ct al. (1968) yielded satisfactory predictions. (The average is within ±21 % for Shah's correlation and 
±22 for Soliman’s correlation.) The effects of air velocity, liquid flow rate, and the liquid viscosity on the initiation of interfacial waves were studied with an 
adiabatic air-liquid flow. The liquid viscosity was varied by varying the mass fraction of glycerine in the water-glycerine mixture. It was observed that the inter- 
facial waves were initiated closer to the leading edge when the air velocity or the liquid flow rate was increased. It was found that with condensation, the flow 
appeared more stable than the adiabatic liquid (low; condensation appears to have a damping effect on the initiation of interfacial waves. Several criteria for the 
instability in two-phase flow arc examined and compared with each other. It is found that none of them can be confidently used in the present experiment. With 
the modification of one of the investigated criteria and the use of the present experimental data of the air-liquid flow, a criterion for predicting the instability of 
the two phase flow is tentatively given. This criterion provides a reasonable prediction for the initiation waves with air-liquid flow, as well as with condensation. 


14. SUBJECT TERMS 

Two phase flow; Condensation; Film condensation; Normal gravity 


15. NUMBER OF PAGES 

241 


16 . 


PRICE CODE 

All 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298 102 


FOURTH CLASS MAIL 


National Aeronautics and 
Space Administration 

Lewis Research Center 

Cleveland, Ohio 44135 

Official Business 

Penalty I of Private Use >300 


ADDRESS CORRECTION REQUESTED 




Postage and F ees Paid 
National Aeronautics and 
Space Administration 
NASA 45t 




NASA 


I 


k 


t 

l 


